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Abstract: Different from the system in classical mathematics, a Smarandache system is 
a contradictory system in which an axiom behaves in at least two different ways within the 
same system, i.e., validated and invalided, or only invalided but in multiple distinct ways. 
Such systems exist extensively in the world, particularly, in our daily life. In this paper, we 
discuss such a kind of Smarandache system, i.e., non-solvable ordinary differential equation 
systems by a combinatorial approach, classify these systems and characterize their behaviors, 
particularly, the global stability, such as those of sum-stability and prod-stability of such 
linear and non-linear differential equations. Some applications of such systems to other 
sciences, such as those of globally controlling of infectious diseases, establishing dynamical 
equations of instable structure, particularly, the n-body problem and understanding global 


stability of matters with multilateral properties can be also found. 


Key Words: Global stability, non-solvable ordinary differential equation, general solution, 
G-solution, sum-stability, prod-stability, asymptotic behavior, Smarandache system, inherit 


graph, instable structure, dynamical equation, multilateral matter. 


AMS(2010): 05C15, 34430, 34A34, 37C75, 70F10, 92B05 


§1. Introduction 


Finding the exact solution of an equation system is a main but a difficult objective unless some 
special cases in classical mathematics. Contrary to this fact, what is about the non-solvable 
case for an equation system? In fact, such an equation system is nothing but a contradictory 
system, and characterized only by having no solution as a conclusion. But our world is overlap 
and hybrid. The number of non-solvable equations is much more than that of the solvable 
and such equation systems can be also applied for characterizing the behavior of things, which 
reflect the real appearances of things by that their complexity in our world. It should be noted 
that such non-solvable linear algebraic equation systems have been characterized recently by 
the author in the reference [7]. The main purpose of this paper is to characterize the behavior 


of such non-solvable ordinary differential equation systems. 


1Received November 16, 2012. Accepted March 1, 2013. 
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Assume m, n > 1 to be integers in this paper. Let 
X = F(X) (DES") 
be an autonomous differential equation with F :R” — R” and F(0) = 0, particularly, let 
X=AX (LDES') 
be a linear differential equation system and 
a + aj¢?D 4..-4a,7 =0 (LDE”) 


a linear differential equation of order n with 


G11 G12 +: Ain r1(t) filt, X) 
sue - ¢ t,X 

Peale ee bie | 2) ons F(t,X) = falt,X) ; 
an aAn2 aa Ann In (t) Flt xX) 


where all a;, aij, 1 < 2,7 <n are real numbers with 
7 . : - \T 
X = (%1,42,-++ ,¥n) 


and f;(¢) is a continuous function on an interval [a,b] for integers 0 < i < n. The following 
result is well-known for the solutions of (LDES') and (LDE”) in references. 


Theorem 1.1([13]) Jf F(X) is continuous in 
U(Xo): |t-—tol <a, ||X — Xol| <b (a>0, b> 0) 


then there exists a solution X(t) of differential equation (DES') in the interval |t — to| < h, 
where h = min{a,b/M}, M = ax EGA) 
0 


m 
(t,X)EU (to, X 
Theorem 1.2([13]) Let A; be the k;-fold zero of the characteristic equation 
det(A — AInxn) = |A- AInxn| = 0 
or the characteristic equation 


A” + ay AP + e+ + an1A + an = 0 


with ky +ko+-+--+k, =n. Then the general solution of (LDES') is 


> ci, (t)e™", 
i=1 
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where, c; is a constant, 3;,(t) is an n-dimensional vector consisting of polynomials in t deter- 
mined as follows 















































ti 
Bye =| ™ 
tri 
tit + tie 
Bat) = | ON ™ 
tnit + tne 
ett + batt? +o + tay 
Bi, (e) = | etm + eat? ++ tae 
Tt? a pitt th? Ee 
| Ci (ky 41) ] 
Bets |e 
tn(k1 +1) 
tuit+tig 
S205) ere 
tnit + tna 
- i gag Mt et 
By (t) = Beckett phe} 4 ambi) ps2 + 2. + ton 








En(n=ks +1) 4 


ks—1 1 tn(n—ks+2) 
(eo) a t 


(ks—2)! hace +++ +tran 


with each ty; a real number for 1 <i,7 <n such that det((tij],,.,) A 0; 


Mt, tf l<i<kh; 
A2, tf ky +1 <i< ky; 
a= 





Ag; if ky ko t-+++kgyp +1 <i 
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The general solution of linear differential equation (LDE") is 


s 
YS fea + cia? eee oP Ci(k;—1)t + Cik, er", 
i=1 


with constants qj, 1<i<s,l<j<k. 


Such a vector family 3;(t)e%*, 1 <i <n of the differential equation system (LDES') and 
a family t'e*!, 1 <1 < kj,1 <i <-s of the linear differential equation (L DE") are called the 
solution basis, denoted by 





B={G,(the* |1<i<n} or GH={ tle’ |1<i<s,1<i< kh}. 


We only consider autonomous differential systems in this paper. Theorem 1.2 implies that 
any linear differential equation system (LDES") of first order and any differential equation 
(LDE”) of order n with real coefficients are solvable. Thus a linear differential equation system 
of first order is non-solvable only if the number of equations is more than that of variables, and 
a differential equation system of order n > 2 is non-solvable only if the number of equations 
is more than 2. Generally, such a contradictory system, i.e., a Smarandache system [4]-[6] is 
defined following. 


Definition 1.3((4]-[6]) A rule R in a mathematical system (U;R) is said to be Smarandachely 
denied if it behaves in at least two different ways within the same set %, 1.e., validated and 
invalided, or only invalided but in multiple distinct ways. 

A Smarandache system (;R) is a mathematical system which has at least one Smaran- 
dachely denied rule R. 


Generally, let (51;R1) (H2;Re), ++: ,(2mjRm) be mathematical systems, where R; is a 
rule on 4; for integers 1 <i < m. If for two integers i,j, 1 < i,j <m, 4; A 4; or Ly = Uj; but 





Ri # Rj, then they are said to be different, otherwise, identical. We also know the conception 
of Smarandache multi-space defined following. 


Definition 1.4([4]-[6]) Let (X1;Ri), (L2; Ra), +--+, (Um; Rm) be m > 2 mathematical spaces, 


~ m ed m 
different two by two. A Smarandache multi-space 4 is a union (J X; with rules R = U Ri on 
i=1 


t=1 
S i.e., the rule Ry on %; for integers 1 <1<m, denoted by (E:R). 


A Smarandache multi-space (5; R) inherits a combinatorial structure, i.e., a vertex-edge 


labeled graph defined following. 


Definition 1.5((4]-[6]) Let (5:2) be a Smarandache multi-space with S = US; and R = 
i=1 


U R;. Its underlying graph G =, R| is a labeled simple graph defined by 
i=1 


V (G [=. R]) = {X1, Me,::- , Em}, 


B(G|5,]) = { (D:,2,) | DD; 40,1< 4,5 <m} 
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with an edge labeling 


lead : (Sy, U,) ECE (c [S. R|) — Ra) =W (2 =:) 5 
where w is a characteristic on X; (|X; such that U;()X; is isomorphic to Xp () Xi if and only 
if m(Xi (2X5) = w (Lef)X1) for integers 1 <i, j, k,l <m. 


Now for integers m, n > 1, let 
X =F,(X), X = F(X),---,X = Fy(X) (DES? ) 


be a differential equation system with continuous F; : R"™ — R” such that F;(0) = 0, particu- 
larly, let 
X = A X,-+-,X = ApX,---, X = AmX (LDES}) 


be a linear ordinary differential equation system of first order and 











o™ 4 ler) 4...4a2=0 
n 0 n— o 
al + allel Ny..-+a@le=0 (LDE”) 





g®™ pl all g(™—-2) + aie + ala = 0 


a linear differential equation system of order n with 


k k k 

a di. af tt) 
k k k 

ea or ere 
ol al os all wall 


[k] 


where each a;;" is a real number for integers 0 << k <m, 1<i,j <n. 





Definition 1.6 An ordinary differential equation system (DES},) or (LDES},) (or (LDE™)) 
are called non-solvable if there are no function X(t) (or x(t)) hold with (DES},) or (LDES),) 
(or (LDE®.)) unless the constants. 


The main purpose of this paper is to find contradictory ordinary differential equation 
systems, characterize the non-solvable spaces of such differential equation systems. For such 
objective, we are needed to extend the conception of solution of linear differential equations in 


classical mathematics following. 


Definition 1.7 Let S° be the solution basis of the ith equation in (DES},). The V-solvable, A- 
solvable and non-solvable spaces of differential equation system (DES) are respectively defined 
by 


Ls? (8? and (Js? 1) 8°, 
t=1 t=1 t=1 t=1 


where S° is the solution space of the ith equation in (DES}). 
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According to Theorem 1.2, the general solution of the ith differential equation in (LDES} ) 
or the ith differential equation system in (ZLDE?) is a linear space spanned by the elements 
in the solution basis A; or G; for integers 1 < i < m. Thus we can simplify the vertex-edge 
labeled graph G [> R| replaced each }/; by the solution basis 4; (or @) and >7;()>0, by 
Bi(\B; (or GG) if BiB; ZO (or GG FO) for integers 1 < i,j < m. Such a vertex- 
edge labeled graph is called the basis graph of (LDES},) ((LDE™,)), denoted respectively by 
G[LDES}] or G[LDE®] and the underlying graph of G[LDES},] or G[LDE™], i.e., cleared 
away all labels on G[LDES1] or G[LDE®] are denoted by G[LDES}] or G[LDE®). 

Notice that a i U S°, ie., the non-solvable space is empty only if m = 1 in 


i=1 w=1 
(LDEq). Thus G[LDES'| ~ Kk, or G[LDE"| ~ Kk, only if m = 1. But in general, the 


basis graph G[LDES},] or G[LDE®] is not trivial. For example, let m = 4 and 4? = 
{err erat ersth Bs = fersh, oat, erst). BS = {errt erst erst} and BY = {erst erst eret} 
where \;, 1 < i < 6 are real numbers different two by two. Then its edge-labeled graph 
G([LDES}] or G[LDE™] is shown in Fig.1.1. 


> {erst} 
@ @) 









{errt, ersth {erat ersth 











2 eo oe 
@) @) 


Fig.1.1 





If some functions F;(X), 1 < i < m are non-linear in (DES},), we can linearize these 
non-linear equations X = F;(X) at the point 0, i.e., if 


F(X) = F/(0)X + R(X), 


where F/(0) is an n x n matrix, we replace the ith equation X = F;(X) by a linear differential 


equation 


X= F/O)X 


in (DES},). Whence, we get a uniquely linear differential equation system (LDES},) from 
(DES}.) and its basis graph G[LDES!]. Such a basis graph G[LDES}] of linearized differen- 
tial equation system (DES}_ ) is defined to be the linearized basis graph of (DES},) and denoted 
by G[DES}]. 

All of these notions will contribute to the characterizing of non-solvable differential equation 
systems. For terminologies and notations not mentioned here, we follow the [13] for differential 
equations, [2] for linear algebra, [3]-[6], [11]-[12] for graphs and Smarandache systems, and [1], 
[12] for mechanics. 
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§2. Non-Solvable Linear Ordinary Differential Equations 


2.1 Characteristics of Non-Solvable Linear Ordinary Differential Equations 


First, we know the following conclusion for non-solvable linear differential equation systems 
(LDES}) or (LDE®). 


Theorem 2.1 The differential equation system (LDES},) is solvable if and only if 
(|A1 — AInxn; |A2 — MInxn|s+++ 5 |Am — MInxnl) #1 
i.e., (LDEgq) is non-solvable if and only if 
(|Ai — AInxn,|A2 — AInxn|,+++ 5 [Am — Alnxn|) = 1. 
Similarly, the differential equation system (LDE?.) is solvable if and only if 
(Pi(A), Pa(a),-++ » Pm(A)) #1, 
i.e., (DDE?) is non-solvable if and only if 
(Pi(A), PaQa),-++ » Pm(A)) = 1, 


where P;(X) = A" + all n—1 at dora | 


i(n—1)> + al’! for integers 1<i<m. 


Proof Let Ai1, Ai2,-++ , Ain be the n solutions of equation |A; — AInxn| = 0 and &; the 
solution basis of ith differential equation in (LDES},) or (LDE™,) for integers 1 < i < m. 
Clearly, if (LDES},) ((LDE™)) is solvable, then 


() Bi #9, Le., (Pris diz. saint BO 
t=1 t=1 


by Definition 1.5 and Theorem 1.2. Choose Ao € ia) {\i1,Ai2,°°* Ain}. Then (A — Xo) is a 


common divisor of these polynomials |A1 — AInxn, Ao — Mnxn|)-++ ;|Am — AInxn|. Thus 
(\Ai — AInxn; |A2 — AInxnls-++ y|Am — AInxn|) 4 1. 
Conversely, if 
(|A1 — AInxns|A2 — AInxnl,+++ »|Am — nxn) #1, 


let (A—Ao1), (A—Ao2), +++ , (A—Aor) be all the common divisors of polynomials | A, —AInxn,|A2— 
Mnxnl,-++ ;|Am — MInxn|, where Aoi A Aoz if i # 7 for 1 < i,j <1. Then it is clear that 


Cye*? + Cye*? +--+. + Cher 














is a solution of (LEDq) ((LDE®.)) for constants Cy,C2,--- ,C). 


For discussing the non-solvable space of a linear differential equation system (LEDS1,) or 
(LDE”.) in details, we introduce the following conception. 
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Definition 2.2 For two integers 1< i,j <m, the differential equations 








dX; 
= A,X 
dt 
(LDES}.) 
dX; 
—i — A,X 
dt a 
in (LDES})) or 
o™ + alem—d +...4 ae =0 (EDER) 
2 4 alla 4... 4 ale =0 - 














in (LDE”.) are parallel if Bi(\B; =9. 
Then, the following conclusion is clear. 


Theorem 2.3 For two integers 1 <i,j < m, two differential equations (LDES};) (or (LDE}.)) 
are parallel if and only if 


(|Ai] — nxn |Agl — nxn) = 1 (or (Pi(A), PQA) = 1), 


where (f(x), g(x)) is the least common divisor of f(x) and g(a), P(A) =A" + al’) AME ep 


Gadgs + all fork =1,9. 


Proof By definition, two differential equations (LEDSj,) in (LDESj/,) are parallel if and 
only if the characteristic equations 


|A; —Alnxn| =O and |Aj; —AInxn| =0 


have no same roots. Thus the polynomials |A;| — AInxn and |A;| — AInxn are coprime, which 
means that 


(|Ai — AInxn, |Aj — Alnxn) = 1. 


Similarly, two differential equations (LED7,) in (LDE7,) are parallel if and only if the 
characteristic equations P;(\) = 0 and P;(A) = 0 have no same roots, i.e., (P;(A), Pj(A)) = 1. 














Let f(x) =agz™ +4,07 1 +--+ +Gm12+4m, g(x) = box” + bya" 1 +--+ + bya t by 
with roots 71,%2,°++,%m and y1, y2,°+: , Yn, respectively. A resultant R(f,g) of f(x) and g(x) 
is defined by 


Rf, g) = ag°bo | |e — ys)- 
a,j 


The following result is well-known in polynomial algebra. 


Theorem 2.4 Let f(x) = agxz™ + a,2™ 1 +--+ +@m-1t+am, g(x) = box” + bya" 1 Hee + 
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by—12 + by with roots 71,%2,°°: ,Lm and y1,Y2,°°* ,Yn, respectively. Define a matrix 
ao ay eae Am 0 ee 0 0 
0 ao ay eae Am 0 eee 0 
0 0 0 ao ay am 
Vif,g) = 
ba Bio 3% By OC ote 0 
OG) be: Bg 965: bg. OP ae 
0 0 0 bo by bn 
Then 


R(f, 9) = detV(f, g)- 


We get the following result immediately by Theorem 2.3. 


Corollary 2.5 (1) For two integers 1 < i,j < m, two differential equations (LDES;;) are 
parallel in (LDES},) if and only if 


R(|Ai — Mnxn|,|Aj — nxn) #9; 
particularly, the homogenous equations 
V(\Ai — AInxnl, [Aj — Mnxnl)X =0 
have only solution (0,0,---,0)7 if |Aj — AInxn| = oA” + aA"! +++) + Gn-1A + an and 
a 


|A; — Alnxn| = BoA” + bart +++ + dniA + bn. 
(2) For two integers 1 < i,j < m, two differential equations (LDE};) are parallel in 
(LDE”) if and only if 
R(Pi(A), Py(A)) 0, 
particularly, the homogenous equations V(P;(A), P;(A))X =0 have only solution (0,0,-++ ,0)7. 
eae 


Proof Clearly, |A; — Inxn| and |A; — AInxn| have no same roots if and only if 
R(|Ai — AInxn|,|Aj — Inxnl) #9, 


which implies that the two differential equations (LEDS/,) are parallel in (LEDS,,,) and the 
homogenous equations 
WA Aiea ihe 


have only solution (0,0,--- ,0)7. That is the conclusion (1). The proof for the conclusion (2) 
— 


2n 











is similar. 
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Applying Corollary 2.5, we can determine that an edge (4, ;) does not exist in G[LDES, | 
or G[LDE*,] if and only if the ith differential equation is parallel with the jth differential equa- 
tion in (LDES?,) or (LDE”,). This fact enables one to know the following result on linear 


non-solvable differential equation systems. 


Corollary 2.6 A linear differential equation system (LDES},) or (LDE®,) is non-solvable if 
G(LDES},) & Km or G(LDE®,) & Km for integers m,n > 1. 


2.2. A Combinatorial Classification of Linear Differential Equations 


There is a natural relation between linear differential equations and basis graphs shown in the 


following result. 


Theorem 2.7 Every linear homogeneous differential equation system (LDES},) (or (LDE™)) 
uniquely determines a basis graph G[LDES},] (G[LDE™]) inherited in (LDES}) (or in(LDE”)). 
Conversely, every basis graph G uniquely determines a homogeneous differential equation system 
(LDES}) (or (LDE®)) such that GILDES},] ~ G (or G[LDE®] ~ G). 


Proof By Definition 1.4, every linear homogeneous differential equation system (LDES} ) 
or (LDE”.) inherits a basis graph G[LDES}] or G[LDE™], which is uniquely determined by 
(LDES1) or (LDE®). 

Now let G be a basis graph. For Vu € V(G), let the basis 4, at the vertex v be By, = 
{ B,(t)e* | 1<i<n,} with 


Ai, tf 1<Stsky; 


ro, if ky +1<i< ky; 
a= 





We construct a linear homogeneous differential equation (LDES"') associated at the vertex v. 
By Theorem 1.2, we know the matrix 


tir ota2 tes tan, 
T= ta. toa +++ tan, 
tnyl tny2 mene tnyny 


is non-degenerate. For an integer 7, 1 <i < s, let 


Global Stability of Non-Solvable Ordinary Differential Equations With Applications 11 


be a Jordan black of k; x k; and 


Then we are easily know the solution basis of the linear differential equation system 


= = AX (LDES') 


with X = [21 (t), r2(t),--- ,an,(t)]7 is nothing but 4, by Theorem 1.2. Notice that the Jordan 
black and the matrix T are uniquely determined by #,. Thus the linear homogeneous differen- 
tial equation (LDES") is uniquely determined by &,. It should be noted that this construction 
can be processed on each vertex v € V(G). We finally get a linear homogeneous differential 
equation system (LDES1},), which is uniquely determined by the basis graph G. 

Similarly, we construct the linear homogeneous differential equation system (ZLDE” ) for 
the basis graph G. In fact, for Vu € V(G), let the basis Z,, at the vertex u be B, = { t'e** | 1 < 
i<s,1<Il<k;}. Notice that A; should be a k;-fold zero of the characteristic equation P(A) = 0 
with kj + kg +--- +k, =n. Thus P(\;) = P’(\y) =--- = P@ YA) = 0 but PMH) (A) A 0 
for integers 1 <i < s. Define a polynomial P,,(A) following 





s 


P(A) = ]Ja-A* 


i=1 


associated with the vertex u. Let its expansion be 
Py (A) = A” + Gut"! H+ + Gu(n—1)A + Gun- 
Now we construct a linear homogeneous differential equation 
Praag ee Sr Wee, cima ee Ciera Qu(n—1)2 + dunt = 0 (LP DE?) 


associated with the vertex u. Then by Theorem 1.2 we know that the basis solution of (LDE”) 
is just G,. Notices that such a linear homogeneous differential equation (LDE”) is uniquely 


constructed. Processing this construction for every vertex u € V(G), we get a linear homoge- 














neous differential equation system (LDE”,). This completes the proof. 
Example 2.8 Let (LDE*) be the following linear homogeneous differential equation system 


Z— 34+ 2x = 0 1 
z—5x%+ 6x = 0 2 
&—7z+12% =0 3 
x — 9x + 20x = 0 
Z—11e+ 302 = 0 5 
xz—7x+6x = 0 


y 





Ne SA Re SE 





aD 
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2 


d d 
where ¢ = <* and ¢ = —. Then the solution basis of equations (1) — (6) are respectively 


fete), fe eft), {ert ce}, fet eF4}, fe5, e61 {e% ce} and its basis graph is shown in 
Fig.2.1. 
{e"} {e**} 
{cSt et} {c3t, At} 


{e*, erry {er} Lee”, et} 








Fig.2.1 The basis graph H 


Theorem 2.7 enables one to extend the conception of solution of linear differential equation 
to the following. 


Definition 2.9 A basis graph GILDES},] (or G[LDE®]) is called the graph solution of the 
linear homogeneous differential equation system (LDES},) (or (LDE™)), abbreviated to G- 


solution. 


The following result is an immediately conclusion of Theorem 3.1 by definition. 


Theorem 2.10 Every linear homogeneous differential equation system (LDES},) (or (LDE™)) 
has a unique G-solution, and for every basis graph H, there is a unique linear homogeneous 
differential equation system (LDES},) (or (LDE®,)) with G-solution H. 


Theorem 2.10 implies that one can classifies the linear homogeneous differential equation 
systems by those of basis graphs. 


Definition 2.11 Let (LDES},), (LDES},)’ (or (LDE™), (LDE”)') be two linear homo- 
geneous differential equation systems with G-solutions H, H’'. They are called combinato- 
rially equivalent if there is an isomorphism py : H — H’, thus there is an isomorphism 
yp: H — H' of graph and labelings 0, 7 on H and H' respectively such that p(x) = Ty(x) for 
Va € V(H)J E(H), denoted by (LDES},) & (LDES1)! (or (LDE”) © (LDE®)’). 


{e ve ay {e aby Le 2t e ey 
e 3} 


jek, et 


ea} 


{e Ste BE te PEL de ate PEt 


















Fig.2.2. The basis graph H’ 
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Example 2.12 Let (LDE?,)' be the following linear homogeneous differential equation system 











#+3¢+22=0 1) 
&+5é¢ +62 =0 2) 
£+7é+12c =0 3) 
#494 +202 =0 4) 
#+1lé+302=0 (5) 
£+7¢ + 6c =0 6) 


Then its basis graph is shown in Fig.2.2. 
Let y: H — H’ be determined by y({e>**, e%#*}) = {e~**, es") and 


p({er*, ert} (Here ")) = fee leo ae Ph 


for integers 1 < i,k <6 and j =i+1 = 6(mod6), 1 =k +1 = 6(mod6). Then it is clear that 
H & H'. Thus (LDE”)' is combinatorially equivalent to the linear homogeneous differential 
equation system (LDE?.) appeared in Example 2.8. 


Definition 2.13 Let G be a simple graph. A vertex-edge labeled graph 6: G— Z* is called 
integral if 0(uv) < min{@(u), O(v)} for Vuv € E(G), denoted by G’°. 

Let Gie and ens be two integral labeled graphs. They are called identical if Gy z G2 and 
O(x) = r(y(x)) for any graph isomorphism y and Vx € V(G1) UE(G1), denoted by Gi? = GP. 


For example, these labeled graphs shown in Fig.2.3 are all integral on K4—e, but Gis = ere 
GrAtGe 
V # G3. 


3 2 4 4 BD 3 3 1 3 
1 2 2 1 2 2 
4 2 3 3 4 4 1 4 
Ge Gi Gy 
Fig.2.3 


Let G[LDES}] (G[LDE®]) be a basis graph of the linear homogeneous differential equa- 
tion system (LDES},) (or (LDE™)) labeled each v € V(G[LDES}]) (or v € V(G[LDE™])) 
by &,. We are easily get a vertex-edge labeled graph by relabeling v € V(G[LDES}]) (or 
v € V(G[LDE™])) by |B,| and uw € E(G[LDES}]) (or uv € E(G[LDE™])) by |@uf\ Fol. 
Obviously, such a vertex-edge labeled graph is integral, and denoted by G/[LDES}] (or 
G'|LDE®™]). The following result completely characterizes combinatorially equivalent linear 
homogeneous differential equation systems. 


Theorem 2.14 Let (LDES},), (LDES},)' (or (LDE®,), (LDE”)') be two linear homogeneous 
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differential equation systems with integral labeled graphs H, H’. Then (LDES},) x (LDES.) 
(or (LDE”) © (LDE®,)') if and only if H = H’. 


Proof Clearly, H = H' if (LDES1,) © (LDES},) (or (LDE") © (LDE®)') by defini- 
tion. We prove the converse, i.e., if H = H’ then there must be (LDES?,) © (LDES},)’ (or 
(LDE”) & (LDE®)’). 

Notice that there is an objection between two finite sets 5), So if and only if |S,| = |So|. 
Let 7 be a1 —1 mapping from ¥, on basis graph G[LDES})] (or basis graph G[LDE™]) to 
By on basis graph G[LDES}]' (or basis graph G[LDE™]') for v,v' € V(H’). Now if H = H’, 
we can easily extend the identical isomorphism zdy on graph H to a 1— 1 mapping id}, : 
G[LDES}] — G[LDES}]' (or id}, : G[LDE™] — G[LDE®]’) with labelings 0: v — @, and 
0, : vw — By on G[LDES}], GILDES}])’ (or basis graphs G[LDE™], G[LDE®]'). Then 
it is an immediately to check that id3,0(x) = 6’r(x) for Vz € V(G[LDES},]) E(G[LDES,,]) 
(or for Vz € V(G[LDE),]) UU E(G[LDE®,])). Thus id}, is an isomorphism between basis graphs 
G[LDES}] and G[LDES1J (or G[LDE®] and G[LDE®)). Thus (LDES1) ‘“ (LDES1,)' 
(or (LDE”) eo (LDE”)'). This completes the proof. 














According to Theorem 2.14, all linear homogeneous differential equation systems (LDES} ) 
or (LDE?.) can be classified by G-solutions into the following classes: 


Class 1. G[LDES!] ~ Km or G[LDE®] ~ Ky for integers m,n > 1. 


The G-solutions of differential equation systems are labeled by solution bases on K,, and 
any two linear differential equations in (LDES},) or (LDE”.,) are parallel, which characterizes 
m isolated systems in this class. 


For example, the following differential equation system 


i+ 3x + 2x =0 
«—5x+ 6x =0 
@+ 2% —- 3x =0 





is of Class 1. 
Class 2. G[LDES1] ~ Km or GILDE®] ~ Ky, for integers m,n > 1. 


The G-solutions of differential equation systems are labeled by solution bases on complete 
graphs Ky, in this class. By Corollary 2.6, we know that GILDES1,] ~ Km or GILDE”,] ~ Km 
if (LDES),) or (LDE”.) is solvable. In fact, this implies that 


() B= [) fA) 46. 


veV (Km) u,vEV (Km) 


Otherwise, (LDES},) or (LDE”,) is non-solvable. 


For example, the underlying graphs of linear differential equation systems (A) and (B) in 
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the following 


Z-— 32+ 27 =0 
Z-—34+27 =0 
r-x=0 


(A) (B) # — 5+ 6x =0 
# —4¢ +32 =0 
#—4¢ +32 =0 


&+ 2x — 3x =0 








are respectively K4, K3. It is easily to know that (A) is solvable, but (B) is not. 


Class 3. G[LDES!] ~ G or G[LDE”] ~ G with |G| = m but G # Km, Km for integers 
m,n > 1. 


The G-solutions of differential equation systems are labeled by solution bases on G and all 
linear differential equation systems (LDES},) or (LDE®,) are non-solvable in this class, such 
as those shown in Example 2.12. 


2.3 Global Stability of Linear Differential Equations 


The following result on the initial problem of (LDES') and (LDE”) are well-known for differ- 
ential equations. 


Lemma 2.15([13]) Fort € [0,00), there is a unique solution X(t) for the linear homogeneous 
differential equation system 


= AX (L' DES") 


with X(0) = Xo and a unique solution for 
a™ + aya D+...4a,7 =0 (L2 DE) 
with «(0) = ao, 2'(0) = 2,--- ,a("-Y(0) = ce), 


Applying Lemma 2.15, we get easily a conclusion on the G-solution of (LDES},) with 
X,(0) = X® for Vu € V(G) or (LDE) with x(0) = x, 2'(0) = 2),--- ,a@-Y(0) = a") by 
Theorem 2.10 following. 


Theorem 2.16 For t € [0,00), there is a unique G-solution for a linear homogeneous dif- 
ferential equation systems (LDES},) with initial value X,(0) or (LDE™) with initial values 
ay(0),2/,(0),--- , a6" (0) for Vu € V(G). 


For discussing the stability of linear homogeneous differential equations, we introduce the 
conceptions of zero G-solution and equilibrium point of that (LDES},) or (LDE”,) following. 


Definition 2.17 A G-solution of a linear differential equation system (LDES},) with initial 
value X,(0) or (LDEP.) with initial values x,(0), 2/,(0),--- 06" (0) for Vu € V(G) is called 
a zero G-solution if each label B; of G is replaced by (0,--- ,0) (|\B| times) and Bi (|B; by 
(0,--- ,0) (Ai A;| times) for integers 1 < i,j <m. 
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Definition 2.18 Let dX/dt = A,X, oO) + aya") +--+ aynx = 0 be differential equations 
associated with vertex v and H a spanning subgraph of G[LDES},| (or G[LDE™]). A point 
X* € R” is called a H-equilibrium point if A,X* = 0 in (LDES},) with initial value X,(0) 
or (X*)" + yi (X*)""1 + ++ + Gun X* = 0 in (LDE”,) with initial values x,(0), x/,(0),---, 
a0) for Vu € V(A). 


We consider only two kind of stabilities on the zero G-solution of linear homogeneous 


differential equations in this section. One is the sum-stability. Another is the prod-stability. 


2.3.1 Sum-Stability 


Definition 2.19 Let H be a spanning subgraph of GILDES},] or G[LDE™] of the linear 
homogeneous differential equation systems (LDES},) with initial value X,(0) or (LDE”) with 
initial values vy(0),2/,(0),--+, a Yo). Then G[LDES}] or GILDE®] is called sum-stable 
or asymptotically sum-stable on H if for all solutions Y,(t), v € V(H) of the linear differential 
equations of (LDES},) or (LDE™) with |Y,(0)—Xv(0)| < dy exists for allt >0,| SS Y,(t)- 
veV (HA) 
So Xy(t)| <¢, or furthermore, lim| S> Y(t)- SS Xp(t)| =0. 
veV(H) $0 veV(H) veV(H) 
Clearly, an asymptotic sum-stability implies the sum-stability of that G[LDES},] or G[LDE” ]. 
The next result shows the relation of sum-stability with that of classical stability. 


Theorem 2.20 For a G-solution G[LDES}] of (LDES},) with initial value X,,(0) (or GILDE™] 
of (LDEP.) with initial values x,(0),x/,(0),--- , as" (0)), let H be a spanning subgraph of 
G[LDES}] (or G[LDE™]) and X* an equilibrium point on subgraphs H. If G[LDES}] (or 
G[LDE™]) is stable on any Vu € V(H), then GILDES}] (or G[LDE®]) is sum-stable on H. 
Furthermore, if G[LDES}] (or G[LDE®]) is asymptotically sum-stable for at least one vertex 
v€ V(H), then GILDES}] (or G[LDE™}) is asymptotically sum-stable on H. 


Proof Notice that 
| SS pKwO- SS pH) < So plV(t) - Xd) 
a (A) veEV(H) veV(H) 


and 


lim| SY) pY(t)- DY) peXOlS DY) po lim|¥.(t) -— Xv(@)]. 


veV (HA) veV(H) veV(H) 














Then the conclusion on sum-stability follows. 
For linear homogenous differential equations (LDES") (or (LDE")), the following result 
on stability of its solution X(t) =0 (or x(t) = 0) is well-known. 


Lemma 2.21 Let y = max{ Red] |A — Alnxn| = 0}. Then the stability of the trivial solution 
X(t) = 0 of linear homogenous differential equations (LDES') (or x(t) = 0 of (LDE™)) is 


determined as follows: 


(1) ify <0, then it is asymptotically stable; 
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(2) ify > 0, then it is unstable; 
(3) if y = 0, then it is not asymptotically stable, and stable if and only if m'(A) = m(A) 
for every with ReX = 0, where m(A) is the algebraic multiplicity and m’(A) the dimension of 


eigenspace of X. 


By Theorem 2.20 and Lemma 2.21, the following result on the stability of zero G-solution 
of (LDES},) and (LDE®) is obtained. 


Theorem 2.22 A zero G-solution of linear homogenous differential equation systems (LDES} ) 
(or (LDE™.)) is asymptotically sum-stable on a spanning subgraph H of GILDES},] (or G[LDE”]) 
if and only if Rea, < 0 for each B,,(t)e%* € By, in (LDES") or Redy <0 for each t'-e* € G, 
in (LDE®.) hold for Vu € V(#). 


Proof The sufficiency is an immediately conclusion of Theorem 2.20. 
Conversely, if there is a vertex v € V(H) such that Rea, > 0 for B,(t)e%* € Ay in 
(LDES') or Red, > 0 for t'e*! € G, in (LDE®,), then we are easily knowing that 
jim B,(t)e%"* — 00 


if a, > 0 or B,(t) Aconstant, and 


lim t' et 5 00 
t—-00 


if Ay > 0 or ly > 0, which implies that the zero G-solution of linear homogenous differential 














equation systems (LDES') or (LDE") is not asymptotically sum-stable on H. 


The following result of Hurwitz on real number of eigenvalue of a characteristic polynomial 
is useful for determining the asymptotically stability of the zero G-solution of (LDES},) and 
(LDE”). 


Lemma 2.23 Let P(A) =A" +a,A\"-1 +++» +@n-1A + an be a polynomial with real coefficients 
aj, 1<i<nand 


ay 1 0 vee 0 
a3 ag ay, 0 sane 0 
ay, 1 
Ai = |ai|, Ag = pot An = a5 a4 a3 a2 ay 0 tee 0 
a3 «a2 
0 ame Qn 








Then ReX < 0 for all roots X of P(A) if and only if A; > 0 for integers 1<i<n. 


Thus, we get the following result by Theorem 2.22 and lemma 2.23. 


Corollary 2.24 Let Aj, A$,--- ,A? be the associated determinants with characteristic polyno- 
mials determined in Lemma 4.8 for Vv € V(G[LDES}]) or V(G[LDE™]). Then for a spanning 
subgraph H < G[LDES}] or G[LDE®], the zero G-solutions of (LDES},) and (LDE™,) is 
asymptotically sum-stable on H if AY > 0,A3 >0,---,A® > 0 for Vu € V(A). 
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Particularly, if n = 2, we are easily knowing that Red < 0 for all roots A of P(A) if and 
only if a, > 0 and ag > 0 by Lemma 2.23. We get the following result. 


Corollary 2.25 Let H < G[LDES}] or G[LDE®] be a spanning subgraph. If the characteristic 
polynomials are \7 + a¥\ + a3 for v € V(H) in (LDES},) (or (L"DE?,)), then the zero G- 
solutions of (LDES},) and (LDE2,) is asymptotically sum-stable on H if at > 0, a} > 0 for 
Vue V(#). 


2.3.2 Prod-Stability 


Definition 2.26 Let H be a spanning subgraph of G[LDES},] or G[LDE™] of the linear 
homogeneous differential equation systems (LDES},) with initial value X,(0) or (LDE®,) with 
initial values xy(0),2/,(0),---, af Yo). Then GILDES}] or G[LDE™] is called prod-stable 


or asymptotically prod-stable on H if for all solutions Y,(t), v € V(A) of the linear differential 
equations of (LDES},) or (LDE™) with |Y,(0)—X.v(0)| < 5, exists for allt >0,| [T] Y.(t)- 
veV (A) 
II Xv(t)| <, or furthermore, lim| [] Yi(t)- JT] Xv(t)| =0. 
veV(H) t>0  veV(H) veV (H) 
We know the following result on the prod-stability of linear differential equation system 
(LDES1) and (LDE"). 


Theorem 2.27 A zero G-solution of linear homogenous differential equation systems (LDES} ) 
(or (LDE™ )) is asymptotically prod-stable on a spanning subgraph H of GILDES},] (or G[LDE”]) 


if and only if S > Reay <0 for each B,(t)e%* € By in (LDES') or SY > Redry <0 for 
veEV(H) veV(H) 


each tert € ©, in (LDE™). 
Proof Applying Theorem 1.2, we know that a solution X,(t) at the vertex v has the form 


n 


X(t) = D> eB, (te 


i=l 


Whence, 


I] *@) = 


I> 
veV (HA) EV(A) t=1 
_ II, eve" = STL cde] etn, 
i=1 veV(H 


i=1 veV(H) 








Whence, the zero G-solution of homogenous (LDES?,) (or (LDE®.)) is asymptotically sum- 
stable on subgraph H if and only if >> Rea, < 0 for VG,(t)e%* € ZB, in (LDES') or 
veEV (A) 
>> Red, < 0 for Vtlver? € @, in (LDE™). 
veV (HA) 














Applying Theorem 2.22, the following conclusion is a corollary of Theorem 2.27. 


Corollary 2.28 A zero G-solution of linear homogenous differential equation systems (LDES} ) 
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(or (LDE”)) is asymptotically prod-stable if it is asymptotically sum-stable on a spanning 
subgraph H of G[LDES}] (or G[LDE®]). Particularly, it is asymptotically prod-stable if the 


zero solution 0 is stable on Yu € V(H). 


Example 2.29 Let a G-solution of (LDES},) or (LDE”,) be the basis graph shown in Fig.2.4, 
where v, = {e72!, e734, 8*}, ug = {e734 e-4*}, ug = Lew 4t, 7 BL, vg = {e754 8, 84}, 


—2t 
’ 


vs = f{e*,e ©}, ug = {e *,e 7, e *}. Then the zero G-solution is sum-stable on the triangle 


V4U5 U6, but it is not on the triangle v;v2v3. In fact, it is prod-stable on the triangle v,v2v3. 


ve {e**} U1 





{e“} 
U5 fe} {e**} 


{e~% 








v4 {e*} us 


Fig.2.4 A basis graph 


§3. Global Stability of Non-Solvable Non-Linear Differential Equations 


For differential equation system (DES!},), we consider the stability of its zero G-solution of 
linearized differential equation system (LDES},) in this section. 


3.1 Global Stability of Non-Solvable Differential Equations 


Definition 3.1 Let H be a spanning subgraph of G[DES}] of the linearized differential equation 
systems (DES},) with initial value X,(0). A point X* € R” is called a H-equilibrium point of 
differential equation system (DES}) if fy(X*) =0 for Vu € V(A). 


Clearly, 0 is a H-equilibrium point for any spanning subgraph H of G[DES},] by definition. 
Whence, its zero G-solution of linearized differential equation system (LDES1},) is a solution 
of (DES},). 


Definition 3.2 Let H be a spanning subgraph of G[DES},] of the linearized differential equation 
systems (DES?) with initial value X,(0). Then G[DES}] is called sum-stable or asymptoti- 
cally sum-stable on H if for all solutions Y,(t), v € V(H) of (DES}) with ||Y,(0) — Xv(0)|| < 6 
exists for all t > 0, 


+ KwO- YL XO] <a 


veV(H) veV(H) 


or furthermore, 
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lim S> YwWO- SS xX] =0, 


veV(H) veV(H) 
and prod-stable or asymptotically prod-stable on H if for all solutions Y,(t), v € V(H) of 
(DES}.) with ||Y,(0) — X1(0)|| < 6, exists for allt > 0, 


I] “®- TI xl <e 


veV(H) veV(H) 








or furthermore, 


lim|| J] Y@®- J] x] =o. 


t—0 
veEV (A) veEV (A) 














Clearly, the asymptotically sum-stability or prod-stability implies respectively that the 
sum-stability or prod-stability. 


Then we get the following result on the sum-stability and prod-stability of the zero G- 
solution of (DES}). 


Theorem 3.3 For a G-solution G[DES}] of differential equation systems (DES},) with initial 
value X,(0), let Hy, Hz be spanning subgraphs of G[DES}.]. If the zero G-solution of (DES},) 
is sum-stable or asymptotically sum-stable on H, and Hp, then the zero G-solution of (DES},) 
is sum-stable or asymptotically sum-stable on Hy U Ho. 

Similarly, if the zero G-solution of (DES},) is prod-stable or asymptotically prod-stable on 
Hy and X,(t) is bounded for Vv € V(H2), then the zero G-solution of (DES},) is prod-stable 
or asymptotically prod-stable on Hy \) Ho. 


Proof Notice that 
|X1 + X2|| < |[Xa|] + [|X|] and ||X1Xal] < || X1||||-Xell 


in R”. We know that 


S- X,(t)| = S> X(t)+ SS Xd) 


veV (1) U V(H2) veV (1) veV (H2) 





IA 
4 
= 
aS 
ts 
5 








and 


Il X,(t) 


veEV(H1)U V(A2) veV (Ai) veEV(H2) 


lA ll 
cs ee 
Le Le 
| 
= 
i. 
x S 
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Whence, 


>» 


vE€V(A1)U V(A2) 


pee: 


if € =€, + €2 with 
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<e€ or lim 
veEV(Ai)U V (Ao) 








S> X(t] <a: and | S> X,(t)|| < e2 
veV (1) veV (Ae) 
or 
lim S> X,(t)]] =0 and lim S> X,(t)|] = 0. 
veV(A1) veV (Az) 


This is the conclusion (1). For the conclusion (2), notice that 


I 


vEV(A1)U V(A2) 


X(t) 


I] xt 


veV(H1) 


<) [[ ~@] [[ %@]) < me 
veV (1) veV (H2) 
)|| <e and Il X,(t)|| <M. 


vEV (He) 


Consequently, the zero G-solution of (DES},) is prod-stable or asymptotically prod-stable on 


H,U Ap. 


Theorem 3.3 enables one to get 


oO 





the following conclusion which establishes the relation of 


stability of differential equations at vertices with that of sum-stability and prod-stability. 


Corollary 3.4 For a G-solution G[DES},] of differential equation system (DES},) with initial 
value X,(0), let H be a spanning subgraph of G[DES},]. If the zero solution is stable or 
asymptotically stable at each vertex v € V(A), then it is sum-stable, or asymptotically sum- 
stable and if the zero solution is stable or asymptotically stable in a vertex u € V(A) and X(t) 
is bounded for Vu € V(H) \ {u}, then it is prod-stable, or asymptotically prod-stable on H. 


It should be noted that the converse of Theorem 3.3 is not always true. For example, let 


S- X,(t)|| <a+e and S- Xy(t)|| < -at+e. 


veEV(H1) 


Then the zero G-solution G[DES}, 
on subgraphs H; and Ho, but 


B 


veV(A1U Ae) 


X(t) 





Thus the zero G-solution G[DES}, 
subgraphs A; |) H2. Similarly, let 


vEV (He) 


of differential equation system (DES},) is not sum-stable 


< 


S> Xo(t)]] + 


veEV(A1) 


Ss) X(t) 


vEV (Ha) 


=€. 


of differential equation system (DES) is sum-stable on 
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I] x @]) < = and || S> X,()|| <¢” 
veV(A1) veV (Ae) 

for a real number r. Then the zero G-solution G[DES},] of (DES},) is not prod-stable on 

subgraphs H, and X,(t) is not bounded for v € V(H2) if r > 0. However, it is prod-stable on 


subgraphs A () H2 for 


I] xls) Tf ©! IT] x @ll=e 


veEV(AMiU A2) veV(H1) veV (H2) 


3.2 Linearized Differential Equations 


Applying these conclusions on linear differential equation systems in the previous section, we 
can find conditions on F;(X), 1 < i < m for the sum-stability and prod-stability at 0 following. 
For this objective, we need the following useful result. 


Lemma 3.5(({13]) Let X = AX + B(X) be a non-linear differential equation, where A is a 
constant nx n matrix and Red; < 0 for all eigenvalues A; of A and B(X) is continuous defined 
ont >0, ||X|| <a@ with 

IB 


lim “" _ 9, 
|x||+0 |X| 


Then there exist constants c > 0, 8 >0 and 6, 0<6<a such that 


|X| <e< implies that ||X(t)|| < cee~%/?. 
Cc 


Theorem 3.6 Let (DES?) be a non-linear differential equation system, H a spanning subgraph 
of G[DES},] and 


F(X) = Fi (0) X + R,(X) 


such that 
|| Ru (X)]| 


im Xt =o 
x0 |X| 


for tv € V(H). Then the zero G-solution of (DES},) is asymptotically sum-stable or asymp- 
totically prod-stable on H if Rea, <0 for each B,,(t)e%* € By, v € V(H) in (DES). 


Proof Define c = max{c,, v € V(H)}, ¢ = min{e,, v € V(A)} and 6 = min{G,, v € 
V(H)}. Applying Lemma 3.5, we know that for Vu € V(A), 


6 
||Xv(0)|| <e < 5 implies that ||Xy(¢)|| < geese, 
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Whence, 
S> Xo) < SD XO < |Hlece*? 
veV (A) veV(H) 
I] ©) < [TL 1e@]) s clleltle tee2, 
veV (A) veV(H) 
Consequently, 
lim S> X,(t)] +0 and lim I] x@]} -o. 
veV(H) veV(H) 





Thus the zero G-solution (DES?) is asymptotically sum-stable or asymptotically prod-stable 
on H by definition. 














3.3 Liapunov Functions on G-Solutions 


We have know Liapunov functions associated with differential equations. Similarly, we introduce 
Liapunov functions for determining the sum-stability or prod-stability of (DES},) following. 


Definition 3.7 Let (DES!},) be a differential equation system, H < G|DES}] a spanning 
subgraph and a H-equilibrium point X* of (DES},). A differentiable function L: @ + R 
defined on an open subset O C R” is called a Liapunov sum-function on X* for H if 


(1) L(X*) =0 ond =) >O0if YS X(t) #4 X*; 


vEeV( veEV (HA) 
Qil + x) <ofr rT xX Ax*, 
veV(H) veV(H) 
and a Liapunov prod-function on X* for H if 
(1) L(X*) =0 and L ( II x10) >Oif JTL X(t) 4x"; 
H) 


vEeV( veV (HA) 


V(A) veV (A) 


2 E( II (0) <0 fr Tl X(t) #X*. 


Then, the following conclusions on the sum-stable and prod-stable of zero G-solutions of 
differential equations holds. 


Theorem 3.8 For a G-solution G[DES},| of a differential equation system (DES},) with 
initial value X,(0), let H be a spanning subgraph of G[DES},] and X* an equilibrium point of 
(DES!) on H. 


(1) If there is a Liapunov sum-function L : @ — R on X*, then the zero G-solution 
G[DES}] is sum-stable on X* for H. Furthermore, if 


ED Sea a0 


veV (HA) 
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for So X(t) # X*, then the zero G-solution G[DES},] is asymptotically sum-stable on X* 
veEV (HA) 
for H. 


(2) If there is a Liapunov prod-function L: @ > R on X* for H, then the zero G-solution 
G[DES} |] is prod-stable on X* for H. Furthermore, if 


LY JJ X@)] <0 


veV (HA) 


for JT] Xv(t) #4 X*, then the zero G-solution G[DES}] is asymptotically prod-stable on X* 
veV(H) 
for H. 


Proof Let ¢ > 0 be a so small number that the closed ball B.(X*) centered at X* with 
radius «€ lies entirely in @ and w the minimum value of L on the boundary of B,(X%*), i-e., 


the sphere S.(X*). Clearly, @ > 0 by assumption. Define U = {X © B.(X*)|L(X) < a}. 
Notice that X* € U and L is non-increasing on > X,(t) by definition. Whence, there are 


veEV (A) 
no solutions X,(t), v € V(H) starting in U such that S> X,(t) meet the sphere S.(X%*). 
veEV (HA) 
Thus all solutions X,(t), v € V(H) starting in U enable 5> X,(t) included in ball B,(X*). 
veEV (HA) 


Consequently, the zero G-solution G[DES},] is sum-stable on H by definition. 


Now assume that 


L{ So X(t) ] <0 


veV (HA) 
for So X,(t) # X*. Thus L is strictly decreasing on 5S > X,(t). If X(t), ve V(A) are 
veEV(H) veV(H) 
solutions starting in U — X* such that > Xy(tn) — Y* for n — oo with Y* € B.(X*), 
veEV (HA) 


then it must be Y* = X*. Otherwise, since 


LS XsG) | She 


veV (HA) 


LE] S° X(t) } <0 
veV (HA) 


forall So X,(t) 4 X* and 
veEV (HA) 

Xy(tn) | > L(Y*) 

vEV a 


by the assumption 


by the continuity of L, if Y* # X*, let Y,(t),v € V(A) be the solutions starting at Y*. Then 


for any 7 > 0, 
L| >> Yn) ] < L(Y"). 
veEV (A) 
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But then there is a contradiction 


L| So Xoltn+n) | < L(Y") 
veEV (A) 
yields by letting Y,(0) = 54> Xux(tn) for sufficiently large n. Thus, there must be Y,* = X*. 
veV (HA) 
Whence, the zero G-solution G[DES}] is asymptotically sum-stable on H by definition. This 
is the conclusion (1). 














Similarly, we can prove the conclusion (2). 


The following result shows the combination of Liapunov sum-functions or prod-functions. 


Theorem 3.9 For a G-solution G[DES},| of a differential equation system (DES},) with 
initial value X_(0), let Hi, Hz be spanning subgraphs of G[DES},], X* an equilibrium point of 
(DES},) on Hi \) He and 

os a,j" y! 


i>0,j>0 
be a polynomial with a;,; > 0 for integers i,j > 0. Then R* (Ly, L2) is a Liapunov sum-function 
or Liapunov prod-function on X* for H, |) H2 with conventions for integers i,j,k, > 0 that 


ag Te ye X(t) | tamLPL} S- X,(t) 
ve€V (Hi U V(A2) veV (Hi U V(A2) 


=ajLi| >> X,(t)| £4 
veEV(A1) Pr. 
+api Lt ‘> X(t) i, 


veV(H1) he 4 


if L1, Le are Liapunov sum-functions and 


aij Li 13 Il X,(t) | +a LiLy 
vEV(A1U V(A2) 


=ajli| |] X(t) | 24 
veV (Hi) ane 


+anlh} JT X@)|25{ [J x@ 


veV(M) veV (H2) 


if L1,L2 are Liapunov prod-functions on X* for Hy, and Ho, respectively. Particularly, if 
there is a Liapunov sum-function (Liapunov prod-function) L on H, and Ho, then L is also a 
Liapunov sum-function (Liapunov prod-function) on Hy |) Ho. 


Proof Notice that 


di — aij (Gemenre: + pTLA Es) 
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if 74,7 > 1. Whence, 


aij LL S- X,(t) | >0 
veV (Hi U V(A2) 
if 
Ii{ S> X,(t)] 20 and Ig2{ S> X(t] 0 
veV(A1) veV (A) 
and 
pe 
Mois TAts) X,(t) < 0 
veV (Hi U V (Ae) 

if 


Ei{ S> X(t)| <0 and bot S> X,(t)} <0. 
veV(A1) veV (Ae) 


Thus R*(L 1, D2) is a Liapunov sum-function on X* for Hi J Ao. 
Similarly, we can know that Rt(L1, L2) is a Liapunov prod-function on X* for Hy U He if 














Ly, L2 are Liapunov prod-functions on X* for H, and Ao. 
Theorem 3.9 enables one easily to get the stability of the zero G-solutions of (DES}.). 


Corollary 3.10 For a differential equation system (DES), let H < G[DES}]] be a spanning 
subgraph. If Ly is a Liapunov function on vertex v for Vu € V(H), then the functions 


Lg= >> L, and LE = JI Ly 
veV(H) veV(H) 


are respectively Liapunov sum-function and Liapunov prod-function on graph H. Particularly, 
if L = Ly forVu € V(A), then L is both a Liapunov sum-function and a Liapunov prod-function 
on H. 


Example 3.11 Let (DES},) be determined by 


dx, /dt = Aw21 dx, /dt = A121 dx, /dt = Ani £1 
dx /dt => A122 dx /dt => A22X2 dix / dt = An2X2 
dxy,/dt = AinXn diy /dt = r2nXn diy /dt = Xnn&n 


where all yj, 1 <i< m,1 <j < nare real and \yj, A Aij. if j1 A Jo for integers 1 <i <m. 
Let L=2?+a3+---+22. Then 





L= Ney + N23 a rin Xt? 


for integers 1 < i <n. Whence, it is a Liapunov function for the ith differential equation if 
rij < 0 for integers 1 < j < n. Now let H < G[LDES},] be a spanning subgraph of G[LDES}.,]. 
Then L is both a Liapunov sum-function and a Liapunov prod-function on H if A,; < 0 for 
Vu € V(#) by Corollaries 3.10. 
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Theorem 3.12 Let L: @ — R be a differentiable function with L(0) = 0 and L ( ee x) > 
vEV (HA) 


0 always holds in an area of its e-neighborhood U(e) of 0 for ¢ > 0, denoted by U*(0,¢) such 
area of €-neighborhood of 0 with L SS X | >0andH < G[DES}] be a spanning subgraph. 
veEV (HA) 


(1) If 


L S- X|||<™M 
veV (HA) 


with M a positive number and 


i S> xX] >0 


veV (HA) 


in U*(0,€), and for Ve > 0, there exists a positive number ci, ¢2 such that 


L S- X |} >c, >0 implies Lb S- X | >c>O0, 
veV(H) veV(H) 
then the zero G-solution G[DES},] is not sum-stable on H. Such a function L: @ > R is 


called a non-Liapunov sum-function on H. 


(2) If 


LY JJ x] <n 
veEV (A) 
with N a positive number and 
L| [J x]>o 
vEV (HA) 


in U*(0,€), and for Ve > 0, there exists positive numbers d,,d2 such that 


L Il X |} >d,>0 implies L Il X | >d2>0, 
veV(H) veV (A) 
then the zero G-solution G[DES},] is not prod-stable on H. Such a function L: 6 > R is 


called a non-Liapunov prod-function on H. 


Proof Generally, if ||L(X)|| is bounded and L (X) > 0 in U+(0,e), and for Ve > 0, there 
exists positive numbers ci, c2 such that if D(X) > c; > 0, then L (X) > c2 > 0, we prove that 
there exists t; > to such that ||X(t1,to)|| > €o for a number €9 > 0, where X(t1, to) denotes 
the solution of (DES) passing through X (to). Otherwise, there must be ||X (£1, to)|| < €0 for 
t > to. By L(X) > 0 we know that L(X(t)) > L(X(to)) > 0 for t > to. Combining this fact 
with the condition L (X) > co > 0, we get that 

t 


L(X(t)) = L(X(to)) + : 


to 


dL(X(s)) 


a > L(X (to)) + co(t — to). 
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Thus L(X (t)) too if t Foo, a contradiction to the assumption that L(X) is bounded. 
Whence, there exists t; > to such that 


|| X(t1, to) || > €0. 


Applying this conclusion, we immediately know that the zero G-solution G[DES}]] is not sum- 











stable or prod-stable on H by conditions in (1) or (2). 





Similar to Theorem 3.9, we know results for non-Liapunov sum-function or prod-function 
by Theorem 3.12 following. 


Theorem 3.13 For a G-solution G[DES}] of a differential equation system (DES},) with 
initial value X,(0), let Hy, Hz be spanning subgraphs of G[DES},], 0 an equilibrium point of 
(DES!) on Hi He. Then Rt(L1,L2) is a non-Liapunov sum-function or non-Liapunov 
prod-function on 0 for Hi |) Hz with conventions for 


CPE EMBs S- XG) ag Ld S- X,(t) 
veV (M1 U V(A2) veV (M1 U V(A2) 
and 
ai; Li 13 Il X(t) | + awLPL} Il X,(t) 
v€V (M1 U V(A2) veV (M1 U V(A2) 


the same as in Theorem 3.9 if Ly, Lz are non-Liapunov sum-functions or non-Liapunov prod- 
functions on 0 for Hy and Ho, respectively. Particularly, if there is a non-Liapunov sum- 
function (non-Liapunov prod-function) L on Hy and Ho, then L is also a non-Liapunov sum- 


function (non-Liapunov prod-function) on Hy VU Ae. 


Proof Similarly, we can show that R*(Lj, D2) satisfies these conditions on H; J H2 for 
non-Liapunov sum-functions or non-Liapunov prod-functions in Theorem 3.12 if DL, D2 are 
non-Liapunov sum-functions or non-Liapunov prod-functions on 0 for H; and He, respectively. 











Thus Rt (Ly, L2) is a non-Liapunov sum-function or non-Liapunov prod-function on 0. 





Corollary 3.14 For a differential equation system (DES}), let H < G[DES}]] be a spanning 
subgraph. If Ly is a non-Liapunov function on vertex v for Vu € V(H), then the functions 


LE= S° L, and B= JI L, 


veV(H) veV(H) 


are respectively non-Liapunov sum-function and non-Liapunov prod-function on graph H. Par- 
ticularly, if L = Ly for Wu € V(A), then L is both a non-Liapunov sum-function and a non- 
Liapunov prod-function on H. 


Example 3.15 Let (DES},) be 


Ly = a7 a di 3 ep) = Aor? a Aor Ly => Ama? = Ama 


4 M1 5 d2 F m 
2 = Zt12%2 t2 = Zt12%2 LQ = VW F142 


2 2 2 
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with constants \; > 0 for integers 1 < i < m and L(a1,a22) = 2? — 242. Then L(21,22) = 
4); L (21,22) for the i-th equation in (DES},). Calculation shows that L(x1,x%2) > 0 if 
xy > V2x or 41 < —V2x9 and L(21, £2) > 4c? for L(x1,22) > cin the area of L(a1, x2) > 0. 
Applying Theorem 3.12, we know the zero solution of (DES},) is not stable for the i-th equation 
for any integer 1 <i < m. Applying Corollary 3.14, we know that LD is a non-Liapunov sum- 
function and non-Liapunov prod-function on any spanning subgraph H < G[DES}]. 


§4. Global Stability of Shifted Non-Solvable Differential Equations 


The differential equation systems (DES},) discussed in previous sections are all in a same 
Euclidean space R”. We consider the case that they are not in a same space R”, i.e., shifted 
differential equation systems in this section. These differential equation systems and their 
non-solvability are defined in the following. 


Definition 4.1 A shifted differential equation system (SDES},) is such a differential equation 


system 





X, = F,(X1), Xo = Fo(X2),--- Xm = Fin(Xm) (SDES}) 
with 

Xy = (21, 22,°°* , 81, %1041), £142))°** Lin), 

X2q = (&1,£2,°°* , £1, ©2141), La42),°** » Lan), 

Xm = (€1, £2, +++ 5 Bt Lm(141)) Um(l42)0°** + Bmn); 
where %1,%2,°°* , 1, % 4,1 <t<m,1l<j7 <n—Ll are distinct variables and Fs: R” > R” 





is continuous such that F,(0) =0 for integers 1<s<m. 

A shifted differential equation system (SDES},) is non-solvable if there are integers i,j, 1 < 
i,j <m and an integer k, 1 <k <1 such that oll (t) # oll] (t), where oll (t), wlll (t) are solutions 
xp(t) of the i-th and j-th equations in (SDES},), respectively. 


The number dim(SDESj;,) of variables #1, %2,-++ ,@1,Tii4j),1 <i<mil<j<n-lin 
Definition 4.1 is uniquely determined by (SDES1,), i-e., dim(SDES1,) =mn—(m—1)l. For 
classifying and finding the stability of these differential equations, we similarly introduce the 
linearized basis graphs G[SDES1,] of a shifted differential equation system to that of (DES}), 
i.e., a vertex-edge labeled graph with 





V(G[SDES,,]) = {Bill <i < m}, 
E(G[SDES;,]) = {(Bi, Bi )|FiN Fi #91 < i,9 < mj, 


where &; is the solution basis of the i-th linearized differential equation X; = F/(0)X; for 
integers 1 < i < m, called such a vertex-edge labeled graph G[SDES},] the G-solution of 
(SDES},) and its zero G-solution replaced 4; by (0,--- ,0) (|@i| times) and (1) 4B; by 
(0,--- ,0) (|Ai (| B;| times) for integers 1 < i,j <m. 

Let (LDES}), (LDES})' be linearized differential equation systems of shifted differential 
equation systems (SDES}) and (SDES},) with G-solutions H, H’. Similarly, they are called 
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combinatorially equivalent if there is an isomorphism y : H — H’ of graph and labelings 
6, 7 on H and H’ respectively such that y6(~) = Ty(x) for Vz € V(H)U E(A), denoted by 
(SDES!) © (SDES1,)'. Notice that if we remove these superfluous variables from G[SDES}], 
then we get nothing but the same vertex-edge labeled graph of (LDES},) in R!. Thus we can 
classify shifted differential similarly to (LDES},) in R!. The following result can be proved 
similarly to Theorem 2.14. 


Theorem 4.2 Let (LDES},), (LDES,,)' be linearized differential equation systems of two 
shifted differential equation systems (SDES},), (SDES1,)' with integral labeled graphs H, H’. 
Then (SDES},) & (SDES1,)’ if and only if H = H'. 


The stability of these shifted differential equation systems (SDES} ) is also similarly to 
that of (DES},). For example, we know the results on the stability of (SDES},) similar to 
Theorems 2.22, 2.27 and 3.6 following. 


Theorem 4.3 Let (LDES},) be a shifted linear differential equation systems and H < G[LDES}] 
a spanning subgraph. A zero G-solution of (LDES},) is asymptotically sum-stable on H if and 
only if Rea, <0 for each 3,,(t)e%* € B, in (LDES") hold for Vv € V(H) and it is asymptot- 
ically prod-stable on H if and only if > Reay <0 for each B,(t)e*' € B, in (LDES?). 
veV(H) 
Theorem 4.4 Let (SDES},) be a shifted differential equation system, H < G[SDES},] a 
spanning subgraph and 
F(X) = F, (0) X + R,(X) 
such that R(X] - 
xo ||| 
for Vv € V(H). Then the zero G-solution of (SDES}) is asymptotically sum-stable or asymp- 
totically prod-stable on H if Rea, <0 for each 3,,(t)e%* € By, v € V(H) in (SDES}). 


For the Liapunov sum-function or Liapunov prod-function of a shifted differential equation 
system (SDES?,), we choose it to be a differentiable function L : @ C R&™(SPZ5m) —, R with 
conditions in Definition 3.7 hold. Then we know the following result similar to Theorem 3.8. 


Theorem 4.5 For a G-solution G[SDES}]] of a shifted differential equation system (SDES} ) 
with initial value X,(0), let H be a spanning subgraph of G[DES},] and X* an equilibrium 
point of (SDES}) on H. 


(1) If there is a Liapunov sum-function L: @ Cc Riim(SDES;,) _, R on X*, then the zero 
G-solution G[SDES}] is sum-stable on X* for H, and furthermore, if 


Ll Se MG <0 


veV (HA) 


for SS X(t) # X*, then the zero G-solution G[SDES},| is asymptotically sum-stable on 
veEV (A) 
X* for H. 
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(2) If there is a Liapunov prod-function L : @ Cc Rim(SDES;,)  R on X* for H, then 
the zero G-solution G[SDES},| is prod-stable on X* for H, and furthermore, if 


LY [J x.@] <0 
veEV (HA) 
for |] Xv(t) 4 X*, then the zero G-solution G[SDES}] is asymptotically prod-stable on 


veV (HA) 
X* for H. 


§5. Applications 


5.1 Global Control of Infectious Diseases 


An immediate application of non-solvable differential equations is the globally control of infec- 
tious diseases with more than one infectious virus in an area. Assume that there are three kind 
groups in persons at time ¢, i.e., infected I(t), susceptible S(t) and recovered R(t), and the 
total population is constant in that area. We consider two cases of virus for infectious diseases: 


Case 1 There are m known virus \,, 42,--+,Vm with infected rate k;, heal rate h; for integers 


1<i<m and an person infected a virus ¥; will never infects other viruses ¥; for 7 #4. 


Case 2 There are m varying \\,%2,--- ,¥%m from a virus V with infected rate k;, heal rate hy 


for integers 1<1<m such as those shown in Fig.5.1. 


Q-O----@ 


Fig.5.1 


We are easily to establish a non-solvable differential model for the spread of infectious 
viruses by applying the SIR model of one infectious disease following: 


S = —kSI S = —k2SI S = —kpSI 
I =k SI — hil T=kSI-hI °-:: I = kmSI — mI (DES}) 
R=hl R=hgl R=hmlI 


Notice that the total population is constant by assumption, i.e., S++ R is constant. 
Thus we only need to consider the following simplified system 


S=-k SI S=-kSI S=-k,,SI 


; ; (DES) 
f=hSI-hyI T= kySI — hel T= kmSI — mI 


The equilibrium points of this system are I = 0, the S-axis with linearization at equilibrium 
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points 


S=-kS S=—-kS S = —kmS 
; : ‘ ve (LDES},) 
T=khS—h, T=kS — he T=knS—hm 


Calculation shows that the eigenvalues of the 7th equation are 0 and k;S—h;, which is negative, 
i.e., stable if 0 < S' < hi/k; for integers 1 <i < m. For any spanning subgraph H < G[LDES}], 
we know that its zero G-solution is asymptotically sum-stable on H if 0 < S < hy/ky for 
v € V(H) by Theorem 2.22, and it is asymptotically sum-stable on H if 


Ss (&S-hm)<0 ie, OxS< Yom / Yb 


veV (A) veEV (A) veV(H) 
by Theorem 2.27. Notice that if [;(t), S;(t) are probability functions for infectious viruses 
¥,, 1 <i < m in an area, then [] Z(t) and [] 5i(t) are just the probability functions for 


t=1 ji 
all these infectious viruses. This fact enables one to get the conclusion following for globally 
control of infectious diseases. 


Conclusion 5.1 Form infectious viruses V,,%2,--- ,¥%m in an area with infected rate k;, heal 
rate hy for integers 1 <i<m, then they decline to 0 finally if 


soe 5h 5 
1=1 t=1 


i.e., these infectious viruses are globally controlled. Particularly, they are globally controlled if 
each of them is controlled in this area. 


5.2 Dynamical Equations of Instable Structure 


There are two kind of engineering structures, i.e., stable and instable. An engineering structure 
is instable if its state moving further away and the equilibrium is upset after being moved 
slightly. For example, the structure (a) is engineering stable but (b) is not shown in Fig.5.2, 


A C3 











Fig.5.2 


where each edge is a rigid body and each vertex denotes a hinged connection. The motion of 
a stable structure can be characterized similarly as a rigid body. But such a way can not be 
applied for instable structures for their internal deformations such as those shown in Fig.5.3. 
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Fig.5.3 


Furthermore, let A, Yo,--- , Am be m particles in R? with some relations, for instance, 
the gravitation between particles Y; and #; for 1 < i, 7 < m. Thus we get an instable structure 
underlying a graph G with 

V(G) = {Fi P2,-** eae i 

E(G) = {(A, Y;)|there exists a relation between Y; and F;}. 

For example, the underlying graph in Fig.5.4 is Cy. Assume the dynamical behavior of particle 


F,, at time t has been completely characterized by the differential equations X = F,(X,t), 
where X = (21, 22,23). Then we get a non-solvable differential equation system 


underlying the graph G. Particularly, if all differential equations are autonomous, i.e., depend 
on X alone, not on time t, we get a non-solvable autonomous differential equation system 


X =F,(X), 1<i<m. 


All of these differential equation systems particularly answer a question presented in [3] for 
establishing the graph dynamics, and if they satisfy conditions in Theorems 2.22, 2.27 or 3.6, 
then they are sum-stable or prod-stable. For example, let the motion equations of 4 members 
in Fig.5.3 be respectively 


AB: Xap = 0; CD: Xcp = 0, AC: Xac = aac, BC: Xgc =agc, 


where X4p,Xcp,Xac and Xgc denote central positions of members AB,CD, AC, BC and 
aAc,@pc are constants. Solving these equations enable one to get 


Xap =capt+dap, Xac =aact? +cact + dac, 


Xcp =ccpt+dop, Xgc =azct® + cect t+ dzc, 


where c4p,Cac,ccp,cBc,¢d4B,d4ac,dcp,dgc are constants. Thus we get a non-solvable dif- 
ferential equation system 


xX 0; xX 0, xX GAC; X =asgc, 
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or a non-solvable algebraic equation system 


X =capt+dap, X =a,ct? + cact+dac, 





X=copt+dcp, X= apct? +cpct+ dsc 


for characterizing the behavior of the instable structure in Fig.5.3 if constants c4B, CAC, CCD; CBC; 
dap, dac,dcop,dgc are different. 
Now let X1, X2,--- , Xm be the respectively positions in R? with initial values X?, X9,--- , X° 


9 m? 
X9, X9,---, X° and Mj, Mo,--- , Mm the masses of particles P:, P2,--- , Pm. If m = 2, then 
from Newton’s law of gravitation we get that 


Xo —X +. X,—- X 
2 1 x, =GM, 1 2 


X, = GM, ———,,, wo? 
|X2— X1| |X1 — X9| 


where G is the gravitational constant. Let X = X2— X1 = (#1, 22,23). Calculation shows that 


“ xX 
X = —G(M;+4+ M2) 3° 
|X| 


Such an equation can be completely solved by introducing the spherical polar coordinates 
1 =rcos@cosé 


t= rcos@cosé 


z3 =rsind 


with r > 0,0 < 6 < 7,0 < 0 < 27, where r = ||X||, 6 = Z2X0z,0 = ZX'ox with X' 
the projection of X in the plane xoy are parameters with r = a/(1 + cos@) hold for some 





constants a,¢. Whence, 


x) Gam f ( rae it and X(t) =-ca, f ( | rae ‘is 


Notice the additivity of gravitation between particles. The gravitational action of particles 
Py, P2,--+, Am on FY can be regarded as the respective actions of Wy, Po,---, Am on FY, 
such as those shown in Fig.5.4. 


Fig.5.4 
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Thus we can establish the differential equations two by two, i.e., A, acts on Y, We acts on 


PF, +--+, Am acts on FY and get a non-solvable differential equation system 
i ae 
X =GM;——_,, Bx, 1<i<m. 
[Xi — X| 


Fortunately, each of these differential equations in this system can be solved likewise that of 
m = 2. Not loss of generality, assume X;(t) to be the solution of the differential equation in 
the case of A; 4 YP, 1<i<m. Then 


a X;—X 
X(t)= So X()=E YO wf (fF 5a) dt 
PAP PAP [Xi — X| 

is nothing but the position of particle A at time t in R? under the actions of A, 4 FY for 
integers 1 < i < m, i.e., its position can be characterized completely by the additivity of 


gravitational force. 
5.3 Global Stability of Multilateral Matters 


Usually, one determines the behavior of a matter by observing its appearances revealed before 
one’s eyes. If a matter emerges more lateralities before one’s eyes, for instance the different 
states of a multiple state matter. We have to establish different models, particularly, differential 
equations for understanding that matter. In fact, each of these differential equations can be 
solved but they are contradictory altogether, i.e., non-solvable in common meaning. Such a 
multilateral matter is globally stable if these differential equations are sum or prod-stable in all. 

Concretely, let $,52,---,5m be m lateral appearances of a matter .W@ in R® which are 
respectively characterized by differential equations 


where X; € R°, a 3-dimensional vector of surveying parameters for S;,1 <%i< m. Thus we get 


a non-solvable differential equations 


in R3. Noticing that all these equations characterize a same matter .@, there must be equilib- 
rium points X* for all these equations. Let 


H,(X,t) = Hi(X*)X + R,(X*), 


where 
aC 
nn ah 


is an n X n matrix. Consider the non-solvable linear differential equation system 


X =Hi(X*)X, 1<i<m (LDES} ) 
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with a basis graph G. According to Theorem 3.6, if 


ARGON 
hex xn 


for integers 1 < 7 < m, then the G-solution of these differential equations is asymptotically 
sum-stable or asymptotically prod-stable on G if each Real! < 0 for all eigenvalues all of 
matrix H/(X*), 1<i<m. Thus we therefore determine the behavior of matter @ is globally 
stable nearly enough X*. Otherwise, if there exists such an equation which is not stable at the 
point X*, then the matter WH is not globally stable. By such a way, if we can determine these 
differential equations are stable in everywhere, then we can finally conclude that M is globally 
stable. 

Conversely, let @ be a globally stable matter characterized by a non-solvable differential 
equation 

X = H,(X,t) 


for its laterality S;, 1<%7<m. Then the differential equations 


are sum-stable or prod-stable in all by definition. Consequently, we get a sum-stable or prod- 
stable non-solvable differential equation system. 


Combining all of these previous discussions, we get an interesting conclusion following. 


Conclusion 5.2 Let MS A be respectively the sets of globally stable multilateral matters, 
non-stable multilateral matters characterized by non-solvable differential equation systems and 
QE,DVE the sets of sum or prod-stable non-solvable differential equation systems, not sum or 
prod-stable non-solvable differential equation systems. then 

(1) VM €. A° > A(DES}) € GE; 

(2) VM 6M? > IDES.) ¢ GE. 

Particularly, let @ be a multiple state matter. If all of its states are stable, then @ is 
globally stable. Otherwise, it is unstable. 
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Abstract: In this paper, we introduce a m‘"-root Randers changed Finsler metric as 


L(x, y) = L(a,y) + B(x, y), 


where L = {dizig---im (x)y"ty? Lyi ban is a m*"-root metric and (-is one form. Further 


we obtained the relation between the v- and hy- curvature tensor of m‘”-root Finsler space 


h 


and its m*"-root Randers changed Finsler space and obtained some theorems for its $3 and 


S4-likeness of Finsler spaces and when this changed Finsler space will be Berwald space 


h 


resp. Landsberg space). Also we obtain T-tensor for the m*”-root Randers changed Finsler 
g g 


space F™. 


Key Words: Randers change, m‘”-root metric, Berwald space, Landsberg space, $3 and 
S4-like Finsler space. 
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§1. Introduction 


Let F” = (M”,L) be a n-dimensional Finsler space, whose M” is the n-dimensional differen- 
tiable manifold and L(x, y) is the Finsler fundamental function. In general, L(x, y) is a function 
of point « = (a*) and element of support y = (y’), and positively homogeneous of degree one 
in y. In the year 1971 Matsumoto [6] introduced the transformations of Finsler metric given by 


L'(a,y) = L(a,y) + B(a,y) 
L'?(x,y) = L(x, y) + B?(2,y), 


where, 3 = b;(x)y’ is a one-form [1] and b;(x) are components of covariant vector which is a 
function of position alone. If L(z,y) is a Riemannian metric, then the Finsler space with a 
metric L(x, y) = a(x, y) + G(x, y) is known as Randers space which is introduced by G.Randers 
[5]. In papers [3, 7, 8, 9], Randers spaces have been studied from a geometrical viewpoint and 
various theorem were obtained. In 1978, Numata [10] introduced another (-change of Finsler 
metric given by L(x,y) = u(2,y) + B(a,y) where pw = {ai;(y)y’y9}2 is a Minkowski metric 
and ( as above. This metric is of similar form of Randers one, but there are different tensor 
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properties, because the Riemannian space with the metric a is characterized by Ch, = 0 and 
on the other hand the locally Minkowski space with the metric w by Rnijk = 0, Crash = 90. 

In the year 1979, Shimada [4] introduced the concept of m‘” root metric and developed it as 
an interesting example of Finsler metrics, immediately following M.Matsumoto and S.Numatas 
theory of cubic metrics [2]. By introducing the regularity of the metric various fundamental 
quantities as a Finsler metric could be found. In particular, the Cartan connection of a Finsler 
space with m-th root metric could be discussed from the theoretical standpoint. In 1992-1993, 
the m-th root metrics have begun to be applied to theoretical physics [11, 12], but the results 
of the investigations are not yet ready for acceding to the demands of various applications. 

In the present paper we introduce a m*’-root Randers changed Finsler metric as 


L(z,y) = L(z,y) a B(x, y) 


’_root metric. This metric is of the similar form to 


where L = {ai;i9.-i,, ()yy2°Y" }m isa mt 
the Randers one in the sense that the Riemannian metric is replaced with the m*”-root metric, 
due to this we call this change as m‘"-root Randers change of the Finsler metric. Further we 
obtained the relation between the v-and hv-curvature tensor of m*’-root Finsler space and its 
m*-root Randers changed Finsler space and obtained some theorems for its $3 and S4-likeness 
of Finsler spaces and when this changed Finsler space will be Berwald space (resp. Landsberg 


space). Also we obtain T-tensor for the m*”-root Randers changed Finsler space F™. 


§2. The Fundamental Tensors of F” 


We consider an n-dimensional Finsler space F” with a metric L(a,y) given by 


L(a,y) = L(a,y) + di(x)y’ (1) 
where 
in, i ak 
L = {Gizig--im (L)Y Y? soy laa (2) 
By putting 
(I). Lago y) = Gigs sig, (x)y"? y"* ae yim (3) 
(II). L"~?ai;(x,y) = Gijigige--ig, (ZY? y"* ++ y'™ 


(IIL). LP ass (2, Y) = Gizhigig--im (ayy »--y'™ 
Now differentiating equation (1) with respect to y’, we get the normalized supporting element 
1; = O,L as 


=a, +); (4) 


’root metric. Again differen- 


where a; = 1; is the normalized supporting element for the m* 
tiating above equation with respect to y?, the angular metric tensor his a L0;0;L is given 


as 


=> (5) 
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where h;; is the angular metric tensor of m‘"-root Finsler space with metric L given by [4] 


hag = (m— 1) (aig — a5) (6) 
eo eae 
The fundamental metric tensor 9;; = Vj = hj; +11; of Finsler space F” are obtained from 
equations (4), (5) and (6), which is given by 








Gij = (m 1)Tai; t {1 (m 1)T}aja; t (aib; t a;b;) t bib; (7) 
L 
where T = i It is easy to show that 
dr = Aa Tei + bit oe (m= D(aig = aiay) oe (m — 2) (aigk ~ Gian) 


L , L L 


Therefore from (7), it follows (h)hv-torsion tensor Cijz = On of the Cartan’s connection CT 
are given by 


2LCijx = (m—1)(m— 2)raiye + KL—(m—Ur}(m— 1)l(ayjae (8) 


+A; hai Aki a; ) } (m 1) (aij bx t ajKd; t anid; ) 








(m — 1)(ajajby + ajanb; + ajand;) + (m— 1){(2m — 1)7 — 38} a;aj;ax 
In view of equation (6) the equation (8) may be written as 
Cijk = TCigk + Wagimicy Rautmierr Tass) sisi orm fixate) (9) 


where m; = b; — Pai and C;;, is the (h)hv-torsion tensor of the Cartan’s connection CT of the 


m*"-root Finsler metric L given by 
2LC 5k = (m = 1)(m = 2){ ijk = (Giz Ok + GjRa4 + Aid; ) + 2a;a; Ax} (10) 


Let us suppose that the intrinsic metric tensor a;;(z,y) of the m‘’-root metric L has non- 
vanishing determinant. Then the inverse matrix (a’’) of (a;;) exists. Therefore the reciprocal 
metric tensor g’? of F” is obtain from equation (7) which is given by 
‘ aw ne 
j# = 1 wee b? + (m— 1)r mee = (a'b? + a?b*) (11) 
(m —1)r (m—1)r(1+Q) (m—1)r(1 +4) 

where a’ = aaj, b' = ab; b? =5'b;, .g =a"b; =.0;b' = B/L. 
Proposition 2.1 The normalized supporting element l;, angular metric tensor hij, metric 
tensor gi; and (h)hv-torsion tensor Ci, of Finsler space with m"_root Randers changed metric 
are given by (4), (5), (7) and (9) respectively. 


§3. The v-Curvature Tensor of F” 


From (6), (10) and definition of m; and a‘, we get the following identities 


aa; = i. aijha" = Qjk, Cijpa' = 0, hija’ = 0, (12) 


mia’ = 0, hi jb! = 3Mi, mb" — (b? — q’) 


m*-Root Randers Change of a Finsler Metric 41 


To find the v-curvature tensor of F”, we first find (h)hv-torsion tensor CL = =% PC oie 


Ai 1 i l i i a 
a’ 1 1 ; 
——~{m;mp + — hr} — Cir” 
Fara" * Go on—a)™~ Gaara” 

where LO%, = LCjrpa"" = (m — 1){a%, — (Sjax + da; + atajn) + 2a*ajax}, 

hi = hjpa'™ = (m— 1) (6; —a’‘aj;) (14) 

m' =m,a"™ =b'—qa', and a, = a’ ajrk 
From (12) and (14), we have the following identities 

Cijrhy = Cih pr = =(m _ as Cigr™m™ = Crd", (15) 


mh; =(m—1)m, mim = (b — q), 


From (9) and (13), we get after applying the identities (15) 


2 1 
CijrChr = any Crh + OL —(CijnmMr + Cigk™Mn + Chjrmi + Chik™;) (16) 
1 1 
———(Cj;,h Chrkhi;)0" + —————~ (b? — q?)hizh 
* Shia rn + Carkhij) + Tm DEL! T )hijhne 


1 
aT (2hijmnme + 2henmim,; + hjnmume 


+hj~pmMiMp + hinmjmr + hirmjmn) 


Now we shall find the v-curvature tensor Shiik = Cijr Ct, _ Cine Of 5. The tensor is obtained 


from (16) and given by 





Shigk = Spt 7 Ty Cir Che “ hij™Mnk + Ankmi; } (17) 
= ay + O(n {higmnk + harm; } 
where 2 2) | 
1 -—q m—1)-_4 
Maj m= 1k gro + AL ae 5) mam} (18) 


and the symbol ©(;,){--- } denotes the exchange of j,k and subtraction. 


Proposition 3.1 The v-curvature tensor Shigk of m'"_root Randers changed Finsler space Er 


with respect to Cartan’s connection CT is of the form (17). 


It is well known [13] that the v-curvature tensor of any three-dimensional Finsler space is 

of the form 
L? Snijk = S(haghir — hai) (19) 
Owing to this fact M. Matsumoto [13] defined the S3-like Finsler space F” (n > 3) as such a 
Finsler space in which v-curvature tensor is of the form (19). The scalar S' in (19) is a function 


of x alone. 
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The v-curvature tensor of any four-dimensional Finsler space may be written as [13] 
LD? Shigk = Ogu) {Rng Keni + hirKny} (20) 


where K;; is a (0, 2) type symmetric Finsler tensor field which is such that Ki;y? = 0. A Finsler 
space F'”(n > 4) is called $4-like Finsler space [13] if its v-curvature tensor is of the form (20). 
From (17), (19), (20) and (5) we have the following theorems. 


Theorem 3.1 The m*”-root Randers changed S3-like or S4-like Finsler space is S4-like Finsler 


Space. 


Theorem 3.2 If v-curvature tensor of m‘’-root Randers changed Finsler space F” vanishes, 
then the Finsler space with m*"-root metric F” is S4-like Finsler space. 


If v-curvature tensor of Finsler space with m*’-root metric F” vanishes then equation (17) 
reduces to 
Snigk = hijmne + harmiz — hikmMng — hngmir (21) 
By virtue of (21) and (11) and the Ricci tensor Six = 9" Syijx is of the form 
1 


Six = Gump + (m = 1)(n = 3)mix}, 


where m = m,ja"), which in view of (18) may be written as 





Sik + Hihin + HoCigrd” = H3ymimn, (22) 
where 

m (n — 3)(? — @) 
H. = ———_ Sa =P 
Gale Bone. ° 

Hels 

2Xm—1)L 

(n—3) 
oe 


From (22), we have the following 


Theorem 3.3 If v-curvature tensor of m‘"-root Randers changed Finsler space F" vanishes 
then there exist scalar H, and Hy in Finsler space with m*”-root metric F"(n > 4) such that 
matrix || Six + Hihin + HoCixpb” || is of rank two. 


§4. The (v)hv-Torsion Tensor and hv-Curvature Tensor of F” 


Now we concerned with (v)hv-torsion tensor P;;;, and hv-curvature tensor Pp;;,. With respect 
to the Cartan connection CT, Ly; = 0, li; = 0, hijjx = 0 hold good [13). 
Taking h-covariant derivative of equation (9) and using (4) and J; = a; = 0 we have 
- bith 


Cijhin = TCighin + F—Cijk + 


(higbejn + hyxdijn + heidsin) 


OL (23) 
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From equation (6) and using relation h,,\, = 0 We have 


2LC ijn 
Qijjn =0, and aizgin = aM (24) 


m—1(m—2) 


The (v)hv-torsion tensor P;;, and the hv-curvature tensor Ppi;~ of the Cartan connection CT 


are written in the form, respectively 
Pijk = Cigrjo, (25) 
Phigk = Cigrih — Chgr|i + Pier Ch, — Phar Chi 


where the subscript ’0’ means the contraction for the supporting element y’. Therefore the 
(v)hv-torsion tensor Pi5x and the hv-curvature tensor Prijk of the Cartan connection CT for 
the Finsler space with m‘"-root Randers metric by using (10), (23), (24) and (25) we have 


5 (m — 1)(m — 2) bilo (hijbejo + hjadijo + hridjio) 
Pijk = 97 aiaR|0 an Tz Ciik + ees aes (26) 
and 
Prigk = (m — 1)(m — 2)(2L)* Oxy (Gajajn + PitrChn) (27) 


Definition 4.1([13]) A Finsler space is called a Berwald space (resp. Landsberg space) if 
Cijrin = 9 (resp. Pijz = 0) holds good. 


Consequently, from (24) and (26) we have 


Theorem 4.1 A Finsler space with the m*”-root Randers changed metric is a Berwald space 
(resp. Landsberg space), if and only if ayjrjn = 0 (resp. Ajjnj0 = 0 and by, is covariently 


constant. 


Proposition 4.1 The v(hv)-torsion tensor and hv-curvature tensor Pnigk of m'"-root Randers 
changed Finsler space F” with respect to Cartan’s connection CT is of the form (26) and (27). 


§5. T-Tensor of F” 


Now, the T-tensor is given by [11,13] 

Thigk = LCnijle t+ liCnjr + ly Crin + laCniy + lnCijr 
The above equation for m‘”-root Randers changed Finsler space F'” is given as 

Thijgk = LCnij|k + iCnjz + Gy Cnin + leCnij + nCijx (28) 
The v-derivative of h;; and L is given by [13] 


1 
hig | => ~ 7 (hanks + hykli), and Li = l; (29) 
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Now using (29), the v-derivative of Ci;, is given as 


Lbp — Bl = 1 
en Py + Lr Cin |n - Top (hindj7mx + hjnlime (30) 


+hjnlymi Tr Aenljymi + hinlkm; + hrenlim,; + hijlame 


LCigklh = 





i 
+hyjglnmi + heilam;) + 5 (his mela + hypmiln + hesm|n) 


Using (4), (9) and (30), the T-tensor for m*”-root Randers changed Finsler space F” is given 
by 


as T 
Thijk = T(Thijk + Bhajk) + ap (hainmali + hirmal; + hijmnle (31) 
1 
+hgimyln) + ap (hnamndi + hjrmnds + hen; by + higmnd; 
+hinmnd; + henmid; + higmnbr + hinmybr + hjnmibs + higmrbn 
i 
+h bn + hjpmibn) + 5 (higmaln + hypmiln + hesmy|n) 


(Lon = Bln) 
T i C; k 





where Bnijk = BChij e+ biCnjk + b;Chik + beChnij + brCijr- Thus, we know 


Proposition 5.1 The T-tensor Thijk for m*"-root Randers changed Finsler space F” is given 
by (31). 
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Abstract: The object of this paper is to introduce a quarter-symmetric metric connec- 
tion in a pseudosymmetric Lorentzian a-Sasakian manifold and to study of some proper- 
ties of it. Also we shall discuss some properties of the Weyl-pseudosymmetric Lorentzian 
a—Sasakian manifold and Ricci-pseudosymmetric Lorentzian a—Sasakian manifold with re- 
spet to quarter-symmetric metric connection. We have given an example of pseudosymmetric 


Lorentzian a-Sasakian manifold with respect to quarter-symmetric metric connection. 


Key Words: Lorentzian a-Sasakian manifold, quarter-symmetric metric connec- 
tion, pseudosymmetric Lorentzian a—Sasakian manifolds, Ricci-pseudosymmetric, Weyl- 
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§1. Introduction 


The theory of pseudosymmetric manifold has been developed by many authors by two ways. 
One is the Chaki sense [8], [3] and another is Deszcz sense [2], [9], [11]. In this paper we 
shall study some properties of pseudosymmetric and Ricci-symmetric Lorentzian a— Sasakian 
manifolds with respect to quarter-symmetric metric connection in Deszcz sense. The notion 
of pseudo-symmetry is a natural generalization of semi-symmetry, along the line of spaces of 
constant sectional curvature and locally symmetric space. 

A Riemannian manifold (M,g) of dimension n is said to be pseudosymmetric if the Rie- 


mannian curvature tensor R satisfies the conditions ((1]): 


for all vector fields X, Y,U, V,W on M , whereLp € C*(M), R(X, Y)Z = Vix,yjZ—-[V x, Vy]Z 
and X AY is an endomorphism defined by 


(X AY)Z = g(Y, Z)X — g(X, Z)Y (2) 


1Received October 6, 2012. Accepted March 6, 2013. 
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2. (R(X, Y).R)(U,V,W) = R(X,Y)(R(U,V)W) — R(R(X,Y)U, V)W 

—RU, R(X,Y)V)W — RU, V)(R(X,Y)W) (3) 
3. ((XAY).R\(U,V,W) = (XAY)(R(U,V)W) — R(X AY)U,V)W 

—RU(X AY)V)W — RU,V)\(X AY)W). (4) 


M is said to be pseudosymmetric of constant type if DL is constant. A Riemannian manifold 
(M, g) is called semi-symmetric if R.R = 0, where R.R is the derivative of R by R. 


Remark 1.1 We know, the (0, + 2) tensor fields R.T and Q(g,T) are defined by 


(RDS Ap eV) = CACY) PF) age Xe) 
=. SPURAGY Kiyo hg) eT Aagt RG VY XB) 
Og; T)( Migr pA AY oS AY) DOG, ee 5 2G) 


= T((XAY)X,-++,Xn) +--+ T(X, +++ (KX AY) Xx), 


where T is a (0,k) tensor field ([4],[5]). 


Let S and r denote the Ricci tensor and the scalar curvature tensor of M respectively. The 
operator Q and the (0,2)—tensor S$? are defined by 


S(X,Y) = 9(QX,Y) (5) 


and 
S?(X,Y) = S(QX,Y) (6) 


The Weyl conformal curvature operator C' is defined by 
1 
O(X,Y) = R(X,Y) - —5IXAQY +QX AY -—_XAyY]. (7) 
n— n— 


If C = 0, n > 3 then M is called conformally flat. If the tensor R.C and Q(g,C) are linearly 
dependent then M is called Weyl-pseudosymmetric. This is equivalent to 


R.CUU,V,W; X,Y) = Lol[((X AY).C)(U,V)W], (8) 
holds on the set Uc = {x € M:C #0 at x}, where Lc is defined on Uc. If R.C = 0, then M 


is called Weyl-semi-symmetric. If VC = 0, then M is called conformally symmetric ([6],[10]). 


§2. Preliminaries 


A n-dimensional differentiable manifold M is said to be a Lorentzian a—Sasakian manifold if 
it admits a (1,1)—tensor field ¢, a contravariant vector field €, a covariant vector field 7 and 
Lorentzian metric g which satisfy the following conditions, 


o =14+7@®E, (9) 
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IX, OY) = G(X, Y) + (X)n(¥), (11) 
WX, £) = (X) (12) 

and 
(Vx@)(Y) = atg(X, YE + n(V)X} (13) 


for VX,Y € x(M) and for smooth functions a on M, V denotes covariant differentiation 
operator with respect to Lorentzian metric g ([6], [7]). 
For a Lorentzian a—Sasakian manifold,it can be shown that ([6],[7]): 


Vxé = adX, (14) 
(Vxn)Y ag(ox,Y). (15) 


I 


Further on a Lorentzian a—Sasakian manifold, the following relations hold ((6]) 


n(R(X,Y)Z) = a? [g(¥, Z)n(X) — 9(X, Z)n(V)], (16) 
R(EX)Y = a°[g(¥, Z)E—n(¥)X], (17) 
R(X,Y)E = a*[n(¥)X — n(X)Y], (18) 

S(E,X) = S(X,£) =(n— 1)a7n(X), (19) 
S(E,€) = —(n—-1)o’, (20) 
QE = (n—1)o°%E. (21) 


The above relations will be used in following sections. 


§3. Quarter-Symmetric Metric Connection on Lorentzian a—Sasakian Manifold 


Let M bea Lorentzian a-Sasakian manifold with Levi-Civita connection V and X,Y, Z € x(M). 
We define a linear connection D on M by 


Dx¥ =Vx¥ +n(¥)(X) (22) 


where 7 is 1—form and ¢ is a tensor field of type (1,1). D is said to be quarter-symmetric 


connection if 7’, the torsion tensor with respect to the connection D, satisfies 


T(X,Y) = (VY )bX — 7(X) oY. (23) 
D is said to be metric connection if 
(Dxg)(Y, Z) = 0. (24) 


A linear connection D is said to be quarter-symmetric metric connection if it satisfies (22), (23) 
and (24). 
Now we shall show the existence of the quarter-symmetric metric connection D on a 


Lorentzian a—Sasakian manifold M. 
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Theorem 3.1 Let X,Y,Z be any vectors fields on a Lorentzian a—Sasakian manifold M and 


let a connection D is given by 
2g(DxY,Z) = Xg(Y,Z)+YG(Z,X) -— Zg(X,Y) + o([X, Y], Z) 
—9(¥, Z], X) + 914, X],Y) + 9(n(V)oX — n(X)6Y, Z) 
+9(MX)oZ — 0(Z)bX,Y) + g(n(¥ )bZ — 0(Z)@Y, X). (25) 
Then D is a quarter-symmetric metric connection on M. 


Proof It can be verified that D: (X,Y) — DxY satisfies the following equations: 





Dx(¥+Z) = DxY+DxZ, (26) 
DxsyZ = DxZ+DyZ, (27) 
DyxY = fDxyY, (28) 
Dx(fY) = f(DxY)+(Xf)Y (29) 


for all X,Y, Z € x(M) and for all f, differentiable function on M. 
From (26), (27), (28) and (29), we can conclude that D is a linear connection on M. From 
(25) we have, 


g(DxY, Z) — g(Dy X, Z) = g([X, Y], Z) + nV )a(bX, Z) — 1(X)g9 (GY, Z) 


or, 
DxY — DyX — [X,Y] = (Y)bX — (X) oY 


or, 


T(X,Y) =n(¥ )oX — (X)o¥ (30) 
Again from (25) we get, 


2g(DxY, Z) + 2g(DxZ,Y) =2Xg(Y,Z), or, (Dxg)(Y,Z) =0. 











This shows that D is a quarter-symmetric metric connection on M. 





§4. Curvature Tensor and Ricci Tensor with Respect to Quarter-Symmetric 


Metric Connection D in a Lorentzian a—Sasakian Manifold 


Let R(X,Y)Z and R(X,Y)Z be the curvature tensors with respect to the quarter-symmetric 
metric connection D and with respect to the Riemannian connection V respectively on a 
Lorentzian a— Sasakian manifold M. A relation between the curvature tensors R(X,Y)Z and 
R(X, Y)Z on M is given by 


R(X,Y)Z = R(X,Y)Z+alg(X, Z)6Y 
—9 (PY, 2)6X] + an(Z)[n(¥)X — n(X)Y]. (31) 


Also from (31), we obtain 


S(X,Y) = S(X,Y) + alg(X,¥) + nn(X)n(¥)), (32) 
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where S$ and S are the Ricci tensors of the connections D and V respectively. 
Again 
S7(X,Y) = S?(X,Y)—-a(n—2)8(X,Y) —a?(n—-1)9(X,Y) 
ta?n(n — 1)(a — 1)n(X)n(Y). (33) 





Contracting (32), we get 
Pat: (34) 


where 7 and r are the scalar curvature with respect to the connection D and V respectively. 
Let C be the conformal curvature tensors on Lorentzian a— Sasakian manifolds with 
respect to the connections D. Then 


1 


n—- 





C(X,Y)Z = R(X,Y)Z- 5 [S(Y, Z)X — g(X, Z)QY + 9({¥, Z)QX 


(n — 1)(n — 2) 


where Q is Ricci operator with the connection D on M and 


—S(X, Z)Y] + [o(Y, Z)X — g(X, Z)Y], (35) 
S(X,Y) = 9(QX,Y), (36) 


SCY) =S(OX, VY). (37) 


Now we shall prove the following theorem. 


Theorem 4.1 Let M be a Lorentzian a—Sasakian manifold with respect to the quarter- 
symmetric metric connection D, then the following relations hold: 


REE X)¥Y = a®[g(X,Y)E- (VY) X] + an(V)[X + (XQ, (38) 
M(R(X,Y)Z) = a? [g(¥, Z)n(X) — 9({X, Z)n(¥)], (39) 
R(X,Y)E = (a? —a)[n(¥)X — (X)Y], (40) 
S(X,£) = S(E,X) =(n—1)(e? — a)n(X), (41) 
S*(X,€) = S*(€,X) =a*(n—1)*(a— 1)’n(X), (42) 
568 = —(n—-1)(e’—a), (43) 

QX = QX-a(n-1)X, (44) 

QE = (n-1)(a? —a)é. (45) 


Proof Since M is a Lorentzian a—Sasakian manifold with respect to the quarter-symmetric 
metric connection D, then replacing X = € in (31) and using (10) and (17) we get (38). Using 
(10) and (16), from (31) we get (39). To prove (40), we put Z = € in (31) and then we use 
(18). Replacing Y = € in (32) and using (19) we get (41). Putting Y = € in (33) and using (6) 
and (19) we get (42). Again putting X = Y = € in (32) and using (20) we get (43). Using (36) 
and (41) we get (44). Then putting X = € in (44) we get (45). 
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§5. Lorentzian a—Sasakian Manifold with Respect to the Quarter-Symmetric 
Metric Connection D Satisfying the Condition C.S = 0. 


In this section we shall find out the characterization of Lorentzian a—Sasakian manifold with 
respect to the quarter-symmetric metric connection D satisfying the condition C.S = 0. We 
define C.S = 0 on M by 

where X,Y,Z,W € x(M). 


Theorem 5.1 Let M be an n—dimensional Lorentzian a— Sasakian manifold with respect to 


the quarter-symmetric metric connection D. If C.S = 0, then 








i a2 2 ‘a 
+——S a(n — 1)(a—n +1) - Flg(X,Y) 
—a(n — 1)(a* — a)n(X)n(¥). (47) 


Proof Let us consider M be an n-dimensional Lorentzian a—Sasakian manifold with 
respect the quarter-symmetric metric connection D satisfying the condition C.$ = 0. Then 
from (46), we get 


S(C(X,Y)Z,W) + S(Z,C(X, Y)W) =0, (48) 


where X,Y,Z,W € y(M). Now putting X = € in (48), we get 
5(C(E, XY, Z) + BY, OE, X)Z) =0. (49) 


Using (35), (37), (38) and (41), we have 




















S(CEX)Y.Z) = (n= 1)(0? ~a)[a? —- SI) 4, inayat x.) 
2, (nla? — a) : : 
+la -—-avct+ =e = GoDMsa ss Z) 
+a(a? = a)(n = 1)n(X)n(V)n(Z) 
—n—1)(6 ~ a)n(Z)5(X,Y) ~ 8X, Z)n(V) (50) 
and 
SW.C(,X)Z) = (n= 1a? - ofa? - PE) into x, 2) 
2, (n= Ila? — a) ; E 
+[a—-—a*+ ae “Gotan oe 
+a(a? = a)(n = 1)n(X)n(V)n(Z) 
~—1 [n= 1)(0? — a) nl ¥)5(X, Z) — $°(X, Y)n(Z)] (51) 


n—2 


52 C.Patra and A.Bhattacharyya 


Using (50) and (51) in (49), we get 














2 2_ (n—1)(e? —a) . 
(n— 1)(a* —a)| = in Iy(n — ay lo Vn) 
+9(X, Z)n(Y)] + 2a(a? — a)(n — 1)n(X)n(Y )n(Z) 
2, (n—1)(* -a) r i é 
SO. gas mopman is, Z) +(Z)S(Y, X)| 
ppm (es 1)(a? — a){n(Z)S(X,Y) + n(¥)S(X, Z)} 
—{5°(X, Z)n(Y) + $?(X, Y)n(Z)}] = 0. (52) 


Replacing Z = € in (52) and using (41) and (42), we get 


SX oY) = le) a5) 





n—2 ) 
——~ [a(n — 1)(a— n+ 1) ~ FIg(X,Y) 


—a(n — 1)(a* — a)n(X)n(Y). 


+ 




















An n—dimensional Lorentzian a—Sasakian manifold M with the quarter-symmetric metric 
connection D is said to be n—Einstein if its Ricci tensor $ is of the form 


S(X,Y) = Ag(X,Y) + Bn(X)n(Y), (53) 


where A, B are smooth functions of M. Now putting X = Y = e;, i = 1,2,--- ,m in (53) and 


taking summation for 1 <i <n we get 
An— B=#. (54) 
Again replacing X = Y = € in (53) we have 
A-—B=(n-1)(a? —- a). (55) 
Solving (54) and (55) we obtain 


_ ey acen op pte 
ear (a“ —a) and B aa n(a“ — a). 





Thus the Ricci tensor of an 7—Einstein manifold with the quarter-symmetric metric connection 
D is given by 


5(X,¥) =[— - (a? — a)lg(X,¥) + [5 — na? — a)n(X)n(¥). (56) 


n—-1 n—- 








§6. »1—Einstein Lorentzian a—Sasakian Manifold with Respect to the 
Quarter-Symmetric Metric Connection D Satisfying the Condition C.S = 0. 


Theorem 6.1 Let M be an n—Einstein Lorentzian a—Sasakian manifold of dimension. Then 
C.S =0 iff 
na — 2a 


In R(X, Y)Z)n(W) + (R(X, Y)W)n(Z)] = 0. 
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where X,Y,Z,W € x(M). 
Proof Let M be an 7—Einstein Lorentzian a—Sasakian manifold with respect to the 
quarter-symmetric metric connection D satisfying C.S = 0. Using (56) in (48), we get 


mC(X,Y)Z)n(W) + n(C(X, Y)W)n(Z) = 0, 


or 


[n(R(X,Y)Z)n(W) + n( R(X, Y)W)n(Z)] = 0. 


Conversely, using (56) we have 


r 2 




















(C(X,Y).8)\(Z,W) = -[— - nla? — a)]ln(C(X,Y)Z)n(W) + n(E(X,Y)W)n(Z) 
= MOS nl RUX,Y)Z)n(W) + (R(X, Y)W)n(Z)] = 0. 


§7. Ricci Pseudosymmetric Lorentzian a—Sasakian Manifolds with 


Quarter-Symmetric Metric Connection D 


Theorem 7.1 A Ricci pseudosymmetric Lorentzian a- Sasakian manifolds M with quarter- 
symmetric metric connection D with restriction Y = W = € and Lg = 1 ts an n—Einstein 


manifold. 


Proof Lorentzian a—Sasakian manifolds M with quarter-symmetric metric connection D 


is called a Ricci pseudosymmetric Lorentzian a—Sasakian manifolds if 
(R(X, Y).S)(Z, W) = Lg[((X AY).8)(Z,W)], (57) 
or 
S(R(X,Y)Z,W) + 8(Z, R(X, Y)W) = La[S((X AY)Z,W) + S(Z,(X AY)W)). (58) 
Putting Y = W = € in (58) and using (2), (38) and (41), we have 


Lg[$(X, Z) — (n— 1)(a? — a)g(X, Z)| 


= (a* — a)S(X, Z) — a*(a* — a)(n — 1)9(X, Z) — a(a? — a)(n — 1)n(X)n(Z). (59) 








Then for Lz = 1, 





(a? — a —1)8(X, Z) = (a? — a)(n — 1)[(a? — 1)9(X, Z) + an(X)n(Z)]. 











Thus M is an 7—Einstein manifold. 





Corollary 7.1 A Ricci semisymmetric Lorentzian a-Sasakian manifold M with quarter-symmetric 


metric connection D with restriction Y =W = € is an n—Einstein manifold. 
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Proof Sine M is Ricci semisymmetric Lorentzian a-Sasakian manifolds with quarter- 


symmetric metric connection D, then Le = 0. Putting Le = 0 in (59) we get 











S(X, Z) = a?(n— 1)g(X, Z) + a(n — 1)n(X)n(Z). 





§8. Pseudosymmetric Lorentzian a—Sasakian Manifold and Weyl-pseudosymmetric 


Lorentzian a—Sasakian Manifold with Quarter-Symmetric Metric Connection 


In the present section we shall give the definition of pseudosymmetric Lorentzian a—Sasakian 
manifold and Weyl-pseudosymmetric Lorentzian a—Sasakian manifold with quarter-symmetric 


metric connection and discuss some properties on it. 


Definition 8.1 A Lorentzian a—Sasakian manifold M with quarter-symmetric metric connec- 
tion D is said to be pseudosymmetric Lorentzian a—Sasakian manifold with quarter-symmetric 


metric connection if the curvature tensorR of M with respect to D satisfies the conditions 


(R(X, Y).R)(U,V,W) = Lal((X AY).R)(U,V,W)}, (60) 
where 
(R(X,Y).R)(U,V,W) = R(X,Y)(R(U,V)W) — R(R(X,Y)U,V)W 
—R(U, R(X, Y)V)W — R(U,V)(R(X,Y)W), (61) 
and 
((X AY).R)(U,V,W) = (XAY)(R(U,V)W) — R((X AY)U,V)W 
—R(U,(X AY)V)W — R(U,V)((X AY)W). (62) 


Definition 8.2 A Lorentzian a—Sasakian manifold M with quarter-symmetric metric con- 
nection D is said to be Weyl- pseudosymmetric Lorentzian a—Sasakian manifold with quarter- 


symmetric metric connection if the curvature tensorR of M with respect to D satisfies the 


conditions 
(R(X, Y).C)(U,V,W) = Le[((X AY).C)U,V,W)], (63) 
where 
(R(X, Y).C)\(U,V,W) = R(X,Y)(C(U,V)W) — C(R(X,Y)U, V)W 
—C(U, R(X, Y)V)W — C(U,V)(R(X,Y)W) (64) 
and 
((X AY).C)\U,V,W) = (XAY)(C(U,V)W)—C((X AY)U,V)W 
—C(U,(X AY)V)W — C(U,V)((X AY)W). (65) 


Theorem 8.1 LetM be ann dimensional Lorentzian a—Sasakian manifold. If M is Weyl- 


pseudosymmetric then M is either conformally flat and M is n—Einstein manifold or Le = a?. 
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Proof Let M be an Weyl-pseudosymmetric Lorentzian a—Sasakian manifold and X,Y, 
U,V,W € x(M). Then using (64) and (65) in (63), we have 


R(X, Y)(G(U, V)W) — G(R(X, Y)U, VW 


—C(U, R(X, Y)V)W — C(U,V)(R(X,Y)W) 
= Lel(X AY)(G(U, V)W) — G(X AY)U,V)W (66) 
_O(U,(X AY)V)W — G(U,V)((X AY)W)]. (67) 


Replacing X with € in (66) we obtain 


R(E,Y)(CU,V)W) — C(R(E,Y)U,V)W 

—C(U, RE, Y)V)W - CU,V)(RE,Y)W) 

= Lel(ENY)(C(UU,V)W) — C(EAY)U, VW 
—CUU,(EAY)V)W —CU,V)(EAY)W)).- (68) 


Using (2), (38) in (67) and taking inner product of (67) with €, we get 


a*|—C(U, V,W,Y) — n(CU,V)W)n(¥) — 9(¥, U)n(C(E, VW) 





+n(U)n(C(¥, VW) — g(¥, V)n(CU, EW) + nV )n(CU,Y)W) 
+(W)n(CU, V)Y)] — alnU)n(C*Y, V)W) 

+(V)n(CU, YW) + o(W)n(CU,V)¢*Y)] 

= Le[-C¥, U,V, W) — n(¥)n(CU,V)W) — g{¥, U)n(C(E, VW) 
+n(U)n(C(¥, VW) — g(¥, V)n(C(U, EW) + nV (CU, Y)W) 
+(W)n(C(U, V)Y)]. (69) 


Putting Y = U, we get 
[Le — a” |[9(U, U)n(C(E, VW) + 9U,V)n(CU,€)W)] + an(V)n(C(#U,V)W) =0. (70) 
Replacing U = € in (68), we obtain 
[Le — a*|n(C(E,V)W) =0. (71) 


The formula (69) gives either n(C(€,V)W) = 0 or Le — a? = 0. 
Now Le — a? £ 0, then n(C(é,V)W) = 0, then we have M is conformally flat and which 
gives 
S(V,W) = Ag(V,W) + Bn(V)n(W), 


where (0? 
and ‘ i 
ee ee Ce 


which shows that M is an n—Einstein manifold. Now if n(C(é,V)W) 40, then Le = a?. 
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Theorem 8.2 LetM be an n dimensional Lorentzian a—Sasakian manifold. If M is pseu- 
dosymmetric then eitherM is a space of constant curvature and ag(X,Y) = n(X)n(Y), for 
a#0 or Lp =a’, for X,Y € y(M). 


Proof Let M be a pseudosymmetric Lorentzian a—Sasakian manifold and X,Y,U,V,W © 
x(/). Then using (61) and (62) in (60), we have 


R(X, Y)(R(U, V)W) — R(R(X,Y)U, V)W 


—R(U, R(X, Y)V)W — R(U,V)(R(X, Y)W) 
= Lal(X AY)(RU,V)W) — R(X AY)U,V)W 
—R(U,(X AY)V)W — RU,V)(X AY)W)). (72) 


Replacing X with € in (70) we obtain 
R(E,Y)(R(U,V)W) — R(R(E,Y)U,V)W 

—R(U, R(E,Y)V)W — RU,V)(R(E,Y)W) 

= Lal(EAY)(RU,V)W) — R(EAY)U,V)W 

—RUU,(EAY)V)W — RU,V)(EAY)W)]. (73) 
Using (2), (38) in (71) and taking inner product of (71) with €, we get 

a*|—R(U,V,W,Y) — (RU, V)W)n(¥) — g(¥, U)n( RE, V)W) 

+n(U)n(R(Y,V)W) — oY, V)n( RU, 


g)W) + (V)n(RU,Y)W) 
+n(W)n(RU, V)¥)] — a[n(U)n( RY, VW) 
+n(V)n(RU, oY )W) + n(W)n(RU, V)¢°Y)| 
= Lal-R(Y,U,V,W) — n(¥ )n(R(U, VW) — oY, U)n(R(E,V)W) 
+n(U)n(R(Y, VW) — oY, V)n( RU, EW) + n(V)n(RU, Y)W) 
+n(W)n(RU,V)Y))- 


Putting Y = U, we get 
[La — o°][9(U, U)n(R(E,V)W) + g(U, V)n( RU, E)W)] + an(V)n( RGU, V)W) =0. (74) 
Replacing U = € in (72), we obtain 
[La — a7 |n(R(E,V)W) = 0. (75) 


The formula (73) gives either n(R(€,V)W) = 0 or Lp — a? = 0. Now Lp — a? F 0, then 
n(R(E,V)W) = 0. We have M is a space of constant curvature and n(R(€,V)W) = 0 gives 
ag(V,W) = 7(X)n(Y) for a 4 0. If n(R(E,V)W) 4 0, then we have Lp = a?. 














§9. Examples 


Let us consider the three dimensional manifold M = {(x1, 72,73) € R°® : x1, x2, 73 € R}, where 
(a1, 22,23) are the standard coordinates of R?. We consider the vector fields 


3) a) 7) 
ey = e§ —, eg =e*3(— + —) and e3 =a—, 


Ox Ox 0x3 
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where a is a constant. 
Clearly, {e1, e2,e3} is a set of linearly independent vectors for each point of M and hence 
a basis of y(M). The Lorentzian metric g is defined by 


(e1, €3) = 0, 


(es, €3) =-l. 


g(e1, €2) (€2, €3) 


= G9 = 9 
g(e1, €1) = g(e2,e2) =9 


Then the form of metric becomes 


eS 


aes ae 
g= (eray2 (2) 2 (dx3)", 


which is a Lorentzian metric. 
Let 7 be the 1—form defined by n(Z) = g(Z,e3) for any Z € x(M) and the (1, 1)—tensor 
field @ is defined by 


ge, =—€1, beg =—e2, ez =0. 
From the linearity of ¢ and g, we have 
n(es) = —1, 
o?(X)=X+n(X)ez and 
G(X, OY) = g( X,Y) + 0(X)n(V) 
for any X € y(M). Then for e3 = €, the structure (¢,,7,g) defines a Lorentzian paracontact 


structure on M. 


Let V be the Levi-Civita connection with respect to the Lorentzian metric g. Then we have 
[e1,e2] =0, [e1,e3] =—@e1, [e2, e3] = —aer. 
Koszul’s formula is defined by 


2g(VxY,Z) = XgGg(Y,Z)+YG(Z, xX) — Zg(X,Y) 
—9(X, ba Z)) _ GY, [X, Z)) + g(Z, [X, |): 


Then from above formula we can calculate the followings, 


Vex €1 = —a€3, Vex €2 = 0, V e163 = —aeéq, 
Ve.é1 =0, Ve,e2 = —ae3, Ve,e3 = —ae2, 


Vex€1 = 0, V eg €2 => 0, Ve3€3 = 0. 


Hence the structure (¢,£,7,g) is a Lorentzian a—Sasakian manifold [7]. 
Using (22), we find D, the quarter-symmetric metric connection on M following: 


De,e1 = —ae3, De,eg =0, De,e3 = e1(1— a), 
D.,€1 =0, De,e2 = —ae3, De,e3 = €2(1— a), 


Deze = 0, Dez €2 = 0, Dez €3 = 0. 
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Using (23), the torson tensor J’, with respect to quarter-symmetric metric connection D as 


follows: 
T (ei, €:) = 0, Vi= 1,2,3 
T(e1, €2) = 0, T(e1, €3) =€1, T (ea, €3) = €2. 
Also (De,g)(€2,e€3) = (Deog)(e3, €1) = (Desg)(e1, €2) = 0. Thus M is Lorentzian a—Sasakian 


manifold with quarter-symmetric metric connection D. 
Now we calculate curvature tensor R and Ricci tensors S as follows: 


R(e1,e2)e3 = 0, R(e1,e3)e3 = —(a” — a)ei, 
R(e3, €2)e2 = a%e3, R(e3,€1)e1 = aes, 
R(e2,€1)e1 = (a? — a)ez, Rez, e3)e3 = —a7e2, 
R(e1, €2)e2 = (a? — ajey 

(e1, €1) 


= —a and S(e3,e3) = —2a7 + (n—1)a. 


wi 
o 
os 
® 
AS: 
ll 
wi 
ae 
@ 
dN 
® 
NO 
Sy 


Again using (2), we get 


(e1, e€2)e3 — 0, (e A €4)e; = 0, Vij = 1, 2, 3, 
(e1 A e2)e2 = (€1 A e3)e3 = —€1, (€2 A €1)€1 = (€2 A €3)€3 = —e2, 


(e3 A €2)e€2 = (e3 A €1)e1 = —€3. 


R(ex e2)(R R(e3, €1)€2) = = 0, R(R(ex, €2)e3, €1)€2 =0, 
R(es, R(e1, e2)e1)e2 = —a*(a? — a)es, 
(Rees, e1)(R(e1, e2)e2) = a?(a? — a)e3. 
Therefore, (R(e1, e2).R)(e3, €1, €2) = 0. 
Again, 
(e1 A e2)(R(e3, €1)€2) = 0, R((ex \ €2)e3, €1)€2 — 0, 
Res, (e1 A €2)e€1)e2 = aes, R(es, e1)((e1 \ €2)e2) — —a’es. 
Then ((e1, e2).R)(e3, €1,e2) = 0. Thus (R(e1, e2).R)(e3, e1, €2) = Lpl((e1, e2).R)(e3, 1, €2)] for 
any function = Lz € C~°(M). 


Similarly, any combination of e1,e2 and e3 we can show (60). Hence M is a pseudosym- 


metric Lorentzian a—Sasakian manifold with quarter-symmetric metric connection. 
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Abstract: This paper is motivated by the skew energy of a digraph as vitiated by C.Adiga, 
R.Balakrishnan and Wasin So [1]. We introduce and investigate the skew energy of a Cayley 
digraphs of cyclic groups and dihedral groups and establish sharp upper bound for the same. 
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§1. Introduction 


Let G be anon trivial finite group and S be an non-empty subset of G such that for x € S,x~! ¢ 
S and Ig ¢ S, then the Cayley digraph T = Cay(G, S) of G with respect to S is defined as a 
simple directed graph with vertex set G and arc set E([) = {(g,h)|hg~! € S}. If S is inverse 
closed and doesn’t contain identity then Cay(G, S) is viewed as undirected graph and is simply 
the Cayley graph of G with respect to S. It easily follows that valency of Cay(G, S) is |S| and 
Cay(G,S) is connected if and only if (S) = G. For an elaborate literature on Cayley graphs 
one may refer [5]. A dihedral group D2, is a group with 2n elements such that it contains an 
element ‘a’ of order 2 and an element 'b’ of order n with a~tba = b~!. Thus Don, = (a, bla? = 
b” = 1,a~*ba = db!) = (a, bla? = b” = 1,a7'ba = b°,a# 1(mod_ n),a2 =1(mod_ n)). 

If n = 2, then Dy, is Abelian; for n > 3, Don is not abelian. The elements of diheral group 
can be explicitly listed as 


Dan = {1, a, ab, ab, pare Jab}, b, b?, Pos be}, 


In short, its elements can be listed as a’b* where i = 0,1 and k = 0,1,--- ,(n —1). It is easy to 
explicitly describe the product of any two elements a’b*aJb! = a”b’ as follows: 


1. If 7 =0 then r =7 and s equals the remainder of & + 1 modulo n. 


2. If j =1, then r is the remainder of i+ 7 modulo 2 and s is the remainder of ka + | 
modulo n. 


1Received October 16, 2012. Accepted March 8, 2013. 
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The orders of the elements in the Dihedral group Dz, are: o(1) = 1, o(ab') = 2, where; 
0<i<n-1, o(b') =n, where; 0 <i<n-—1 and if n is even than o(b?) = 2. 

Let I’ be a digraph of order n with vertex set V([) = {vi,--- ,Un},and are set A(T) C 
V(T) x V(T). We assume that [’ does not have loops and multiple arcs, i.e.,(v;,u;) ¢ A(T) for 
alli , and (v,,v,;) € A(L) implies that (v;,v;) ¢ A([). Hence the underlying undirected graph 
Gr of I is a simple graph. The skew-adjacency matrix of I is the n x n matrix S(T) = [s,,], 
where s;; = 1 whenever (v;,v;) € A(I),si; = —1 whenever (v;,u;) € A(T), and si; = 0 
otherwise. Because of the assumptions on I’, S(T) is indeed a skew-symmetric matrix. Hence 
the eigenvalues {A1,--+ , An} of S(T) are all purely imaginary numbers, and the singular values 
of S([) coincide with the absolute values {|\i|,--+ ,|An|} , of its eigenvalues. Consequently, the 
energy of S(T) , which is defined as the sum of its singular values [6], is also the sum of the 
absolute values of its eigenvalues. For the sake of convenience, we simply refer the energy of 
S(T) as the skew energy of the digraph I. If we denote the skew energy of I by ¢,(I) then, 


es(T) = > [Ai]. 

The ie of a vertex in a digraph T is the degree of the corresponding vertex of the 
underlying graph of [. Let D(T) = diag(d1, d2,--- , dn), the diagonal matrix with vertex degrees 
dy, dz,+-+ ,dy of v1,V2,-++,Un and S(T) be the skew adjacency matrix of a simple digraph I, 
possessing n vertices and m edges. Then L(T) = D(T) — S(L) is called the Laplacian matrix of 
the digraph T. If \;,7 = 1,2,--- ,n are the eigenvalues of the Laplacian matrix L(T) then the 


. 2 
skew Laplacian energy of the digraph T is defined as SLE(T) = y [Ay — cel 
n 
i=1 


An nxn matrix S is said to be a circulant matrix if its entries satisfy s;; = s1,;-1+1 , where 
the subscripts are reduced modulo n and lie in the set {1,2,...,n}. In other words, ith row of S$ 
is obtained from the first row of S by a cyclic shift of 1 — 1 steps, and so any circulant matrix is 

n 


determined by its first row. It is easy to see that the eigenvalues of S are A, = S- 81jw9-DR 


j=l 
k =0,1,---,n—1. For any positive integer n, let t, = {w* :0 < k < n} be the set of all nth 
roots of unity, where w =e = cos(#=) + isin(=*) that i? =—1. t is an abelian group with 


respect to multiplication. A circulant graph is a graph I whose adjacency matrix A(T) is a 
circulant matrix. More details about circulant graphs can be found in [3). 

Ever since the concept of the energy of simple undirected graphs was introduced by Gutman 
in [7], there has been a constant stream of papers devoted to this topic. In [1], Adiga, et al. have 
studied the skew energy of digraphs. In [4], Gui-Xian Tian, gave the skew energy of orientations 
of hypercubes. In this paper we introduce and investigate the skew energy of a Cayley digraphs 
of cyclic groups and dihedral groups and establish sharp upper bound for the same. 


§2. Main Results 


First we present some facts that are needed to prove our main results. 


Lemma 2.1({2]) Let I is disconnected graph into the X components T1,T2,---,Ty, then 
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x 
Sipec(T) = U Spec(T). 


Lemma 2.2 Let w =e7 = cos(24=) + isin(#8*) for 1 <k <n, where n is a positive integer 


and i? = —1. Then 
(i) we + wt = 2cos(2#) for 1 <k <n, 


(ii) we —w"* = 2isin(?2@) forl<k <n. 


Lemma 2.3([1]) Let n be a positive integer. Then 


n-1 


= 2k 1 
(i) sin — = <cot—, n = 1(mod2), 
ae 2  2n 
n—2 
2k 
(it) sin ——~ = cot = n = 0(mod2), 
n 
k=1 


2kr 1 T 1 
iii) S cet | = = ese — — =, n = 1(mod2 
(iit) | cos a 5 eae igs tt (mod2), 


n—-1 
2k 
(iv) S- | cos a = 2cot - n = 0(mod4), 
= n n 
eS Qhr T 
(iv) S- | cos ——| = 2csc—, n = 2(mod4), 
—F n n 
2 
. (2k-1)r T 
MOT ~esct n= 2 
(vi) 2 - csc n 0(mod2), 


n—-1 
k 
(vii) s sin — = cot ae n = 1(moa2), 
nm 
k=1 


2n 
eS Qh T 
Ltd —|= —-l1,n=1 d2). 
(vitit) Ne = | = csc rms i (mod2) 


Lemma 2.4 Let n be a positive integer. Then 

oa 
: 2kn T = 

(i) a | cos ie = cot ie 1, n = 0(mod4), 
k=1 
a2 

= 2k 
(it) » | cos =| a ese = —1, n = 2(mod4). 
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Proof The proof of (i) follows directly from Lemma 2.3(iv), and (7) is a consequence of 
Lemma 2.3(v). 














Lemma 2.5 Let n be a positive integer. Then 


n—1 no-1 
(i) S> sin —* = $7 |sin *], n = 1(mod2), 
k=0 k=0 ne 
a 7 
2k; 4k 
(it) sin = | sin a n = 2(mod4), 
k=0 ec” Sipe 4 
pleaded ne 
7 2k; f 2k; 
(ii) S> sin —* = 14.257 sin, n = 0(moad), 
k=1 wy k=1 ie 
a2 nea 
z 4k S 4k 
(iv) S~ |sin —*| = 2 >| sin |, n = 0(mod4), 
k=1 - k=1 a 
a 
4k 4 4 
(v) sin —— = ese — + cot au n = 0(mods), 
kel n nm n 
wee 
a 4k 2 
(vi) S- sin —~ = cot el n = 4(mod8). 
an n n 








2 4 -1 
Proof (i) Let n = 1(mod2), f(k) = sin aE and g(k) = sin le where k € {0,1,2,---, de 
n n 
Then it is easy to check that 
f(2k) f Osk<|"F], 
g(k) = = =1 
—f(n— 2k) if [=| <ks eS 
This implies (2). 
2 4 2 
GN T= ined. Hasta Bad Gay Se RE ON Soe =}. 
n n 


Then it follows that 


f (2k) if 1<k< %, 
—f(-Z+2k) if SB ck< ay. 


g(k) = 


This implies (7). Proofs of (77) and (iv) are similar to that of (i) and (di). 
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(v) Let n = 0(mod8). Then n = 8m, m € N and 


nee 
4 Aka 5 ang vas TE i 30 bio lias (2m — 1)r 
2. sin aS sin =, Le sin nae sin om 
20 _ 4a . (2m —2)r 
+ | sin — + san=— +--+ + s1n ——— 
2m 2m 2m 


= esp +cots 
4 4 

= csc— +cot— (Using Lemma2.3(vi) and (#i)). 
n n 


(vi) Let n = 4(mod8). Then n = 8m+ 4, m € N and 














n—4 
wo. 4k 
sin — - > sin ——— el 
k=1 
2 
= cot ea = cot ao (using Lemma 2.3(vit)) 


Lemma 2.6 Let n be a positive integer. Then 


n—-1 
(i) So jsin 7 = = cot a n = 1(mod2), 
k=1 


n—-1 
(it) ) YJ sin | = 20067, m = o(mod2). 
k=1 


Now we compute skew energy of some Cayley digraphs. 


Theorem 2.7 Let G = {v1 = e,v2,-+- ,Un} be a group, S = {v;} CG with uy, Au; ", vu; Fe 
and T = Cay(G,S) be a Cayley digraph on G with respect to S. Suppose H = (S), |H| =m, 
|G: H| =X. Then 

2rcotx— if m=1(mod2), 


e,([) = 
4rcot= if m=0(mod2). 
Proof Let G = {v1 = €, v2, U3,°++ , Un}, S = {ui}, v; © G with v; re uv, ', vi # e and suppose 
H = (S), |H| =m, |G: H| = 2. If \ = 1 then G = {e,v;, v?,--- , v1} and hence the the 


’ U; 


skew-adjacency matrix of f = Cay(G, S$) is a circulant ue Its first row is [0,1,0,--- ,0,—1]. 
So all eigenvalues of I are Ag = w* — wk-* = wk — wk = 2isin2&®, k =0,1,--- ,n—1 where 
w=em andi? =—1. Applying Lemma 2.2(2i), we obtain Ax = isin 2k, k=0,1,---,n-1. 


Now by Lemma 2.6 we have 


n-1 


2k 2cot+ if n=1(mod2), 
= Finn =e in—| = ie ee) 
Acot= if n= 0(mod2). 
If A > 1, then T is disconnected graph in to the [j,i = 1,---,A components and all 


components are isomorphic with Cayley digraph T,, = Cay(H,S) where H = (uv; : v}” = 1) 
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and m|n, S = {v;}. Since T is not connected, by Lemma 2.1, its energy is the sum of the 


energies of its connected components. Thus 


es(P) =) es (Vi) = Aes(Cay(H, 3)) = ee 


» . 
2rAcots— if m 
j=l 4rAcot= if m 
m 














This completes the proof. 


Theorem 2.8 Let G = {e,b,b?,--. ,b"~1} be a cyclic group of order n, and T = Cay(G, S$) be 
a Cayley digraph on G with respect to S = {b',W)},0<i,j<n—-1,i4j H =(S), |H| =m, 
and |G: H| =X. Then 

(i) e.([) < 4A cot & ifm = 1(mod2), 


(ii) es([T) < 8Acot 5 if m = 2(mod4), 


(iti) e.(T) < 4X\(cot = + 2csc #2 + 2cot ) if m = 0(mods), 





(iv) es([) < 4A(cot = + 2cot 2*)) if m = 4(mod8). 


Proof Let G = {e,b,b?,--- ,b"~1} be a cyclic group of order n and T = Cay(G,$) be a 
Cayley digraph on G with respect to S = {b',b/},0 <i,j <n—-1,i4j, H =(S),|H| =m, and 
|G: H| =. If \=1, then G = A and hence the the skew-adjacency matrix of [ = Cay(G,S) 


is a circulant matrix. So all eigenvalues of I areAy = w* —w7* +w?* —w-?*, k = 0,1,--- ,n—1, 
where w =e and i? = —1. Hence 
2k 2khr . Ak _,. 2kr . Aka 


dp = wh — wk 4 w?* — WH ?* = 2isin— 4+ 2isin—— = 2i(sin— + sin—) 
n n n n 
fork = 0,1,---,n-—1. 


(t) Suppose n = 1(mod2). Then 


= ie Qkr Akr 
eC) = S- ral =) l2i(sin— + sin—)| 
k=0 k=0 
n—-1 
2k Ak 
= \- [2i(sin— + sin—)| 
= n n 
nek 
u 2k Ak 
= A ag! 
mou plots 
a 2k 2 Ak 
< A(S > |sin—*] + $7 |sin—)) 
k=1 ue k=1 ‘ 
cleats net 
a 2kn a 2kn 
a _ 2kT _ 2kT li 
- sin— + sin— ) (using Lemma 2.5(2)) 
k=1 k=1 
as . . 
= 4cot— (applying Lemma 2.3(i)). 


2n 
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Thus (i) holds for \ = 1. 


Suppose A > 1. Then [ is disconnected graph in to the [j,i = 1,--- , A components and all 
components are isomorphic with Cayley digraph T, = Cay(H,S) where H = (uv; : vy’ = 1) and 
mln, S' = {u;}. Since [ is not connected, its energy is the sum of the energies of its connected 


components. Thus 


Tv 
< —. 
= Death i )AEs(Cay(H, S)) < 4A cot Sa 


Now we shall prove (#i),(i7) and (iv) only for A = 1. For \ > 1, proofs are similar to that 
of (i). 


(iz) If n = 2(mod4), then 


n-1 
Ak 
e(T) = 2 el = a (sin +5 sin—)| 


4k 
- ¥ |2i( sin + sin——_)| 
n 


Ak 
- Ds |sin = - sin 


nS 


us sin + |sin=)) 


k=1 


IA 


ws 2h oN ane Teng S08) 
= tn — tn —- usin. mm. -O(2 
—~ * A ~ * A g v 
= 8cot—. 

nm 


Here we used the Lemma 2.3(iz). 


(itt) If n = 0(mod8), then 


nck 
e,([T) = ‘S isin + si —| 
an . 2kr a . Akt 
< 4(5 |sin—| + ) |sin—]) 
k=1 mn k=1 ne 
noe not 
2k 4k 
= A S sin—— +2 s sin——) (using Lemma, 2.5(iv)) 
n n 
k=1 k=1 


4 4 
= Acot~ + 8csc— + 8cot—. 
n n n 


To get the last equality we have used Lemma 2.3(i7) and 2.5(v). 
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(iv) If n = 4(mods8), then 


n-—2 n—-2 

2kr a Akn 
4 ee jay et 
een = Boe = |) 


es(T) < 
noe nes 
2ka 4kr 
= 4 L—— + 2 L—— sing L 2.5(4 
Oo sin— + S sin— ) (using Lemma 2.5(iv)) 


2 
= A(cot— + 2cot—)). 
n n 











To get the last equality we have used Lemma 2.3(i7) and 2.5(vi). 





Lemma 2.9 Let G = (b: b” = 1) be a cyclic group and T = Cay(G, S;), t € {1,--- , [45+] },be 
a Cayley digraph on G with respect to S; = {b!',b?',--. ,b"}, wherel € U(n) ={r:l<r< 
n,gcd(n,r) =1}. Then the eigenvalues of T are 


[Sel Qhjn 
Ae =) Pisin, k=0,1,--: 0-1, 
j=0 





n 


where i? = —1. 


Proof The proof directly follows from the definition of cyclic group and is similar to that 
of Theorem 2.7. 














Lemma 2.10 Let G = (b: b” = 1) be a cyclic group and TV = Cay(G, S;), t € {1,--- , [454] },be 
a Cayley digraph on G with respect to S, = {b!,b7!,--- ,b"} wherel € U(n) ={r:1l<r< 
n,gcd(n,r) = 1}. Also suppose e,((), SLE(T) denote the skew energy and the skew Laplacian 
energy of T’ respectively. Then e,([) = SLE(T). 


Proof The proof directly follows from the definition of the skew energy and the skew 











Laplacian energy. 





Lemma 2.11 Let n be a positive integer. Then 


n-1 n-1 
(i) cos — = Se cos a n = 1(mod2), 
k=l | " 
ns 
= 4k 
(it) S*|cos | = ese 4 — 1, n = 2(moda), 
k=l n n 
a 
Ak 
(iit) cos — = —1, n = 0(mod4). 
nm 
k=l 


2k 4k —1 
Proof (i) Let n = 1(mod2), f(k) = cos ao g(k) = cos —, where k € {1,2,---, e 5 \. 
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It is easy to verify that 


This implies (2). 


(it) Let n = 2(mod4). Then n = 4m + 2 for some m € N. We have 


: | Beene ana 


I 


k 
5 eo |= 5a Too = 


2k 
= 3 | cos Daa - 7 = csc Tee —1 (using Lemma 2.3(viiz) 
= csco—1. 
n 


(iit) Suppose n = 4m,m € N. Then 


n-2 


“2 2m— —-1 2m—-1 
Akan mT 
cos —— = -> cos = = cos — + cost + ) cos 
k=1 k= k=1 k=m+1 


n—2 














= Ak 
Changing k to k +m in the last summation we get S- cos — =—1. 
k=1 
Theorem 2.12 Let G = (b: 6" =1) be a cyclic group and T = Cay(G,S),be a Cayley digraph 
on G with respect to S = {b'} where l € U(n) = {r:1< 1 < n,ged(n,r) = 1} and C,(T) be 
the skew-adjacency matrix of [ , D(T) = diag(di,d2,--- ,dn), the diagonal matrix with vertex 
degrees d,d2,-+- ,dn of e,b,b?,--- ,b"~1. Suppose L(T) = D(T) — C,(T) and p1,--+ , un are 


eigenvalues of L(T). We define a(T) = »; u2. Then 


(4) a(L) < 2n + 2esc= if n = 1(mod2), 
(it) a(T) < 2(n—1) + 4esc* ifn = 2(mod4), 


(itt) a(L) = 2(n — 2) ifn = 4(modo). 


Proof Let G = (b: b” = 1) be a cyclic group and I = Cay(G,S),be a Cayley digraph 
on G with respect to S = {b'} where 1 € U(n) = {r: 1 <r < n,ged(n,r) = 1} and C,(I) 
be the skew-adjacency matrix of T. Note that underlying graph of [I is a 2—regular graph. 
Hence D(T) = diag(2,2,--- ae Suppose L(T) = D(T) — C;,(T) then L(L) is a circulant matrix 
and its first row is [2,—1,0,---,0,1]. This implies that the eigenvalues of L(T) are pp, = 
2— wk + wkr-*§ = 2 — ww = 2 (WR — wr) = 2 — 2isin 242, k = 0,1,---,n—1 where 


2Qri te 
w=en andi? =—-1. 
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If n = 1(mod2), then 


n-1 
al) = Song 
k=0 





n—-1 
- S © (2—2isin—)? =4+ S0(2 2i sin —)? 
k=0 k=1 
nok 
‘ 2k — 
= 4+) °(2-2isin—)? + a (2- di sin") 
k=1 . poet 
A 
a = 
= 4+ 50 -2%sin ty Ses any 
k=1 " 
eee no-1 





k=1 k=1 





n-1 
n-1 ~ 1 1 4k n 
= 4+8 ae Gs cee a 














2 
k=1 
net 1 
n—1 2ka : : 
= 444 5 ) +49 cos (using Lemma 2.11(2)) 
k=1 
< 4+4+4 
-1 1 1 
= A >) + 4(5 ese a — 5) (using Lemma 2.3(i#i)) 
= In+2csc—. 
2n 


If n = 0(mod2), then 


n-1 
al) = Som 
k=0 

















a Ika cae, Qkr 
_ 7 ART \2 gy og an, 2ETy2 
- bme 2i sin rF y=4 So (2 2i sin e ) 
k=0 k=1 
n—-2 
a 2k 9(2 = 2k 
= 4+5°Q- isin —*)? + (2—2isin at + >> Q- isin —*)? 
k=1 kat? 
ne = 
2 2khr 
= 8+ 50 (2-2isin 7 £¥ (04 niin 2 
k=1 k=1 
2k —2 2k 
ze 8+ 2 2(4— 4sin?—*) = 44.8(5 )-8 >> sin? 
k=1 k=1 
oe 
n—2 1 1 Ak 
= 44 8( 5 )-8 (5-908) = 4 


= 2k =| 2 2k 
= 4+ > 2(4— 4sin?—*) = 44.8(— )-8 >> sin? 
nr nr 
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If n = 2(mod4), then 


a(L) 


IA 





n—-2 
n—2 ee Akn 
ca ee ) +4} 7 | cos ——| 
k=1 


n—-1 








4+4(—) + 4(ese —1) (using Lemma2.11(éi)) 
n 


l| 


2(n — 1) +4ese~. 
n 


This completes the proof of (it). 














If n = 0(mod4), then 
nee 
n—2 4k 
af) = 444-5 )+455 cos —— 
k=1 
n—2 ; a 
= 4+4( 5 )+4(-1) (using Lemma2.11(ii7)) 
= 22(n-—2). 


This completes the proof of (i#i). 


Lemma 2.13 Let n 


n 


be a positive integer. Then 


2n 2 


n—1 
(i) » | cos PET = ese 5 +=, n = 3(mod6), 


n—-1 


n-1 


s 6kr os 2kr 
bi 5 —| = — =1 dé 
(it) S- | cos . | > | cos |, n = 1(mod6), 


k=1 


n—-1 


WN 


k=1 


n-1 


z 6kr ra Qkr 
Lit ——/| = — =5 d6 
(iit ys | cos 7 | ys | cos = |, n (mod6), 


k=1 k=1 
a. bk 
(iv) cos — = 0, n = 0(mod6), 
k=1 
He? Hee 
- 6kr = 2kn 
(v) cos — = S- cos n= 2(mod6), 
k=1 k=1 
a 6kr a 2ka 
(vi) cos — = cos —, n = 4(mod6), 
k=1 k=1 . 
a 8kr ae Akan 
vii) cos — = cos ——, n = 1(mod2), 
k=1 k=1 
wz 8kr a Akan 
(vii) cos — = cos —, n = 2(mod4), 
n n 


> 
Il 
un 
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n-2 


g 8k 20 
Lit s —|=-1+2csc—,n=4 d8 
(viii) » | cos . | + 2ese—, 0 (mod8), 
n—2 


= 8k 4 
(ix) S- | cos =| =-—1+4cot ails n = 0(mod16), 
ee n n 
a 
k 4 
(x) S- | cos an =—1+4csc ey n = 8(mod16) and n = 2 or 4(mod6). 
n n 
k=1 


Proof (i) Suppose n = 6m + 3, then 


n-1 




















— 6k ian 2kr 
> | eos ——| a S- [cos od! 
k=1 =1 
2m 3m+1 
2kn 2ka 
= S| cos i} + De | cos 5 ra 
k=1 k=2m+1 
2m m 
2k 2 
= Son BE +S em 
(changing k to k + (2m + 1) in the last summation) 
3. 3 1 
= 50s =~ +5 (using Lemma 2.3(viii), (ii). 
n 


2 
(it) Let n= 6m-+1, f(K) = cos = , g(k) = cos , where k € {1,2,--- ee 


we have 





f (3k) if 1<k<m, 
gk) =% f(n—3k) if m<k<2m, 
f (3k — n) if 2m<k< 3m. 


This implies (#7). 
The proofs of (<i2), (tv), (v), (vt), (viz) and (viiz) are similar to the proof of (iz). 


Suppose n = 4m then 


n—2 


=a 2m-1 
8k 2k 
Y | cos | = S- | cos | 
k=1 k=1 
m—1 2m—-1 
k 2k 
= 1+ S- |cos =| + s | cos —| 
k=1 He k=m+1 
m—1 
2k 
= 142.9 foo BE = 1 +25 feos "| 
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—1+2csc5= if m= 1(mod2) 
(mod4) 


—1+4cotS if m= 


—1+4csc= if m= 2(mod4) 





—14+2csc2* if n= 4(mod8), 
= 4-14+4cot if n=0(modl16), 
—1+44csc #2 if n= 8(modl6). 





This completes the proof of (vidi),(¢x),(). 











Theorem 2.14 Let G = (b: 6” =1) be a cyclic group and T = Cay(G,S),be a Cayley digraph 
on G with respect to S = {b',b*} where l € U(n) = {r:1 <r <n,ged(n,r) =1} and C,(T) be 
the skew-adjacency matrix of [ , D(T) = diag(di, d2,--- ,dn), the diagonal matrix with vertex 
degrees d1,d2,-+- ,dn of e,b,b?,--- ,b"~1. Suppose L(T) = D(T) — C,(T) and p1,--+ , Mn are 


eigenvalues of L(T). Define a(T) = S- ps2. Then 
i=1 


(i) a(L) < 4(8n + 2) — 12 csc 32 if =n = 3(mod6). 





(ii) a(L) < 4(8n+2)—-4ese35 if n=1 or 5(mod6). 





(iit) A(T) < 4(8n—2)+ 16 csc if n= 2(mod4) and n = 0(mod6). 

(iv) a(T) < 4(8n—2)+24cesc5 if n= 2(mod4)and n= 2 or 4(mod6). 

(v) a(f) < 4(8n — 2) + 8cot 2+ 8csc2= if n= 4(mod8) and n = 0(mod6). 

(vi) a(L) < 4(8n — 4) + 16 cot 2+ 8cse% if n=4(mod8) and n= 2 or 4(mod6). 
(vit) a(T) < 4(8n — 2) + 8cot 2+ 16 cot = if n=0(mod16), and n = 0(mod6). 





(viii) a(L) 


IA 


2(n — 8) + 16 cot 2 + 16cot4= if n=0(mod16) and n= 2 or 4(mod6). 





TT 
n 


(ix) a(T) < 4(8n — 2) + 8cot 2+ 16csc4= if n= 8(mod16) and n = 0(mod6). 





(x) a(P) < 4(38n — 4) + 16 cot 2+ 16csc4= if n= 8(mod16) and n= 2 or 4(mod6). 

Proof Let G = (b: b” = 1) bea cyclic group and T = Cay(G, $),be a Cayley digraph on 
G with respect to 9 = {b!,b™} where 1 € U(n) = {r:1<r<n,gcd(n,r) = 1} and C,(I) be 
the skew-adjacency matrix of I. Note that underlying graph of I is a 4—regular graph. Hence 
D(L) = diag(4,4,--- ,4). Suppose L(T) = D(T) — C,(T) then L(L) is circulant matrix and its 
first row is [4,—1,—1,--- ,0,1,1]. This implies that the eigenvalues of L(T) are 


2k Ak 
pe =4—w* wk 4 2k 4 wok = 4 2i(sin = + sin —), k=0,1,---,n-1, 
n n 





Qri 


where w =e 





and 7? = —1. It is clear that 


_, 24kr _ Aka 
bin—k = 4+ 2i(sin — + sin ae) and pp + Un—~- = 8 
n 
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n—-1 





for k = 1,2,---,| |. So 
2 2 = 
Ma + Mn—p = 64— 2ukbn—K 

2k 4k 2k 4 

5 Edi din SO a ein ey 
nm nr nr n 
2 Ak 

= 64—2(16 + 4(sin an + sin —)) 





2k 4k 6k 8k 
= 9d Reo Laebg = 8 toa eg 
n nr nr n 
n—-1 
fork =1,2,---,| 5 |. Let n = 1(mod2), then 
n-1 
a(t) =. SP pe 
k=0 


= 164+ 5 (uk + ue_x) 


n-1 


- 2k 4k 6k 8k 
= 16 + $7 (24— 8cos —* + 4cos ——~ — 8cos — + 4cos ——) 
nm nm n nm 


k= 








n=l n-1 n=l n=l 
eae 2kr =e Akt <ae 6kr en 8kr 
= A(3 1)-8 s— +4 —-8 s— +4 —_ 
(3n + 1) a cos — + », cos — », cos — + », cos — 


n—1 n—1 n—1 n—1 
<m 2kr ee 2Qkhr nm 6kr a, Qkr 
= 4(3n+1)-8 ae a8 aed patie 
(3n + 1) », COs ; + 2, cos - », cos . - 2, cos 7 


(using Lemma2.11(2), 2.13(vit)) 


n-1 n—-1 


a 6kr = 6kr 
= 4(8n+1)-8 S cos — < 4(3n+1)-8 s | cos —]. (2.1) 
n n 
k=l k=1 


(i) If n = 3(mod6) then using Lemma 2.13(i) in above inequelity, we get 


3. 30 1 3 
a(L) < 4(3n +1) 8(5 ose — 5) =4(8n + 2) — 12ese 


This completes the proof of (i). 


(it) If n = 1 or 5(mod6) and using Lemma 2.13(i2) and (cit), we get 


a(T) 


IA 


rc 
nm 


n—1 

=s ok 
4(3n + 1)-8) 0 | cos e 

k=1 





1 T 1 T 
= 4(3n + 1) — 8(, ese 5 — 3) = 4(3n + 2) — 4csc an: 
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Let n = 0(mod2). Then 


n-2 


n-1 n-1 > 
aT) = Domp=wp+ue+ D> we=16+164+ >) (oh + u2_«) 
k=0 k=1,k¢® k=1 
n-—2 
zs 2k 4k 6k 8k 
= 32+ S > (24- 8.cos fdieos asagaee + 4cos ——) 
hel nr nr nr nr 


n-2 n-2 n—-2 n-2 

= 2kr Aka e 6kr “ 8kr 
= 43 2)-—8 —+4 5 — — 8 — +4 S—. (2.2 
(3n + 2) » cos — + 2, cos — » cos — + 2, cos — (2.2) 


If n = 2(mod 4), then employing Lemma 2.13(viii) in (2.2), we get 








n—2 n-2 n-2 n-2 
= 2kr —_ Akan 6kr ae 4kn 
a(T) = 4(3n +2) — 8) | cos —— gO Do eres — 8) | cos —— +4 cos 
k=1 k=1 k=1 k=1 
ae oa a 
2ka Akan 6kr 
= PREIS ye DL (2.3) 


(iti) If mn = 2(mod4) and n = O(mod 6), then using Lemma 2.13(iv) in (2.3) we deduce 
that 


a(T) 


IA 


n~2 n-2 
2Qhr a Akn 
A(3 2)+8 y s——|+8 y — 
(3n + 2) Ze = | + Zee 5 | 





= A(3n+2)+ 16(cse ~1) = 4(3n — 2) +16 cso“ 


by using Lemma 2.4(iz) and 2.11(#i). 


(iv) If n = 2(mod4) and n = 2 or 4(mod6), then using Lemma 2.13(v) and (vi) in (2.3) we 
see that 


n—2 n—2 n—2 
2ka = 4kr = 2kr 
a(T) = 4(8n+2)-8 > est 8 » cog 8 » 8 


IA 


= 


n~=2 n-2 

am 2Qhr —s Akr 
4(3n +2) +1 a : 
(3n ) ome " prea ee 








< 4(3n +2) + 24(cse— — 1) = 4(3n — 2) + 24 ese =. 
nr n 











Similarly we can prove (v) to (x). 





We give few interesting results on the skew energy of Cayley digraphs on dihedral groups 
Don. 


Theorem 2.15 Let Dan = (a, bla? = b” = 1,a~'ba = b') the dihedral group of order 2n and 
T = Cay(Doan,S) be a Cayley digraph on Don with respect to S = {b°}, 1 <i<n-—1, and 
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H =(S), |H| =m, |Dj,, : H| =A that,Dj,, is the commutator subgroup of Dyn. Then 


4dr\cots— if m=1(mod2), 


e,([) = 
8Acot= if m=0(mod2). 


Proof The proof of Theorem 2.15 directly follows from the definition of dihedral group 
and Theorem 2.7. 














Theorem 2.16 Let Dan = (a, bla? = b" = 1,a~'ba = b~') the dihedral group of order 2n and 
T = Cay(Dan, $) be a Cayley digraph on Don, with respect to S = {b',b}, 1 <i,j<n-li¥g, 
and H = (S), |H| =m, |D5, : H| = A Then T = Cay(Don,S) is a circulant digraph and its 


skew energy 

(i) e.([) < 8Acot ifm = 1(mod2), 

(ii) es([) < 16Acot= if m = 2(mod4), 
(iti) e,() < 8\(cot= + 2csc + 2cot42) if m = 0(mods), 
(iv) e.() < 8A(cot= + 2cot2®)) if m = 4(mod8). 


Proof The proof of Theorem 2.16 directly follows from the definition of dihedral group 
and Theorem 2.8. 














Theorem 2.17 Let Dan = (a, bla? = b" = 1,a~'ba = b') the dihedral group of order 2n and 
T = Cay(Dan, S) be a Cayley digraph on Don with respect to S = {b'} where | € U(n) = {r: 
1<r<_n,gcd(n,r) = 1} and C,(T) be the skew-adjacency matriz of T, D(T) is then xn 
matrix such that di; = 2 whenever i = j otherwise dj; = 0. Suppose L(T) = D(T) — C,(T) and 


A1,++:;An are eigenvalues of L(T). Define a(T) = Se d?. Then 
i=1 
(i) a(f) < 4n+4cesc $= ifn = 1(mod 2), 
(ii) a(L) < 4(n—1) 4+ 8esc% ifn = 2(mod 4) 


(iti) a([) = 4(n — 2) ifn = O(mod 4) 


Proof The proof of Theorem 2.17 directly follows from the definition of dihedral group 
and Theorem 2.12. 














Theorem 2.18 Let Dan = (a, bla? = b" = 1,a~'ba = b~') the dihedral group of order 2n and 
T = Cay(Don, S$) be a Cayley digraph on Don with respect to S = {b',b7} where 1 € U(n) = 
{r:l<r<n,gcd(n,r) = 1} and C,(T) be the skew-adjacency matrix of T , D(L) is thenxn 
matrix such that di; = 4 whenever i = j otherwise dj; = 0. Suppose L(T) = D(T) — C,(T) and 


A1,++: ;An are eigenvalues of L(T). Define a(T) = S- d?. Then 
i=1 
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(i) a(T) < 8(3n + 2) — 24 csc 3 if nm =3(mod6). 





(ii) a(T) < 8(8n+2)—8esce$ if n=1 or 5(mod6). 


(iii) a(T) < 8(3n — 2) 4 
(iv) a(L) < 8(3n — 2) + 

(v) a(L) < 8(8n — 2) 4 
(vi) a(L) < 8(3n — 4) 4 
(vii) a(P) < 8(3n — 2) 4 


| 32 csc = 
+ 48 csc = 
+ 16 cot = + 16 csc 


+ 32. cot = + 16 csc 





(viié) a(L) 


(ix) 


IA 


a(L) < 8(3n — 2) 4 


+ 16 cot = + 32 cot 


A(n — 8) + 32cot 2 + 32cot 4 


+ 16 cot = + 32 csc 





(x) a(Lf) < 8(38n — 4) 4 


+ 32 cot = + 32 csc 


if n=2(mod4) and n = 0(mod6). 


if n= 2(mod4)and n = 2 or 4(mod6). 


27 


n 


if n=4(mod8) and n = 0(mod6). 


27 


n 


if n =A4(mod8) and n = 2 or 4(mod6). 
42 if n=0(mod16), and n = 0(mod6). 


if n=O0(modl16) and n = 2 or 4(mod6). 


An 


n 


if n= 8(mod16) and n = 0(mod6). 


An 


n 


if n =8(mod16) and n = 2 or 4(mod6). 


Proof The proof of Theorem 2.18 directly follows from the definition of dihedral group 


and Theorem 2.14. 
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Abstract: A change of Finsler metric L(x,y) — L(x,y) is called Kropina change if 
2 


a L 
L(a,y) = FB? where B(x,y) = bi(x)y’ is a one-form on a smooth manifold M". The 
change L — L is called projective change if every geodesic of one space is transformed to a 
geodesic of the other. The purpose of the present paper is to find the necessary and sufficient 


condition under which a Kropina change becomes a projective change. 
Key Words: Kropina change, projective change, Finsler space. 


AMS(2010): 53060, 53B40 


§1. Preliminaries 


Let F” = (M",L) be a Finsler space equipped with the fundamental function L(x, y) on the 
2 


; L 
smooth manifold M”. Let 3 = b;(x)y’ be a one-form on the manifold M”, then L > a is called 
a AEA =" = 

Kropina change of Finsler metric [5]. If we write L = a and F = (M",L), then the Finsler 
space F’’ is said to be obtained from F” by Kropina change. The quantities corresponding to 
F’” are denoted by putting bar on those quantities. 

The fundamental metric tensor g;;, the normalized element of support /; and angular metric 
tensor h;; of F' are given by 

1 0?L? OL OL 


a 5 Oy*Oy) Oy’ i i Oy’ dys ae 





= ils. 


We shall denote the partial derivative with respect to x’ and y’ by 0; and 9; respectively 


and write 
Lp=OL, Liz =0,0;L, Ligh = 0,0j0,L. 
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Thus 
L,= hk, Lot hij _ Li; 


The geodesic of F’” are given by the system of differential equations 


2 ani 
d*x +26! ( =) =o. 


ds? ds 
where G"(z, y) are positively homogeneous of degree two in y’ and is given by 
: oe : L? 
2G" = g" (y"0;0,F — 0;F), F= i 
where g’! are the inverse of gj;. 
The well known Berwald connection BI = {G%,,,Gj} of a Finsler space is constructed from 
the quantity G’ appearing in the equation of geodesic and is given by [6] 


Gi=d;G', Gi, = OnG4. 


The Cartan’s connection CT = {F? 


“nr Gi, C%,} is constructed from the metric function L 


by the following five axioms [6]: 

(4) gigje = 0; (44) gijle = 0; (itt) FY, = Fy; (iv) FQ, = Gi; (v) Ch, = Ci. 
where \k and |, denote h and v—covariant derivatives with respect to CT. It is clear that 
the h—covariant derivative of L with respect to BI and CT are same and vanishes identically. 
Furthermore the h—covariant derivatives of L;, Li; with respect to CT are also zero. 

We denote 


§2. Kropina Change of Finsler Metric 
The Kropina change of Finsler metric L is given by 

— | : 
(2.1) L= om where B(x, y) = bi(x) y’. 
We may put 
(2.2) G =Gi+ Dt 


Then G; = Gi + Dj and GC = Gi, + Di,, where Di = 0; D* and Di. = OpDi. The tensors 
D', Dj and Dj, are positively homogeneous in y’ of degree two, one and zero respectively. 

To find D* we deal with equations Lijjk = 9 [2], where L,,), is the h—covariant derivative 
of Lig = hij /L with respect to Cartan’s connection CT. Then 


(2.3) OnLiz — LigrGy, — Lrj Fi, — Liv Fj, = 0. 


Since 0;3 = b;, from (2.1), we have 
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_ 2L L? 
= 2L 2 _ 2L 4, 2b? 
2 - QE 2 21? 2L L? 
(c) O,.L; = Fy (nbs) + (Gu-= a ss OnE + .(e B by — re) ts) On 3 — RP =; (Onb;); 
CL F (Osby) $ bos Fp(Labj + Lybi) + x bt (nL) 
2 
- | Fit +3 eon +o bibs — Fa (Lib, +L; iby] (4,3) 
2 2L 
+ [Be ee [ju Be | Ort) 
2 2 
arb , as (Dubs) + are -SL | (xj) 
(e) Lijk = a + List + Lj Lin + Ly Liz) — (Labs + Lind; + L5xb:) 
—FlLib; be Teed, jabs) + 5 (Oy Eek bibs 
is bibjbe- 


Since Lj), = 0 in F"’, after using (2.2), we have 
O,Euy — Lige(G. + DE) — Teg(Fh, + °Dh,) — Eee(Fhe + “Dh = 0, 
where F’,, — Fi, = °Di,. 
Using equations (2.3) and (2.4)(b), (d), (e), the above oe may be written as 


2L 2 
Fae BR 


a) 
4L a L 2 6L? 4L 
+ So.)| L,Gi, — [Fly + eliL, + Bi — bibs — Be (bibs y+ Lb] » 


[LijrDy, + Levy “Dig + Liv ©Djy] + — (Lib; + L;bi) 


2 21 2 2L 
Tok + Sok + b, Gi, + ( L; ——b; i) LirGy, + Ly "A + (Su = zh) x 
( Ok k k) a ge? ( k ie) B G2 
217 217 


“pe 9 — 15) (Tik + Sik + b,-F i, ) 


(LirGy + L,. Ha) + ( 














217 21 2L 
(2.5) a0 (Fro, a ats) (rik + 85k + bp FR) + {3 5 (Lijbr + Lirds + Ljrbi) 
B B B 
2 2 
a as OS PN ee Oa oe a LS Ae FO SFP 
: Fa (bible 4 bybp Ly + bibp Ly) + + Oe ide hie + Df) 
oL 2 21? . 
+ {Ebaby + bib) - Gnd; arias} P+ cpt.) 
QL 2, _ 2b ae 
+ {Fe(dab, + L,b;) — 3 Like — bbe (Fh + ik) = 9, 
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where ‘0’ stands for contraction with respect to y’ viz. rox = Tiny’, Too = Tijy'y?- 
Contracting (2.5) with y*, we get 


oh 2L 2 
2 {bin — SS (Ligbp + Lipby + Lyrbi) — Sp (LiLgby + LiLpbj + Ly Lyb;) 





B (? & 
2 
+ 5 (Libin + Ljbin + Leb) + Fy (bibgbe + bj bp Li + bibl) — E-bbjby LD” 
2L 2 6L? 4D 
ape {2 ae hij ary eS LL; + “ge a BB — (L; b; + L; b; )} 10 
(2.6) ; x 
ale Ty (Lr by ee ge ib; Dy 
2 
+ Lip — a” (ibs + Db; Ri Pe ae = bib D} 
F ) B Be 
2L se 2L 2L? 
+ (Ss — , (rio + 5:0) + ( = =) (rj0 + 850) = 0, 
where we have used the fact that Di,y? = °D‘,y? = Dj, [3]. 


>a 


Next, we deal with Li; = = 0, that is 0; y= LirG; —L, F;, = 0. Then 
(2:7) 0;L; — Lir(G5 + D5) — Ly (FG + °Di;) = 0. 
Putting the values of 0;L;, Li, and L,. from (2.4) in (2.7) and using equation 


we get 
L? OF 2 OL 
~ a Pals = oe + Bribe i Be (bibr + Ly-bj) + a ib |B; 
2L QL? aL L? ae. 
+ (Fa = ah) (r50 + 850) + Faz =— S,| Dis; 


where Dik — OnD; —b ri. 
Since 2rij = by; + bj );, 28i7 = bi; — bj;, the above equation ee 





2L? EF 2 2L wo] Ds 
(2.8) ee B B ~ ga ee ae - 
; 2L y) oF, 
+ Et. + go — Be Libr + L,b;) + = b | Br 
2L 2L? 2L 2L? 
+ (Gabe ~ Fart) (oa) + (Get — et) (ro +0) 
2L L? | 
+2|—L,- —b “Di; 
Foe 
and 
2L? 2L 2 2L 
=F = | Ein + bib Be (Libr + Debi) + 7 bir] Dy y 
2L 2 2L 
(2.9) Fea + Ghilr— (Libr + Lrbj) + ae Dy 
2L 21? 
(Fe - rh) (r50 + 850) = (Fa - by) ( (rio + Si0)- 
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Subtracting (2.8) from (2.6) and re-arranging the terms, we get 

















21 21 2 
ta = Be (Peale t Lib; Ljb;) ga bikie or LiL ,b; + L;L,b;) 
— et Lilie Delig y+ 7 car b; Lp + bjbp Ly + bibp Lj) — S bib; bp bor 
(2.10) 
, 3b? 21 i? 
LL; + a — bib; — B — (L; b; a L;b;) Too — Be 
it he a Dh. 
Contracting (2.10) by y*, we obtain 
2L 21 2 21? 
re Lip + (Libr + Lybj) — GL j Ly — arbi D" 
ii) BB B B 
ee L? oD i? “ 
Be ~ gat Too — B25 = ge ger —b,. D, .. 


Subtracting (2.9) from (2.6) and re-arranging the terms, we get 








2L 2L 
{Fb — Fp (Lijbr + Liss + Ljebi) — (ily Eo mene ON ee ae 
2 AL 6L 
+ Btilir + Lj Lir + L, Li) + Be (bibs be + bj bp Li + bb, L;) — BF ol", b; b, 
1 31? 2L 
(2.12) +|3 aby + ay bibey + Fp bib ~ Fy (Libs + Ly ib] on 
2L 2L? 2L 2L 2 
+ (Fe. = =) ee bon — Fy Liby + Lnbi) + GLb 
ice L? 
ag “bib DS j ~~ gr ets: 
Contracting (2.11) and (2.12) by y’, we get 
2L L? L? 
2.13 oo eee ae a 5 eee 
ae, | BB | age’ 
and 
2L 2L 2 L? 
(2.14) bon — B — (L; b, + Ly bj ) + Boe ea 7 ies b .| Dr = = 3 Ee + (Fu - 5 B 1] Too: 


In view of LLi, = gir — L;L,, the equation (2.14) can be written as 


2L 


2 2 
(2.15) Fond? + [FE a (6,D") - Fyb(L-D") = - Fs + (Sb 35 1s) roo 


B BB 


Contracting (2.15) by b’ = gb;, we get 





(2.16) 26° L?(b,D") — 2b°BL(L,D") = —BL? 89 + (L7b? — 8 )roo, 
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where we have written so for 5, 9b". 


Equation (2.13) can be written as 
(2.17) —2L7(b,D") + 4BL(L,D") = —L? roo. 


The equation (2.16) and (2.17) constitute the system of algebraic equations in L,D" and 
b,D". Solving equations (2.16) and (2.17) for b,D" and L,D", we get 


1 
(2.18) b= apple” — 287)rog — 2GL7 89] 
and 

1 
(2.19) Lae — spas lE*s0 + BL rool. 


Contracting (2.15) by g?* and re-arranging terms, we obtain 
. 26L(b,D") — BL L? 26L(L,D") — 2L?(b,D" nae i 
(2.20) Dia BL( ) = BLroo yi 4 |Z too+ BL( ) ( ) 
232 232 23 
Putting the values of b,D" and L,D" from equations (2.18) and (2.19) respectively in (2.20), 
we get 


i (Broo+L?s50\ ,; Brot L*so\ ;  L? ; i_y 
(2.21) D= a b’ — a) a Yy ~ 3B°0 where |’ = Tr 


Proposition 2.1 The difference tensor D* = G -G of Kropina change of Finsler metric is 
given by (2.21). 


§3. Projective Change of Finsler Metric 


The Finsler space F’”’ is said to be projective to Finsler space F'” if every geodesic of F” is 
transformed to a geodesic of F’. Thus the change L — L is projective if G =Gi+ P(x,y)y’*, 
where P(zx,y) is a homogeneous scalar function of degree one in y’, called projective factor [4]. 

Thus from (2.2) it follows that L — LT is projective iff D' = Py’. Now we consider that 


the Kropina change L > L = a is projective. Then from equation (2.21), we have 
(3.1) Pyi = Broo + L? 89 pi Broo t+ L?sp\ ; L? $ 
eee 2623 B02 Yo ~ 9g %0° 


Contracting (3.1) by y; (= gi; y’) and using the fact that s}y; = 0 and y; y’ = L?, we get 


(3.2) io (Broo + L750). 


1 
2b? L? 
Putting the value of P from (3.2) in (3.1), we get 


(3.3) Broo t+ L?so\ « _ (Broo + L*s80 bi — LD? i 
ey 228 2B °0° 
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Transvecting (3.3) by b;, we get 


2 
(3.4) Too = fe, where A= (5) —vP £0. 


Substituting the value of roo from (3.4) in (3.2), we get 


1 
(3.5) P= aA °° 


Eliminating P and roo from (3.5), (3.4) and (3.2), we get 


(3.6) si = Fa = | = 


The equations (3.4) and (3.6) give the necessary conditions under which a Kropina change 
becomes a projective change. 
Conversely, if conditions (3.4) and (3.6) are satisfied, then putting these conditions in 


(2.21), we get 


Di= Sy ie. D'= Py’, where P= =. 


Thus F’’ is projective to F”. 


Theorem 3.1 The Kropina change of a Finsler space is projective if and only if (3.4) and 
2 
(3.6) hold and then the projective factor P is given by P = 5%, where A = (4) — b. 


§4. A Particular Case 


Let us assume that L is a metric of a Riemannian space ie. L = \/aij(x)y’y? = a. Then 


2 
Bi which is the metric of Kropina space. In this case bj; = b;,; where ;7 denotes the 


covariant derivative with respect to Christoffel symbols constructed from Riemannian metric a. 


Thus r;; and s;; are functions of coordinates only and in view of theorem (3.1) it follows that 


the Riemannian space is projective to Kropina space iff ro9 = —F 29 andes, = (Sv - i) 5 
a 
B 


2 
where A = (£) —b? £0. These equations may be written as 
a 


(4.1) (a) ro08? = a?(b?ro9 — B80); (b) 8§(8? — b? a”) = (G?y' — ab") 80. 


From (4.1)(a), it follows that if a? 4 o(mod £) i.e. 3 is not a factor of a7, then there exists 
a scalar function f(x) such that 


(4.2) (a) b?r00 — B80 = 6? f(x); (b) roo = a? f(z). 
From (4.2)(b), we get rj; = f(x)aij and therefore (4.2)(a) reduces to 


Bso = (b’a* — 3) f(x). 
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This equation may be written as 
bis3 + b 58% = 2(b7ai; _ bib;) f(x) 


which after contraction with b’ gives b*s; = 0. If b? 40 then we get s; = 0, ie. s;; =0. 
Hence equation (4.1) holds identically and (4.2)(a) and (b) give 


(ba? — 6?)f(z)=0 ie. f(z) =0 as Wa? — 6 £0. 


Thus roo0 = 0, i.e. r4; = 0. Hence b;,; = 0, ie. the pair (a, 3) is parallel pair. 
Conversely, if 6;,; = 0, the equation (4.1)(a) and (4.1)(b) hold identically. Thus we get the 
following theorem which has been proved in [1], [7]. 


Theorem 4.1 The Riemannian space with metric a is projective to a Kropina space with 
2 


metric a iff the (a, 8) ts parallel pair. 
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Abstract: A function f is called a geometric mean labeling of a graph G(V,£) if f : 
V(G) — {1,2,3,...,q+1} is injective and the induced function f* : E(G) — {1, 2,3,...,q} 
defined as 

f* (wv) = | VFWFO)|, Vuv € E(G), 
is bijective. A graph that admits a geometric mean labeling is called a geometric mean 
graph. In this paper, we have discussed the geometric meanness of graphs obtained from 


some graph operations. 
Key Words: Labeling, geometric mean labeling, geometric mean graph. 
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§1. Introduction 


Throughout this paper, by a graph we mean a finite, undirected and simple graph. Let G(V, EF) 
be a graph with p vertices and q edges. For notations and terminology, we follow [3]. For a 
detailed survey on graph labeling, we refer [2]. 

Cycle on n vertices is denoted by C;, and a path on n vertices is denoted by P,. A tree T 
is a connected acyclic graph. Square of a graph G, denoted by G?, has the vertex set as in G 
and two vertices are adjacent in G? if they are at a distance either 1 or 2 apart in G. A graph 
obtained from a path of length m by replacing each edge by C’;, is called as mC’,-snake, for 
m>ladn> 3. 

The total graph T(G) of a graph G is the graph whose vertex set is V(G) U E(G) and two 
vertices are adjacent if and only if either they are adjacent vertices of G or adjacent edges of 
G or one is a vertex of G and the other one is an edge incident on it. The graph Tadpoles 
T(n,k) is obtained by identifying a vertex of the cycle C,, to an end vertex of the path P;,. The 
H-graph is obtained from two paths uj, u2,...,Un and v1, V2,°++ , Un of equal length by joining 
an edge u nti Ungi when n is odd and u nt2Un when n is even. An arbitrary supersubdivision 
P(m,mz2,--* ,Mn—1) of a path P, is a graph obtained by replacing each i*” edge of P, by 
identifying its end vertices of the edge with a partition of Kom, having 2 elements, where m; is 
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any positive integer. G © Ky is the graph obtained from G by attaching a new pendant vertex 
to each vertex of G. 


The study of graceful graphs and graceful labeling methods was first introduced by Rosa 
[5]. The concept of mean labeling was first introduced by S.Somasundaram and R.Ponraj [6] 
and it was developed in [4,7]. S.K.Vaidya et al. [11] have discussed the mean labeling in the 
context of path union of cycle and the arbitrary supersubdivision of the path P,,. S.K.Vaidya 
et al. [8-10] have discussed the mean labeling in the context of some graph operations. In [1], 
A.Durai Baskar et al. introduced geometric mean labeling of graph. 


A function f is called a geometric mean labeling of a graph G(V,E) if f : V(G) - 
{1,2,3,---,q@+1} is injective and the induced function f* : E(G) > {1,2,3,---,q} defined as 


fur) =| VF) FQe)|, vuv € BG), 


is bijective. A graph that admits a geometric mean labeling is called a geometric mean graph. 


In this paper we have obtained the geometric meanness of the graphs, union of two cycles 
Cm and C,, union of the cycle C,, and a path P,,, Pe. mC;,-snake for m > 1 and n > 8, the 
total graph T(P,,) of P,, the Tadpoles T(n,k), the graph obtained by identifying a vertex of 
any two cycles C,, and C,, the graph obtained by identifying an edge of any two cycles Cy, 
and C;,, the graph obtained by joining any two cycles C;, and C, by a path Py, the H-graph 
and the arbitrary supersubdivision of a path P(1,2,---,n—1). 


§2. Main Results 


Theorem 2.1 Union of any two cycles Cy, and Cy is a geometric mean graph. 


Proof Let ui, U2,-++ ,Um and v1, vV2,°-+,Un be the vertices of the cycles Cy, and C), re- 
spectively. We define f :V(Cm UC,) > {1,2,3,---,m+n-+ 1} as follows: 


i if 1<i<|vm+2|-1 
ee i+1 if |/m+2| <i<m-l1, 
f(um) =m-+ 2 and 

m+n+3-2i if 1<i< || 
f(vi)=4 m+1 if i=|S)+1 
m—n+t 2 if |B)/+2<i<n. 





The induced edge labeling is as follows: 


a if 


1 
Putin) = 9 ; 
i+1 if | m+2|<i<m-l, 
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f*(u, Um) = [Vm + 3 


mtn+1—-2i if 1<i< [2 





. m+1 if ¢=|]24]+1 and nis odd 

FP (vivigi) = L2 
m+2 if i= |] +1 and nis even 
m—n-+ 2i if |B) +2<i<n-1 

and f*(v1vn) =m +n. 


Hence, f is a geometric mean labeling of the graph Cy, U C;,,. Thus the graph Cy, U Cy is 
a geometric mean graph, for any m,n > 3. 














A geometric mean labeling of C7 U Cio is shown in Fig.1. 





Fig.1 


The graph C,, UnT,n > 2 cannot be a geometric mean graph. But the graph C,, UT may 
be a geometric mean graph. 


Theorem 2.2 The graph Ci, U Pp is a geometric mean graph. 


Proof Let uy, u2,-++ ,Um and v1, V2,°++ ,Un be the vertices of the cycle Cy, and the path 
P,, respectively. We define f : V(Cm U Pr) > {1,2,3,---,m+n} as follows: 


m+n+2-21 if 1<i 
f(w)=< n if 7 

n—-m—1+2i if 
f(u) =1, for 1 <i<n-—1 and 
f(vn) =n+1. 
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The induced edge labeling is as follows: 


m+n—2i if 1<is || 
. n if i= +1 and ™ is odd 
f* (uitiga) = 





| +1 and m is even 


2<i<m-l, 


n+1 ce | 
n—m—1+42i if || 

2 
f*(uitm) =m+n-—1 and 


f* (vivi41) = 1, for 1 < 1 < n—1. 


Hence, f is a geometric mean labeling of the graph C,, U P,. Thus the graph Cy, U Py is a 














geometric mean graph, for any m > 3 and n > 2. 


A geometric mean labeling of C12 U P; is shown in Fig.2. 





igis Y ie 

1 

V2 2 

15 2 

14 13 U3 @ 
149U10 ua 413 U4 i 
1 11 ie bb 
12 ul 5 
U6 6 

6 

8 7 9 v7 $8 


Fig.2 


The T-graph T,, is obtained by attaching a pendant vertex to a neighbor of the pendant 


vertex of a path on (n — 1) vertices. 


Theorem 2.3 For a T-graph Ty, Ty U Cm is a geometric mean graph, forn > 2 andm > 38. 


Proof Let ui, u2,++: ,Un—1 be the vertices of the path P,_; and u,, be the pendant vertex 
identified with u2. Let v1, v2,--- ,Um be the vertices of the cycle Cy. 


V(Tn U Cm) = V(Cm) UV (Pn) U {un} and 
E(Ty, U Cm) = E(Cm) U E(Py) U {tin}. 


We define f : V(T, UCm) > {1, 2,3,---,m-+n} as follows: 


f(u;) =t4+1, for l<i<n-2, 
f (Un-1) =n-—l1, 
flun) =1, 
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m+n+2—-2% if 1<is|4| 
Chetan if i=(P|41 
pei a ae |S] +2<i<m. 


The induced edge labeling is as follows: 


f* (uiuigi) =i+1, for l<i<n-2, 





f* (u2tn) =1, 
m+n—2i if isis || 
n if i=|>| +1 and ™ is odd 
f* (vivi4i) = and 
; m 
n+1 if i= |] +1 and mis even 


n—-m—1+26 if |=] 2<i<m-1 
f* (ium) =m+n-—1. 


Hence f is a geometric mean labeling of T,,UC),. Thus the graph T;, UC, is a geometric mean 











graph, for n > 2 and m > 3. 





A geometric mean labeling of T7 U C¢ is as shown in Fig.3. 


8 





Fig.3 


Theorem 2.4 P? is a geometric mean graph, for n > 3. 


Proof Let v1, V2,+++ ,Un be the vertices of the path P,,. We define f : V(P?) — {1,2,3,--+ ,2(n— 
1)} as follows: 


f(u) = 2i-1, for 1 <i<n-Jand 


f(un) = 2(n— 1). 
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The induced edge labeling is as follows: 


f* (vivigi) = 2i-1, for 1 <i<n-—1 and 
f* (vivige) = 2i, for 1<i<n-2. 





Hence, f is a geometric mean labeling of the graph P?. Thus the graph P? is a geometric mean 











graph, for n > 3. 





A geometric mean labeling of P} is shown in Fig.4. 





Fig.4 


Theorem 2.5 mC,,-snake is a geometric mean graph, for anym > 1 andn = 3,4. 


Proof The proof is divided into two cases. 
Casel n=3. 

Let vl) y) and ow) be the vertices of the i” copy of the cycle C3, for 1 < i < m. 
The mC%3-snake G is obtained by identifying ol) and vith) for 1 <i < m—1. We define 
f : V(G) = {1,2,3--+ ,3m + 1} as follows: 

f(v) = 31-2, forl<i<m 
f (vs?) = 31, for 1 <i< mand 
fo) = 314+ 1,forl<i<m. 
The induced edge labeling is as follows: 
fr (vO) = 31-2, forl<i<m, 
fr(v ®) = 3i, for 1<i<m and 
fi(vv) = 31-1, forl <i<m. 


Hence, f is a geometric mean labeling of the graph mC’3-snake. For example, a geometric mean 


labeling of 6C’'3-snake is shown in Fig.5. 
3 6 9 12 15 18 
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Case2 n=4. 


Let vl), us), ol and ul) be the vertices of the i” copy of the cycle C4, for 1 < i < m. 
The mC4-snake G is obtained by identifying yw) and vith) for 1 <i < m—1. We define 
f : V(G) = {1,2,3,---,4m + 1} as follows: 


f(vy’) =4¢-3, for 1<i<m, 
f (vs?) =4i-1, forl<i<m, 
fvo®) =4i, for 1 <i< mand 
flv) =4i+1,forl<i<m. 


The induced edge labeling is as follows: 


fv) = 4i-—3, for 1 <i<m, 
fr (vw) =4i-1, forl<i<m 
fr (vo) = 4i, for 1 <i< mand 
fv) = 41-2, forl<i<m. 














Hence, f is a geometric mean labeling of the graph mC4-snake. 
A geometric mean labeling of 5C4-snake is shown in Fig.6. 


334 ava: 11472 15 15 16 19 19 20 





Fig.6 


Theorem 2.6 T(P,,) is a geometric mean graph, for n > 2. 


Proof Let V(P,) = {v1,v2,°-:,Un} and E(P,) = fe; = vuivigi;1 < i < n— 1} be the 
vertex set and edge set of the path P,. Then 


V(T(Pn)) _ {U1, V2, +++) Un, €1, €2,°°° ,€n—1} and 
E(T(Pr)) = {uivigi, eri, esViga; 1 <i<n— UU feeni;l1<i<n—2}. 


We define f : V(T(P,)) — {1, 2,3,--- ,4(n — 1)} as follows: 


f(y) = 4¢- 3, for 1 <i<n-1, 
f(vn) = 4n —4 and 
f(e;:) =4¢-1, forl<i<n-1. 
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The induced edge labeling is as follows: 





f* (vivigi) = 41-2, for 1 <i<n-1, 
f* (eieiz1) = 44, for 1 <i<n—-2, 

f* (evi) = 44-3, for 1 <<i<n-—1 and 
f(evign) = 4i-—1, for l<i<n-1. 


Hence, f is a geometric mean labeling of the graph T(P,,). Thus the graph T(P,,) is a geometric 











mean graph, for n > 2. 





A geometric mean labeling of T(P5) is shown in Fig.7. 


V1 2 v2 6 U3 to". ¥4 14 U5 
1 16 


12 


el €2 8 €3 €4 


Fig.7 


Theorem 2.7 Tadpoles T(n,k) is a geometric mean graph. 


Proof Let uy,ug,-++ ,Un and v1, v2,°-- , Up be the vertices of the cycle C;, and the path P, 
respectively. Let T(n, k) be the graph obtained by identifying the vertex uy, of the cycle C,, to 
the end vertex v; of the path P;,. We define f : V(T(n,k)) — {1,2,3,---,n+k} as follows: 


and 


f(ui) = 


i if 1<i<|Vn+1]-1 
i+1 if |¥Wn+1|<i<n 


flu) =n+i, for2<i<k. 
The induced edge labeling is as follows: 
i if 1<i<|/nti|-1 
i+1 if |Vn+1| <i<n-1, 
f*(urun) = |vn+1| and 


f* (wivigi) =n +i, forl<i<k-1. 


f* (uit) = 


Hence, f is a geometric mean labeling of the graph T(n, k). Thus the graph T (n, k) is a geometric 











mean graph. 





A geometric mean labeling of the Tadpoles T(7, 5) is shown in Fig.8. 
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Fig.8 


Theorem 2.8 The graph obtained by identifying a verter of any two cycles Cm and Cp, is a 
geometric mean graph. 


Proof Let ui, U2,-++ ,Um and v1, vV2,°++ ,Un be the vertices of the cycles Cy, and C;, re- 


spectively. Let G be the resultant graph obtained by identifying the vertex u,, of the cycle Cm 
to the vertex vu, of the cycle C,,. We define f : V(G) > {1,2,3,---,m+n+1} as follows: 


i iP Lt < pym-a) 1 











f(ui) = and 
i+1 if |Vm+1| <i<m 
m+1+i if 1<i m —m-—2 
ae 1 S <i<|Vmm+DmsnszD| a 
m+2+i% if | (m+DmFn+0| m—-1<i<n-1. 


The induced edge labeling is as follows: 





. i if 1<i<|vm¥i]-1, 
Putin) = 7 | 
i+1. if | m+1|<i<m-1, 
m+1+72 if sis | (m+n +n+D| -m—2, 
f° (viviga) = , ; 
m+2+% if | (mF Dm+n+D| —-m-1<i<n-2, 





Pasa, 
f° @a10x) = | (m+ 1)\(m+n+ D| and 
f*(v10n) = m+ 1. 


Hence, f is a geometric mean labeling of the graph G. Thus the resultant graph G is a geometric 
mean graph. 














A geometric mean labeling of the graph G obtained by identifying a vertex of the cycles 
Cg and C2, is shown in Fig.9. 


94 A.Durai Baskar, S.Arockiaraj and B.Rajendran 





Fig.9 


Theorem 2.9 The graph obtained by identifying an edge of any two cycles Cy, and C,, is a 
geometric mean graph. 


Proof Let ui, U2,-++ ,Um and v1, vV2,°-:,Un be the vertices of the cycles Cy, and C), re- 
spectively. Let G be the resultant graph obtained by identifying an edge um—1Um of cycle Ci, 
with an edge vp_1Un of the cycle C,,. We define f : V(G) > {1,2,3,---,m+n} as follows: 


if l<i<|Vm+41]-1 
f(ui) = : ; S| ie and 
i+1 if |Vm+1|<i<m 


m+1+i if 1<i<|¥mn¥n)|—m-2 


f(vi) = Pare rir it | m{m+n)|—m-1<i<n—2 








The induced edge labeling is as follows: 


et 
*(uwiti1) 
i+1 if |\Vm+1| <i<m-1, 


m+1+i if 1<i<| m{m+n|—m—2 
*(vivi4i) 








m+2+i if | mm +n)|—m-1<i<n—3, 
f*(uitm) = |vm +1], 
f*(vivn) =m-+1 and 

f* n—2Un-1) = | m(m+n)| : 


Hence, f is a geometric mean labeling of the graph G. Thus the resultant graph G is a geometric 
mean graph. 














A geometric mean labeling of the graph G obtained by identifying an edge of the cycles 
Cio and Cis, is shown in Fig.10. 
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14 





Fig.10 


Theorem 2.10 The graph obtained by joining any two cycles Cy, and C,, by a path Py is a 
geometric mean graph. 


Proof Let G be a graph obtained by joining any two cycles C,, and C,, by a path Px. 


Let uy, U2,°°*,Um and v1, vV2,°-:,Un be the vertices of the cycles C,, and Cy respectively. 


..,Wr be the vertices of the path Py, with u,», = w, and we, = vn. We define 
f:V(G) = {1,2,3,---,m+k-+n} as follows: 


Let wi, w2,. 


Flu) i if 1<i<|vmFl|-1 
Ui = 
i+1 if |¥m+1| <i<m, 


f(w;) =m-+i, for2<i<k and 


Pee cae, ~ 1<i<|Vim+honskFn)| " es 
m+t+k+1+i if | (m+ km+k+n)| —m—-k<isn-1. 














The induced edge labeling is as follows: 


: i if 1<i<|vm41|-1 
Pf (uitigi) = 4. 
i+1 |: m+1|<i<m-1 


% 


f*(wiwig1) =mt+i, for l<i<k—-I, 





m+k+i if ee | (mF Rm+k+n)| —m—k-1 
Ff (vivi4a) = oe 
m+k+1+i if | (m+ Bm +k n)| m—-k<i<n-2, 
f" Qi ten) | m+], 











f* (UnUn—1) = | (m + k)(m +h+ n)| and 
f*(vivn) =mM+k. 


Hence, f is a geometric mean labeling of the graph G. Thus the resultant graph G is a geometric 
mean graph. 
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A geometric mean labeling of the graph G obtained by joining two cycles C7 and Cio by 
a path P,, is shown in Fig.11. 





Fig.11 


Theorem 2.11 Any H-graph G is a geometric mean graph. 


Proof Let uy, u2,-++ ,Un and v1, V2,°+: ,Un be the vertices on the paths of length n in G. 
Case 1 n is odd. 
We define f : V(G) > {1,2,3,--- ,2n} as follows: 
f(u;) =%, for 1 <i<nand 
n+2i if 1<i< |Z] 
fw)=4 n+2%-1 if i=[F]41 
3n+1—2i if [Z| +2<i<n. 
The induced edge labeling is as follows: 
f* (uiigi) =%, for 1 <i<n-1, 
f* (ui) =n, for i= =| +1 and 
n+ 2i if 1<is|t| 


f* (wivigi) = 2. oy n ; 
38n—1-—27 if [S| +isisn-1. 


Case 2. n is even. 
We define f : V(G) > {1,2,3,--- ,2n} as follows: 
f(u;) =%, for 1 <i<nand 
n+2i if 1<i< |Z] 


f(vi) = 
Speco: i [5 | +1si<n 
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The induced edge labeling is as follows: 
f* (uimi4i) = 1, for 1 < a < n— 1, 
f* (uizius) =n, fori = 5 | and 
: ; n 
: nt 2i if 1<i<|S|-1 
f° (vivi4i) = ao eee 
38n—1-—2i if [S| <isn-1 





Hence, H-graph admits a geometric mean labeling. 











A geometric mean labeling of H-graphs G and G2 are shown in Fig.12. 


1 eu Cm me 





1 LT 
2 13 
2 13 
3 15 
3 15 
4 17 
4 17 
5 18 
5 16 
6 16 
6 14 
7 14 
7 12 
8 12 
8 10 
9 ug @ 10 
Gy 
Fig.12 
Theorem 2.12 For any n > 2, P(1,2,3,---,2—1) is a geometric mean graph. 
Proof Let v1,v2,-*:,Un be the vertices of the path P, and let u,; be the vertices of 


the partition of Kom, with cardinality m;,1 <7<mn-—1l1and1< j < mj. We define f : 
V(P(1,2,---,n-—1)) > {1,2,3,...,n(n — 1) + 1} as follows: 
f(u;) =i@-1)4+1, for l<i<nand 
f(uiy) =i -—1) +27, forl<j<tandl<i<n-1. 





The induced edge labeling is as follows: 


f (usu) =i@-1)+ 9, forl<j<iandl<i<n-1 
f*(uijvign) =? +9, forl<j<iand1<i<n-1. 
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Hence, f is a geometric mean labeling of the graph P(1, 2,--- 
1) is a geometric mean graph. 
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A geometric mean labeling of P(1,2,3,4,5) is shown in Fig.13. 





Fig.13 
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Abstract: In this paper, we prove that the Complete expansion graph And(k) (k > 6) is 


4-ordered hamiltonian graph by the method of classification discuss. 
Key Words: Andrdsfai graph, complete expansion graph, k-ordered hamiltonian graph. 
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§1. Introduction 


All graphs considered in this paper are finite, simple and undirected. Let C be a cycle with 
given orientation in graph X, Ce = C) with anticlockwise direction and C with clockwise 
direction. If z € V(C), then we use x* to denote the successor of x on C and x~ to denote its 
predecessor. Use C[x, y] denote (x, y)-path on C; C(x, y) denote (a, y)-path missing x,y on C. 
Any undefined notation follows that of [1, 2]. 


Definition 1.1({1]) Let G be a group and let C be a subset of G that is closed under taking 
inverses and does not contain the identity, then the Cayley graph X(G,C) is the graph with 
vertex set G and edge set E(X(G,C)) ={gh: hg! € C}. 


For a Cayley graph G, it may not be a hamiltonian graph, but a Cayley graph of Abelian 
group is a hamiltonian graph. And(k) is a family of Cayley graph, which is named by the 
Hungarian mathematician Andrdsfai, it is a k-regular graph with the order n = 3k — 1 and it 


is a hamiltonian graph. 


Definition 1.2({1]) For any integer k > 1, let G = 23,1 denote the additive group of integer 
modulo 3k—1 and let C' be the subset of Z3~-1 consisting of the elements congruent to 1 modulo 
3. Then we denote the Cayley graph X(G,C) by And(k). 


For convenience, we note Z3,-1 = {uo, U1,..-,U3k—2}. For uj,uy € V[And(k)], ui ~ uy if 





and only if 7 —i = +1mod3. The result are directly by the definition of Andrdsfai graph. 


1Supported by Natural Science Foundation of Inner Mongolia, 2010MS0113; Inner Mongolia Normal Univer- 
sity Graduate Students’ Research and Innovation Fund. CXJJS11042. 
2Received November 7, 2012. Accepted March 15, 2013. 
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Lemma 1.3 Let C be any hamiltonian cycle in And(k)(k > 2). 


(1) IfVu,a € V(And(k)), u~ x is a chord of C, thenu~ ~ a7, ut~art. 
(2) IfVu,x,y € V(And(k)), u~ x, u~y are two chords of C, then x ~ yt. 


The definition of k-ordered hamiltonian graph was given in 1997 by Lenhard as follows. 


Definition 1.4([3]) A hamiltonian graph G of order v is k-ordered, 2 < k < v, if for ev- 
ery sequence (U1, V2,...,Un) of k distinct vertices of G, there exists a hamiltonian cycle that 


encounters (U1,U2,...,Uk) in this order. 


Faudree developed above definition into a k-ordered graph. 


Definition 1.5((4]) For a positive integer k, a graph G is k-ordered if for every ordered set of 


k; vertices, there is a cycle that encounters the vertices of the set in the given order. 


It has been shown that And(k)(k > 4) is 4-ordered hamiltonian graph by in [5]. The 
concept of expansion transformation graph of a graph was given in 2009 by A Yongga at first. 
Then an equivalence definition of complete expansion graph was given by her, that is, the 
method defined by Cartesian product in [6] as follows. 


Definition 1.6([6]) Let G be any graph and L(G) be the line graph of G. Non-trivial component 
of GOL(G) is said complete expansion graph (CEG for short) of G, denoted by 0(G), said the 


map V be a complete expansion transformation of G. 














The proof of main result in this paper is mainly according to the following conclusions. 


Theorem 1.7({1]) The Cayley graph X(G,C) is vertex transitive. 
Theorem 1.8([5]) And(k)(k > 4) is 4-ordered hamiltonian graph. 
Theorem 1.9([7]) Every even regular graph has a 2-factor. 


The notations following is useful throughout the paper. For u € V(G), the clique with 
the order dg(u) in ¥(G) by u is denoted as V(u). All cliques are the cliques in 0(And(k)) 
determined by the vertices in And(k), that is maximum Clique. For u, v € V(G), 0(u) ~ Vv) 
means there exist z € V(0(u)), y €E V(0(v)), s.t.2 ~ y in V(W(G)), edge (x, y) is said an edge 
stretching out from J(u). Use Gy u)[2, y; 8, t] to denote (x, y)-longest path missing s,¢ in J(u), 
where x,y, s,t € V(vV(u)). 


§2. Main Results with Proofs 


We consider that whether 0(And(k)) (k > 4) is 4-ordered hamiltonian graph or not in this 


section. 


Theorem 2.1 J(And(k)) (k > 6) is a 4-ordered hamiltonian graph. 
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The following lemmas are necessary for the proof of Theorem 2.1. 


Lemma 2.2 For any u € V(And(k))(k > 2), Va,y € N(u), there exists a hamiltonian 
cycle C in And(k), s.t. ua € E(C) and uy € E(C). 


Proof Let Co is a hamiltonian cycle ug ~ uy ~ ug ~... ~ Usr—2 ~ Uo in And(k)(k > 2). 
For u € V(And(k))(k > 2), Va, y € N(u), then we consider the following cases. 


Case l r~u~yon Oe: Then C = Cp is that so, since Co is a hamiltonian cycle. 


Case 2 c~uandud yon Co ora £uandu~ yon Cp. Ifa~uandu¢ yon Co, then 
we can find a hamiltonian cycle C in And(k)(k > 2) according to Lemma 1, that is, 


C=uraerOo(ay)~y ~u ~Colu y)~yr~u 


Ife #uandu~ yon Ch; then we can find a hamiltonian cycle C in And(k)(k > 2) according 
to Lemma 1.3, that is, 


C=urarn Co(x,u*) wut natn Co(at,y)~yru. 
Case 3 cure yon Ch. Then we can find a hamiltonian cycle C in And(k)(k > 2) according 
to Lemma 1.3, that is, 


Caunanyt~Colyt,u nu ~a ~ Co(a-,ut) wut ~ Colet,y] ~ u. 














For any u € V(And(k))(k > 2), Lemma 2.2 is true since And(k) is vertex transitive. 





Corollary 2.3 For any two edges which stretch out from any Clique, there exists a hamiltonian 
cycle in 0(And(k)) containing them. 


zcuns 2.4 Ifk is an odd number, then And(k)(k > 3) can be decomposed into one 1-factor 

-—1 
2 

Proof 3k —1 is an even number, since & is an odd number. There exists one 1-factor 

M in And(k) by the definition of And(k). According to Theorem 1.9 and the condition of 

Lemma 2.4 for integers k > 3, And(k) — E(M) is a (k — 1)-regular graph with a hamiltonian 

cycle Ci, And(k) — E(M) — E(C;) is a (k — 3)-regular graph with a hamiltonian cycle C2, ---, 





and 2-factors. 


k-1 
And(k) — E(M) — >> E(C;) is an empty graph. 
i=1 


Assume k = pa 1(r € Z*), since k is an odd number. First we shall prove the result for 
r = 1, and then by induction on r. If r = 1 (k = 3), it is easy to see that And(k) — E(M) is 
a hamiltonian cycle C; by Theorem 1.9 and the analysis form of Lemma 2.4 , so the result is 
clearly true. 

Now, we assume that the result is true if r = n(r > 1,k = 2n+1), that is, And(2n+1) can 
be decomposed into one 1-factor and n 2-factors. Considering the case of r = n+1(k = 2n+3, 
we know And(2n+3)(And[2(n+1)+1]) can be decomposed into one 1-factor and n+1 2-factors 
according to the induction. 
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Thus, if & is an odd number, then And(k)(k > 3) can be decomposed into one 1-factor and 


-1 
‘— 2-factors. 














Proof of Theorem 2.1 


0(And(k)) is a hamiltonian graph, since And(k) is a hamiltonian graph. So there exists 
a hamiltonian cycle Co in 0(And(k)) and a hamiltonian cycle Co’ in And(k), such that Cp = 
0(Co’), without loss of generality 


/ 
Co = Uo... Ugk—2U0, 


then 

Co = U0,1U0,2---U0,hU1,1U1,2-.-U1Lk~.-U3k—2,1U3k—2,2 --- U3k—2,kU0,1; 
where ui; € V(U(ui)), ue € V(And(k)), danacky(u) = k => 6,7 = 0,1,2,...,3k — 2, and 
Un Ur ¢ V(O(us)), us = Uijti (1 < 7 < k-1) and Us) = Uil-1 (2<l1<k). There 
are three cyclic orders Vua,p, c,d; Ue, f,Ug,h © V[U(And(k))] according to the definition of the 
ring arrangement of the second kind, as follows: (Ua,b, Uc,d; Ue, fs Ug,h), (Ua,b, Ue, f; We,d) Ug,h)s 
(Ua,b; Uc,d; Ug,hy Ue, ¢)(see Fig.1). Let S = {(ua,b, Ue,d; Ue, fy Ug,h)s (Ua,bs Ue, f, Ue,d; Ug,h), (Ua,bs Ue,ds 


Ug,hs Ue,f) } 


Ua,b Ua,b Ua,b 
Ue,d Ug,h Ue, f g,h Uc,d Ue, f 
Ue, f Uc,d Ugh 
(Ua,b; Uc,d; Ue, fs Ug,h) (ta,b» Ue, f, Uc,ds Ug,h) (ta,b; Uc,d, Ugh» Ue, f) 


Fig.1 Three cyclic orders 


Now, we show that 4-ordered hamiltonicity of 0(And(k)) (k > 6). In fact, we need to 
prove that a € S, there exists a hamiltonian cycle containing a. Without loss of generality, 
hamiltonian cycle Co encounters (ta,b, Ue,d, Ue, f;Ug,h) in this order. So we just prove: V3 € 
S'\ (ta,b, Uc,d; Ue, f,Ug,n), there exists a hamiltonian cycle containing (. 

According to the Pigeonhole principle, we consider following cases. 


Case 1 If these four vertices ua,p, Uc,d, Ue,f,Ug,h are contained in distinct four Cliques of 
0(And(k)), respectively. And Theorem 2.1 is true by the result in [5]. 


Case 2 If these four vertices Ua,p, Uc,d, Ue, f, Ug,h are contained in a same Clique of ¥(And(k)), 
thena=c=e=g9,b<d<f<h. Let S = {(uap, tad, Ua,f, Ua,h); (Wa,b; Ua, fs Wad, Wa,h)s (Ua,b; 
Ua,d, Ua,h; Ua,f)}- 


(1) For (ua,p, Ua,d; Ua, f; Ua,h) € S. Co is the hamiltonian cycle that encounters (Ua,p, Ua,d; Ua, fs 
Ua,h) in this order, clearly. 
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(2) For (ta,b, Ua, f; Ua,d; Ua,h) € S. We can find a hamiltonian cycle 


C= Ua,bC'o (ta,b, Ug a)Ua,aua,fCo(Ua,f, Ua,d)Ua,dug ¢Co (uz ps Ua,h)Ua,nCo(Ua,h, Ua,b)Ua,b 


that encounters (ta,b, Ua,f, Ua,d; Ua,n) in this order. 


(3) For (ta,b, Ua,d; Ua,h; Ua, f) € S. We can find a hamiltonian cycle 


C = Ua,1Co(Uar, Ua, f)Ua, pUakCo(Ua,k, Ua, f Ua, fP(C1) (Ua, f, Ua,1)Ua,1 


that encounters (Ua,p, Ua,d; Ua,h; Ua,f) in this order by Lemma 2.2 and Corollary 2.3 (see Fig.2). 


v(C1) 





Fig.2 C= ug iCo(ua1; Ug, p) Ua, fUa,kCo(Ua,k, Ua, f Ua, fO(C1) (Ua, f, Ua,1) Ua, 


Case 3 If these four vertices ua,p, Uc,d, Ue,f,Ug,h are contained in distinct two Cliques of 
0(And(k)), without loss of generality, we assume that wa, Uc € V(U(ua)) and wef, tg,n € 
V(0(ue)) in V(And(k)) or Ua,b, Ue,d; Ue, f € V(U(Ua)) and ug, € V(V(u,)) in V(And(k)) accord- 


ing to the notations. Let S = {(ua,b, Ua,d, Ue, f, Ue,h); (Wa,b; Ue, fy Ua,ds Ue,h); (Wa,b; Ua,d Ue,h; Ue, f)}- 
Subcase 3.1 wa,b,Uc,a € V(U(ua)) and te, f, Ug,n € V(U(ue)) in V(And(k)). 


(1) For (ua,p, Ua,d; Ue, f; Ue,h) € S. Co is the hamiltonian cycle that encounters (u,b, Ue,d; 
Ue, f,Ug,n) in this order, clearly. 


(2) For (ta,b, Ue, f, Ua,d, Ue,h) € S. Let Cy is a hamiltonian cycle in And(k) or And(k) — 
E(M), C2 is a hamiltonian cycle in And(k) — E(C1) or And(k) — E(M) — E(C1)(see Fig.3). 
Use A(C1) to denote a cycle that only through two vertices in 0(u;)(i = 1,2,...,3k — 2) and 
related with 0(C1), and use A(C2) to denote the longest cycle missing the vertex on A(C‘) 
in V(And(k)) or Y(And(k)) — M(see Fig.3). We suppose that P, = [x,y], P2 = [p,ua.»] on 
cycle A(C1) in ¥(And(k)) or ¥(And(k)) — M and P3 = [m,n], Ps = [s,t] on cycle A(C2) in 
0(And(k)) — A(C1) or 0(And(k)) — M — A(C1) by Theorem 3!1, the analysis of Lemma 2.4 and 
the definition of CEG (see appendix). Now, we have a discussion about the position of vertex 
x, Y, p, s and n in V(And(k)). 
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C— 


0 
uy USk—2 


whe oo 





Fig.3 


In where, C; in And(k) or And(k) — E(M), C2 in And(k) — E(C;) or And(k) — E(M) — E(C1), 
A(C1) in ¥(And(k)) or 0(And(k))—M, A(C2) in 0(And(k))— A(C)) or 0(And(k))— A(C1) —M. 


Pte peieeeatoton ele meie ees: (5) 
DiS tae, Sania asesrevinsenidens (6 

re pF Ue,f; SS Ue Pyiaiinadaead whi wed biants (7 

1S Nef a tiaa wine set das (8) 

Bi 2 ops ahietneatuedasiwes (9) 

PH Ue fr) S= We hjie sir dias canoes ceeds (10) 

ZF Ua,b; Uad; Lee re Se eee ere ee (11) 
Bie ep ho nye rudadete (12) 

Deny) 6 Ue pyr cerciecrir erage dans (13) 

Ne Ip eh seertaieie a sateen cia (14) 

Y F Ue, f, Ue,hs S= Ue fpr rere (15) 
WSUS Poti ondies date (16) 

Cerne tee ea (17) 

pF Ue, f, Ue,hs WS We fi isiadetecsdaaees (18) 

8 = We FU = Wehr eivens (19) 

8 = UWerny = Ue fy eis ees (20) 

SF Ue,f, 1% F Wenge eveees (21) 
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For cases (1) and (2), we can find a hamiltonian cycle 
Ua bt Pi (a, y)ysPa(s, t)tGyiu,)[t, 5 Uap, Z]MP3(m, N)NGgu,)[N, D3 Y, 8|pP2(p, Ua,b)Ua,b 


that encounters (Ua,b, Ue, f, Ua,d; Ue,n) in this order. 

For cases (3)-(21), we can find a hamiltonian cycle that encounters (ti¢,5, Ue, f, Wad, Ue,h) in 
this order according to the method of (1) and (2). 

(3)For cases (2)-(11) and (15)-(21), we can find a hamiltonian cycle that encounters 
(Ua,b; Ua,d; Ue,h; Ue, f) in this order according to the method of Case3.1(2). 

For case (1), we can find a hamiltonian cycle 


Ua,bG9(u2) [Wa,b; ™, t|mP3 (m, n)nG'g(ue) [n, PY, s]pysP,(s, t)tuap 


that encounters (ta,b, Ua,d, Ue,h, Ue,f) in this order. P! is the path which through the all vertices 
in V(ui)(i =a,...,e) and related with P;(i = 3,4) in 0(And(k)) — A(C) or 0(And(k)) — M — 
A(C;)(see Fig.4). 


hy Jt) = (ua) 


aiugytee ¥ 
at or > 
2 AY 
rN Hu} ) 


d(w;) 


* (ue) 
. o 





Fig.4 


In where, P;, Pz in 0(And(k)) or 0(And(k)) — M, Ps, Py in 0(And(k)) — A(C)) or 0(And(k)) — 
M— A(C;), P3, Py related with P3, Py in 0(And(k)) — A(C1) or 0(And(k)) — M — A(C;). 

For 12-14, we can find a hamiltonian cycle that encounters (Ua,p, Ua,d; Ue,h; Ue, f) in this 
order according to the method of 1. 


Subcase 3.2 tap, Uc,d, Ue, € V(U(ua)) and ug.n € V(V(u,)) in V(And(k)). For all condition , 
we see the result is proved by the method of Subcase 3.1. 


Case 4 If these four vertices Ua,b, Uc,d; Ue, f,Ug,h are Contained in distinct three Cliques of 
0(And(k)). Without loss of generality, we assume that ua,p, Uc,a € V(U(ua)), Ue,f € V(V(ue)) 
and ug.rn € V(Y(ug)) in V(And(k)). 


(1) For (ua,, Ua,d; Ue, f,Ug,h) € S, Co is the hamiltonian cycle that encounters (Ua,p, Ua,d; 
Ue, f,Ug,n) in this order, clearly. 


(2) For (ta,b, Ue, f, Ua,d, Ug,n) € S. Let Cy is a hamiltonian cycle in And(k) or And(k) — 
E(M), C2 is a hamiltonian cycle in And(k) — E(C,) or And(k) — E(M) — E(C1), C3 is a 
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hamiltonian cycle in And(k) — E(C,) — E(C2) or And(k) — E(M) — E(C1) — E(C2)(see Fig.5). 
Use A(C;) to denote a cycle that only through two vertices in ¥(u;)(¢ = 1,2,...,3k — 2) and 
related with ¥(C;)(j = 1,2), and use A(C3) to denote the longest cycle missing the vertex 
on A(C;,) and A(C2) in ¥(And(k)) or B(And(k)) — M(see Figure5). We can suppose that 
P, = [uca, 2], Po = [y, Uap] on cycle A(C1) in ¥(And(k)) or V(And(k)) — M , P3 = [m,n], 
Py = [p,q] on cycle A(C2) in 0(And(k a A(C) or 0(And(k)) — M — A(C)) and Ps = {s, ¢], 
Ps = [w, z] on A(C3) in 0(And(k)) — » A(C;) or 0(And(k)) — M — > A(C;) by Theorem 3!7), 


= 


the analysis of Lemma 3 and the defeition of CEG (see appendix). Now, we have a discussion 





about the position of vertex m, q, x, y, p and n in 0(And(k)). 


L = Ugh Y A Ugh, 51 Mes By oles Rec Mee D eae SO aes Baas (1) 
Y= Ug hy rrr eee (2) 
m, q x Ua,b, Ua,d) P= Ugny erie te (3) 
zr # Ug,h>s 
y # Ug,h3 NM=Ughy crt (4) 
DNF Ugne corre (5) 
Cc; ic &(ue) 


Ma —_ 


— ty 


a: 
“i a iny 9 ei —It—~> 
—- eas . AO a +02 ee U(us) ce >) 


ok 
PsP Muy 
nf agok 


Wi v4 fy 


— 


en 
on aul) 


> 








Ole, } 
Fig.5 


In where, C; in And(k) or And(k) — E(M), C2 in And(k) — E(C1) or And(k) — E(M) — E(C1), 
C3 in And(k) — E(C,) — E(C2) or And(k) — E(M) — E(C,) — E(C2), A(C1) in 0(And(k)) 
0(And(k))— M, A(C2) in 0(And(k))— A(C1) or 0(And(k)) — A(C1) — M, A(C3) in 0(And(k) 
A(C1) — A(C2) or 0(And(k)) — A(C1) — A(C2) — M. 

For case (1), if we, € V(Pi) (¢ = 2,3,4), we can find a hamiltonian cycle that encounters 


or 


= 





(Ua,b; Ue, f; Ua,d; Ugh) in this order according to the method of Subcase 3.1,(2). 


If ue,p € V(P1), we can find a hamiltonian cycle 
Ua,bdP4(q, p)pnP3(n,m)MG'y(u,)[M, 8; Ua,b, WS-P3(s, ttGy(u,)[t, ys p, n, tlyPo(y, Ua,b)Ua,b Or 


Ua,pmP3(m, n)npP, (q, P)QG'9(ua) [q, 5S; Ua,b> m|sPz (s, t)tG9(u,) [t, Y3P, nN, tlyP; (y, Ua,b) a,b 
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that encounters (Ua,p, Ue, f; Ua,d; Ugh) in this order. There exist some vertices which belong to 
a same Clique on P;,P; and P;(t = 3,4;7 = 5,6). And ue € V(P/)(i = 3 or 4). Pl is 


a 


the path which through the all vertices in J(u;)(¢ = a,...,g) and related with P;(i = 5,6) in 
2 2 
0(And(k)) — > A(C;) or 0(And(k)) — M — > A(C;), and missing the vertex on P3, Pj(refers 
i=1 i=1 
to Figure4). 


For cases (2)-(5), we can find a hamiltonian cycle that encounters (Ua,b, Ue, f, Wa,d, Ug,h) in 


this order according to the method of (1). 


(3) For cases (1)-(5), we can find a hamiltonian cycle that encounters (t¢,b, Ua,d, Ug,hs Ue, f) 











in this order according to the method of Case 4(2). 
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Appendix 


By the theorem 1.9, the analysis of Lemma 2.4, the definition of CEG, And(k) and the parity 
of k(s € Z*) , we know that 


And(4) — E(C1) = C2 0(And(4)) — 








s=2 
And(6) — E(C1) — E(C2) = Cs | 0(And(6)) 





5] Wang lei,A Yongga, 4-Ordered Hamiltonicity of some Cayley graph , Int. J. Math. Comb.,1(2007), 


108 Lian Ying, A Yongga, Fang Xiang and Sarula 





s=3 


If k is odd, it should be illustrated that the M’s selection method, that is, M satisfy 
condition Ua,b, Ue,d, Ue, f,Ug,n £ V(M) in V(And(k)). It can be done, because k > 7. 


And(k) B(And(k)) 

And(3) — E(M) = C1 9(And(3)) — M = B(C21) 

And(5) — E(M) — E(C1) = C2 0(And(5)) — E(M) — B(C1) = B(C2) 

And(7) — E(M) — E(C1) — E(C2) = C3 | 8(And(7)) — M — B(C1) — B(C2) = B(Cs) 
i BC) =Cn 














) 
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Abstract: In this article, we present algorithms for equitable weak product graph of cycles 
Cm and Cn, Cm Xx Cn such that it has an equitable chromatic value, x=(Cm x Cn) = 3, with 


mn odd and m or n is not a multiple of 3. 
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§1. Introduction 


Let G be a graph with vertex set V(G) and edge set E(G). A k-coloring on G is a function 
f :V(G) = [1,4] = {1,2,---,&}, such that if ww € E(G),u,v € V(G) then f(u) 4 f(v). A 
value x(G) = k, the chromatic number of G is the smallest positive integer for which G is 
k-colorable. G is said to be equitably k-colorable if for a proper k-coloring of G with vertex 
color class V|,V2---Vz, then |(|Vi| — |V;|)| < 1 for all 7,7 € [2,k]. Suppose n is the smallest 
integer such that G is equitably k-colorable, then n is the equitable chromatic number, y=(G), 
of G. 

The notion of equitable coloring of a graph was introduced in [6] by Meyer. Notable work on 
the subject includes [7] where outer planar graphs were considered and [8] where general planar 
graphs were investigated. In [1] equitable coloring of the product of trees was considered. Chen 
et al. in [2] showed that for m,n > 3, y=(Cm x Cy = 2) if mn is even and y=(Cm x Cr = 3) 
if mn is odd. Recent work include [4], [5]. Furmanczyk in [3] discussed the equitable coloring 
of product graphs in general, following [2], where the authors separated the proofs of mn into 


various parts including the following: 
1. m,n odd with n = 0 mod 3 
2. m,n odd, with 


(a) either m or n, say n satisfying n —1=0mod 3 
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(b) either m or n, say n satisfying n — 2 =0 mod 3. 
In this paper we present equitable coloring schemes which 
1. improve the proof in (b) above and 


2. can be employed in developing the equitable 3-coloring for Cy, x C, with mn odd. 


§2. Preliminaries 


Let G; and G2 be two graphs with V(G) and E(G,) as the vertex and edge sets for G; respec- 
tively and V(G2) and E(G2) as the vertex and edge sets of G2 respectively. The weak product 
of G; and G2 is the graph G; x G2 such that V(G1 x G2) = {(u,v) =u € V(G)and u € V(G2)} 
and E(G, x G2) = 
{(u1v1)(ugv2) : uug € E(Gy)and viv2 € E(G2)}. A graph Py», = upuru2-+:Um-_i is a path of 
length m— 1 if for all ui,v; € V(Pm),t Aj. A graph Cy, = uouiu2--+Um-1 is a cycle of length 
m if for all w,v; € V(Cm),i #7 and upttm—1 € E(Cm). 

The following results due to Chen et al gives the equitable chromatic numbers of product 
of cycles. 


Theorem 2.1([2]) Let m,n > 3. Then 


2 if mn is even 


xX=(Cm Xx Cn) = 
3 if mn is. odd. 


We require the following lemma in the main result. 


Lemma 2.2 Let n be any odd integer and letn—1=0mod 3. Thenn—1=0mod 6. 


Proof Since n is odd, then there exists a positive integer m, such that n = 2m+1. Now 
since n is odd then, n — 1 is even. Let 2m =0 mod 3. Clearly, n > 3. Now 2m = 3k where k 
is an even positive integer. Thus 2m = 3(2k’) for some positive integer k’ and thus 2m = 6k’. 
Hence n — 1 = 6k’. 














§3. Main Results 


In this section, we present the algorithms for the equitable 3-coloring of Cy, x Cy with where 
m and n are odd with say n -1=0mod 3 andn—2=0mod 3. 


Algorithm 1 Let Ci, x C;, be product graph and let mn be odd, with n—1=0mod 3. 
Step 1 Define the following coloring for ujvj € V(Cm x Cn). 
aq for {ujv;: 7 € [n—1];7 > 5;7 +1=0 mod 3} 
f(uivj) = 4 ay for {uyvj: 7 € [n-1];7 +2 =O mod 3}U {ujve : i € [m— 1} 
a3 for {ujv;: 7 € [n—1];7 > 6;7 =O mod 3}U {uu1,7 € [m—1]}. 
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Step 2 For all u;vo;7 € [2], define the following coloring: 


(a) 
f(ujvo) = a, for =f 
az fori =0,2 
(b) 
f (ujv3) = a. fori=1,2 


a3 fori =O 


Step 3 Repeat Step 2(a) and Step 2(b) for all u;vp and u,v3 for each i € [x, 2 + 2] where 
x =Omod 3. 


Proof of Algorithm 1 Suppose n is odd and n—1=0 mod 3. From Lemma 2.2 above, 


n—1=0 mod 6 and consequently, n — 4 = 0 mod 3. Suppose n—4 =n’, where n is a 


positive integer. Let P,, x P,—4 be a subgraph of P,, x P,, where P,_4 = v4v5+++Un—1.. For 
all uv; € V(Pm X Ppa), let 





a, for {uv;:7 € [n—-1],7+2=0mod 3} 
f(uivj) = 4 ag for {usvy: 7 € [n—-1],7 > 5; +1=0 mod 3} 
a3 for {ujv;:7 € [n—1];7 > 6;7 =O mod 3}. 


From f(u;v;) defined above, we see that P,, x Pra is equitably 3-colorable with color 





set {a1,Q2,a3} = [1,3], where |Va,| = |Va.| = |Va,| = mn’. Next we show that there exists 





a 3-coloring of P,, x P, that merges with P,, x Py whose 3-coloring is defined by f(ujv;) 
above. First, let F(P3 x P,) be the 3- coloring such that 


a2 a3 A, ag 
F(P3 x Ps) = a1 a3 ay ay 
a2 a3 At a3 


From F'(P3 x Ps) we observe for all 7 € [3], that for F(uov;) C F(P3 x P),|Va,| = 
1,|Voo| = 1,|Vag| = 2; for F(urv;) C F(P3 x P1),|Vo,| = 2,|Voo| = 1,|Vag]| = 1; and for 
F(ugv;) C F(P3 x P4),|Voy| = 1, |Vas| = 2, [Vos] = 1. 

We observe, over all, that for F'(P3 x P1),|Vo,| = |Va.| = |Voz| = 4. These confirm that 
Ps x Py, is equitably 3-colorable at every stage of 7 € [2] and that F'(P2 x Py) C F(P3 x Py) is 
an equitable 3-coloring of P2 x Py for both Pp x Py C P3 x P,. Now the equitable 3-coloring of 








Py, x Py is now obtainable by repeating F'(P3 x P,) at each interval [z,7+2], where « = 0 mod 3, 
until we reach m. Clearly, F(ujv3) F(uiva) = @ since ay ¢ F(ujv3). Thus Pi, x Py is equitably 
3-colorable based on the colorings defined earlier. Likewise, F(uiv9) N F(uitn—1) = since 
a3 ¢ F(ujvo). Thus Py x Cy is equitably 3-colorable based on the coloring defined above for 
Pm X Pn. 

Finally, for any m > 3, the equitable 3-coloring of P,, x P,—4 with respect to F'(Pm x Pn—) 
above is equivalent to the equitable 3-coloring of Cm x Cn—a since ujvjuivj+1 ¢ E(Pm xX Pm-—a) 
for all 7 € [n — 5]. Also, for m > 3 the equitable 3-coloring of P,, x P, with respect to 
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F (Pm X Ps) above is equivalent to the equitably 3- coloring of C,, x C4 by mere observation. 
Thus, Cm X Cy is equitably 3- colorable or all positive integer m and odd positive integer n 
such that n —1=0Omod 3. 


Algorithm 2 Let m or n, say n be odd such that n — 2 =0mod 3. 
Step 1 Define the following coloring: 
a, for {uv;:7 € [n—-1],7+1=0mod 3} 
f(uivj) = 4 ag for {uyv;: 7 € [n—1],7 =0 mod 3} 
a3 for {uv;:7 € [n—-1],7-1=0mod 3} 
Step 2(a) For all € [2], let f(uisvo) = a1, a2, a1 respectively a1, a2 € [2]. 
Step 2(b) For all 7 € [2], let f(ujv1) = a3, a2, a3 respectively, a3 € [2]. 


Step 3 Repeat step 2(a) and Step 2(b) above for all 7 € [x,a + 2], where x is a positive 
integer and « = 0 mod 3. 


Proof of Algorithm 2 Let n be odd and let n—2 = 0 mod 3. By f(ujv;) instep 1, Pr x Ph—2, 
where P,-2 = v2U3°++Un—1, is equitably 3-colorable with |Va,| = |Va.| = |Va3| = mn” where 
nl = 25% and F(ujv2) N F(usvn—1) = 0) for all i € [m— 1]. Now, let 








a1 a3 
F(P3 x P2) = a2 ag 
Qa, a3 


It is clear that F'(P3 x P2) above follows from the coloring defined in step 2 of the algorithm 
and that F'(P3 x P2) is an equitable 3-coloring of P3 x Pz where |Va,| = |Va.| = |Vag| = 2. It is 
also clear that F'(P3 x P2) has an equitable coloring at P; x P2 with |Vo,| = 1,|Va.| = 0,|Vo3| = 1 
and at P, x P2 with |Vo,| = 1,|Vo.| = 2,|Vo,| = 1. Now, let with « = 0 mod 3. For all 
x € [m— 1], let f(usv;) = a1,a3 for both 7 = 0,1 respectively; for +1 € [m— 1], let 
f (tx410;) = a2, for 7 = 0,1 and for 7+ 2 € [m— 1], let f(us+2v;) = a1, a3 for 7 = 0,1. With 
this last scheme, we have P, x Pz that has an equitable 3- coloring for any value of m. 








Finally, we can see that P,, x P2, for any m, so equitably, 3-colored merges with Py, x Py—2 
that is equitably 3-colored earlier by f(u;v;), such that F(ujv1) OF (uyve) = 0 for alli € [m—-1 
(by a similar argument as in the proof of Algorithm 1) and F(ujvp) A F(uivn—1) = @ for all 











i € [m— 1] (by a similar argument as in the proof of Algorithm 1). 





Likewise C,, x C;, is equitable 3-colorable (by a similar argument as in the proof of Al- 
gorithm 1). Therefore, C,, <x C, is equitably 3-colorable for any m > 3 and odd n, such that 
n—2=0Omod 3. 


§4. Examples 


In Fig.1, we demonstrate how our algorithms equitably color graphs C5 x C5 and C5 x Cy, 
which are two cases that illustrate n —2=0 mod 3 andn—1=0mod 8 respectively. In the 
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first case, we see that y=(Cs x Cs) = 3, with |Vi| = 8 |V2| = 9 and |V3| = 8 and in the second 
case, x=(Cs5 x C7) = 3, with |Vj| = 12 |V2| = 11 and |V3| = 12. (Note that the first coloring 
takes care of the third instance in subcase 2.4 of [2] where it is a special case.) 





Fig.1 Equitable coloring of graphs C5 x C5 and C5 x C7 
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In this short communication we rectify certain errors which are in the paper, On Set-Semigraceful 
Graphs, International J. Math. Combin., Vol.2(2012), 59-70. The following are the correct ver- 


sions of the respective results. 


Remark 3.2 (5) The Double Stars ST (m,n) where |V| is not a power of 2, are set-semigraceful 
by Theorem 2.13. 


Remark 3.5 (3) The Double Stars ST(m,n) where m is odd and m+n+ 2 = 2!, are not 
set-semigraceful by Theorem 2.12. 

Delete the following sentence below Remark 3.9: ”In fact the result given by Theorem 3.3 
holds for any set-semigraceful graph as we see in the following”. 


Theorem 4.8([3]) Every graph can be embedded as an induced subgraph of a connected set- 
graceful graph. 


Since every set-graceful graph is set-semigraceful, from the above theorem it follows that 


Theorem 4.8A Every graph can be embedded as an induced subgraph of a connected set- 


semigraceful graph. 


However, below we prove: 


Theorem 4.8B Every graph can be embedded as an induced subgraph of a connected set- 


semigraceful graph which is not set-graceful. 


Proof Any graph H with o(H) <5 and s(H) < 2 and the graphs Py, PyU Ki, P3U K2 and 
Ps are induced subgraphs of the set-semigraceful cycle C9 which is not set-graceful. Again any 
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graph H’ with 3 < o(H’) < 5 and 3 < s(H’) < 9 can be obtained as an induced subgraph of 
Hy, V Ky for some graph Hy with o(H,) = 5 and 3 < s(H,) < 9. Then 3 < log, (|E(Mi V K4)|+ 
1) <4, since 8 < s(Hi V K,) < 15 and hence H, V K;, is not set-graceful. By Theorem 2.4, 


4 [log,(|E(A1 V K1)| + 1)| < yi V 1) 
(Ke) (by Theorem 2.5) 


4 (by Theorem 2.19) 


x 


I 


So that H, V Kk, is set-semigraceful. Further, note that Ks is set-semigraceful but not set- 
graceful. 

Now let G = (V,E); V = {v1,...,Un} be a graph of order n > 6. Consider a set- 
indexer g of G with indexing set X = {x,...,%n} defined by g(u;) = {ai}; 1 <a <n. Let 
S = {g(e):e € E}U{g(v): ue V}. Note that |S| = |E£|+n. Now take a new vertex u and join 
with all the vertices of G. Let m be any integer such that 2"~! < m < 2”—(|E|+n+1). Since 


-1 
|E| < n= )) and n > 6, such an integer always exists. Take m new vertices u1,...,Um and 
join all of them with wu. A set-indexer f of the resulting graph G’ can be defined as follows: 


f(u) =9, f(vi)= gui); 1<i<n. 


Besides, f assigns the vertices u1,...,Um with any m distinct elements of 2* \ (SU). 
Thus, 7(G’) < n. But we have 2” > |E|}+n+m+1>m > 27! so that 7(G’) > n, by 
Theorem 2.4. Hence, 


logs (|E(G")| + 1) < [logs (|E(@‘)| + I] =n = 7(G"). 











This shows that G’ is set-semigraceful, but not set-graceful. 





Corollary 4.16 The double fan P, V Kz where k = 2" —m and 2” > 3m; n > 3 is set- 


semigraceful. 


Proof Let G = Py V Ko; Ko = (u1,u2). By Theorem 2.4, 7(G) > flog,(|E|+1)] = 


[log,(3(2” — m) +1)] =n +2. But, 3m < 2” > m < 2.-1—1. Therefore, 


27 — (2+ —2)) <2" -—-m<2"-1 
SSO ee Oy A OD 

= 2)-)41<k-1<2"-2; k=2"%-—m 
=> 27-141 < |E(P,)| < 2" 

= [log,(|E(P.)| + 1)] =n 

=> 7(Pr) =n 











since P;,, is set-semigraceful by Remark 3.2(3). 





Let f be a set-indexer of P;, with indexing set X = {21,...,2,}. Define a set-indexer g of 
G with indexing set Y = X U {2n41, 2n+2} as follows: 


g(v) = f(w) for every v € V(Px), gui) = {angi} and g(u2) = {an+2}. 
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Corollary 4.17 The graph Ky,9n_1 V Kg ts set-semigraceful. 





Proof The proof follows from Theorems 4.15 and 2.33. 











Theorem 4.18 Let Cy where k = 2" —m and 2” +1 > 3m; n > 2 be set-semigraceful. Then 
the graph Cy V Ka is set-semigraceful. 


Proof Let G = Cy V Ko; Ko = (u1,u2). By theorem 2.4, 7(G) > [log,(|E|+1)] = 


[log,(3(2” — m) + 2)] =n +2. But, 3m < 2”"+1=>m < 2"-!. Therefore, 


POR) a1) ye 2st oF 
S22 416k <2"; k= 2" =m 
=> 2-141 < |E(Cy)| < 2” 

= [logs (|E(Ck)| + 1)] =n 

=> 7(Ce) =n 














since Cy is set-semigraceful. 





Let f be a set-indexer of Cy, with indexing set X = {x1,...,2n}. Define a set-indexer g of 
G with indexing set Y = X U {2n41, 2n+2} as follows: 


g(v) = fw) for every v € V(Cx), g(ur) = {ansi} and g(u2) = {tn+42}- 


Corollary 4.21 W,, where 2™—-1<n<2™+42™-!_2:m> 3 is set-semigraceful. 





Proof The proof follows from Theorem 3.15 and Corollary 4.20. 











n-1 





2 
Theorem 4.22 If Wo, where 


of order 2k +1 is set-semigraceful. 


<k <2": n> 4 is set-semigraceful, then the gear graph 


Proof Let G be the gear graph of order 2k + 1. Then by theorem 2.4, 


[logs (3k + 1)] < (G) < y(Wax) (by Theorem 2.5) 








= [log,(4k+1)] (since Wo, is set — semigraceful) 
= [log,(3k + 1)| 
since 
gn-1 
5 <k<2 2 32% 1< 3k < 4k < 2” 
= 21 41< 38kh+1<4k+1< 2". 
Thus 


V(G) = flog, (|E] + 1). 











So that G is set-semigraceful. 





Your time is limited, so don’t waste it living someone else’s life. 


By Steve Jobs. 


First International Conference 


On Smarandache Multispace and Multistructure 


Organized by Dr.Linfan Mao, Academy of Mathematics and Systems, Chinese Academy of 
Sciences, Beijing 100190, P.R.China. In American Mathematical Society’s Calendar website: 


http://www.ams.org/meetings/calendar/2013_jun28-30_beiging100190.html 


June 28-30, 2013, Send papers by June 1, 2013 to Dr.Linfan Mao by regular mail to the above 
postal address, or by email to maolinfan@163.com. 


A Smarandache multispace (or S-multispace) with its multistructure is a finite or infinite 
(countable or uncountable) union of many spaces that have various structures. The spaces 
may overlap, which were introduced by Smarandache in 1969 under his idea of hybrid science: 
combining different fields into a unifying field, which is closer to our real life world since we live 
in a heterogeneous space. Today, this idea is widely accepted by the world of sciences. 

The S-multispace is a qualitative notion, since it is too large and includes both metric 
and non-metric spaces. It is believed that the smarandache multispace with its multistructure 
is the best candidate for 21st century Theory of Everything in any domain. It unifies many 
knowledge fields. A such multispace can be used for example in physics for the Unified Field 
Theory that tries to unite the gravitational, electromagnetic, weak and strong interactions. 
Or in the parallel quantum computing and in the mu-bit theory, in multi-entangled states 
or particles and up to multi-entangles objects. We also mention: the algebraic multispaces 
(multi-groups, multi-rings, multi-vector spaces, multi-operation systems and multi-manifolds, 
geometric multispaces (combinations of Euclidean and non-Euclidean geometries into one space 
as in Smarandache geometries), theoretical physics, including the relativity theory, the M-theory 
and the cosmology, then multi-space models for p-branes and cosmology, etc. 

The multispace and multistructure were first used in the Smarandache geometries (1969), 
which are combinations of different geometric spaces such that at least one geometric axiom 
behaves differently in each such space. In paradoxism (1980), which is a vanguard in literature, 
arts, and science, based on finding common things to opposite ideas, i.e. combination of con- 
tradictory fields. In neutrosophy (1995), which is a generalization of dialectics in philosophy, 
and takes into consideration not only an entity < A > and its opposite < AntiA > as dialec- 
tics does, but also the neutralities jneutA; in between. Neutrosophy combines all these three 
<A>, < AntiA >, and < neutA > together. Neutrosophy is a metaphilosophy, including 
neutrosophic logic, neutrosophic set and neutrosophic probability (1995), which have, behind the 
classical values of truth and falsehood, a third component called indeterminacy (or neutrality, 
which is neither true nor false, or is both true and false simultaneously - again a combination 
of opposites: true and false in indeterminacy). Also used in Smarandache algebraic structures 
(1998), where some algebraic structures are included in other algebraic structures. 


All reviewed papers submitted to this conference will appear in itsProceedings, published 
in USA this year. 
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